Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERCEPTION”

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 523 records · Page 29

Leaf angle as a leaf and canopy trait: Rejuvenating its role in ecology with new technology

Abstract Life on Earth depends on the conversion of solar energy to chemical energy by plants through photosynthesis. A fundamental challenge in optimizing photosynthesis is to adjust leaf angles to efficiently use the intercepted sunlight under the constraints of heat stress, water loss and competition. Despite the importance of leaf angle, until recently, we have lacked data and frameworks to describe and predict leaf angle dynamics and their impacts on leaves to the globe. We review the role of leaf angle in studies of ecophysiology, ecosystem ecology and earth system science, and highlight the essential yet understudied role of leaf angle as an ecological strategy to regulate plant carbon–water–energy nexus and to bridge leaf, canopy and earth system processes. Using two models, we show that leaf angle variations have significant impacts on not only canopy‐scale photosynthesis, energy balance and water use efficiency but also light competition within the forest canopy. New techniques to measure leaf angles are emerging, opening opportunities to understand the rarely‐measured intraspecific, interspecific, seasonal and interannual variations of leaf angles and their implications to plant biology and earth system science. We conclude by proposing three directions for future research.

Yang, Xi↗

Constraints to nitrogen acquisition of terrestrial plants under elevated CO 2

Abstract A key part of the uncertainty in terrestrial feedbacks on climate change is related to how and to what extent nitrogen (N) availability constrains the stimulation of terrestrial productivity by elevated CO 2 ( eCO 2 ), and whether or not this constraint will become stronger over time. We explored the ecosystem‐scale relationship between responses of plant productivity and N acquisition to eCO 2 in free‐air CO 2 enrichment ( FACE ) experiments in grassland, cropland and forest ecosystems and found that: (i) in all three ecosystem types, this relationship was positive, linear and strong ( r 2 = 0.68), but exhibited a negative intercept such that plant N acquisition was decreased by 10% when eCO 2 caused neutral or modest changes in productivity. As the ecosystems were markedly N limited, plants with minimal productivity responses to eCO 2 likely acquired less N than ambient CO 2 ‐grown counterparts because access was decreased, and not because demand was lower. (ii) Plant N concentration was lower under eCO 2 , and this decrease was independent of the presence or magnitude of eCO 2 ‐induced productivity enhancement, refuting the long‐held hypothesis that this effect results from growth dilution. (iii) Effects of eCO 2 on productivity and N acquisition did not diminish over time, while the typical eCO 2 ‐induced decrease in plant N concentration did. Our results suggest that, at the decennial timescale covered by FACE studies, N limitation of eCO 2 ‐induced terrestrial productivity enhancement is associated with negative effects of eCO 2 on plant N acquisition rather than with growth dilution of plant N or processes leading to progressive N limitation.

Feng, Zhaozhong↗

Biomass increases attributed to both faster tree growth and altered allometric relationships under long‐term carbon dioxide enrichment at a temperate forest

Abstract Increases in atmospheric carbon dioxide (CO 2 ) concentrations are expected to lead to increases in the rate of tree biomass accumulation, at least temporarily. On the one hand, trees may simply grow faster under higher CO 2 concentrations, preserving the allometric relations that prevailed under lower CO 2 concentrations. Alternatively, the allometric relations themselves may change. In this study, the effects of elevated CO 2 (eCO 2 ) on tree biomass and allometric relations were jointly assessed. Over 100 trees, grown at Duke Forest, NC, USA, were harvested from eight plots. Half of the plots had been subjected to CO 2 enrichment from 1996 to 2010. Several subplots had also been subjected to nitrogen fertilization from 2005 to 2010. Allometric equations were developed to predict tree height, stem volume, and aboveground biomass components for loblolly pine ( Pinus taeda L.), the dominant tree species, and broad‐leaved species. Using the same diameter‐based allometric equations for biomass, it was estimated that plots with eCO 2 contained 21% more aboveground biomass, consistent with previous studies. However, eCO 2 significantly affected allometry, and these changes had an additional effect on biomass. In particular, P. taeda trees at a given diameter were observed to be taller under eCO 2 than under ambient CO 2 due to changes in both the allometric scaling exponent and intercept. Accounting for allometric change increased the treatment effect of eCO 2 on aboveground biomass from a 21% to a 27% increase. No allometric changes for the nondominant broad‐leaved species were identified, nor were allometric changes associated with nitrogen fertilization. For P. taeda , it is concluded that eCO 2 affects allometries, and that knowledge of allometry changes is necessary to accurately compute biomass under eCO 2 . Further observations are needed to determine whether this assessment holds for other taxa.

Kim, Dohyoung↗

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗

Radiation capture and conversion efficiencies of Miscanthus sacchariflorus , M. sinensis and their naturally occurring hybrid M . × giganteus

Abstract Miscanthus is a rhizomatous C4 grass of great interest as a biofuel crop because it has the potential to produce high yields over a wide geographical area with low agricultural inputs on marginal land less suitable for food production. At the moment, a clonal interspecific hybrid Miscanthus × giganteus is the most widely cultivated and studied in Europe and the United States, but breeding programmes are developing newer more productive varieties. Here, we quantified the physiological processes relating to whole season yield in a replicated plot trial in Wales, UK . Light capture and conversion efficiency were parameterized for four carefully selected genotypes ( M. sinensis , M. sacchariflorus and Miscanthus × giganteus ). Differences in the canopy architecture in mature stands as measured by the extinction coefficient ( k ) were small (0.55–0.65). Sensitivity analysis on a mathematical model of Miscanthus was performed to quantify the accumulative intercepted photosynthetically active radiation ( iPAR ) in the growing season using (i) k, (ii) variation in the thermal responses of leaf expansion rate, (iii) base temperature for degree days and (iv) date start of canopy expansion. A 10% increase in k or leaf area per degree day both had a minimal effect on iPAR (3%). Decreasing base temperature from 10 to 9 °C gave an 8% increase in iPAR . If the starting date for canopy expansion was the same as shoot emergence date, then the iPAR increases by 12.5%. In M . × giganteus, the whole season above ground and total (including below ground) radiation‐use efficiency ( RUE ) ranged from 45% to 37% higher than the noninterspecific hybrid genotypes. The greater yields in the interspecific hybrid M. × giganteus are explained by the higher RUE and not by differences in iPAR or partitioning effects. Studying the mechanisms underlying this complex trait could have wide benefits for both fuel and food production.

Davey, Christopher Lyndon↗

Quantifying the effects of varietal types × management on the spatial variability of sorghum biomass across US environments

Abstract Regional‐scale estimations of sorghum biomass production allow identification of optimum genotype×environment×management (G×E×M) combinations for bioenergy generation. The objective of this study was to determine the degree of contributions of G, E, and M toward variability in sorghum biomass in the United States. Using the Agricultural Production Systems sIMulator in a grid computing platform, biomass was simulated for irrigated and rainfed conditions for 30 years across the United States for four sorghum varietal types (grain—GS, sudangrass—SS, photosensitive—PS, and photo‐insensitive—PI). Simulated biomass was assessed by environments clustered using the sum of intercepted solar radiation ( ir ), mean of temperature stress factor ( tp ) and water stress factor ( sw ). Simulated biomass ranged from 5.8 t ha −1 (GS‐rainfed) to 27.5 t ha −1 (PI‐irrigated). Under high‐temperature environments (mean annual temperature = 25°C), rainfed biomass between 40 and 80 days after planting (DAP) was strongly correlated with sw ( r = 0.64–0.86) and irrigated biomass with ir ( r = 0.68–0.81). Under low‐temperature environments (mean annual temperature = 18°C) after 40 DAP, tp and ir had greater effects than sw ( r = 0.55–0.82). Biomass variance was mainly explained by varietal type (50%–76%) in all environments×irrigation combinations, except in the high‐ and mid‐temperature environments under rainfed conditions where rainfall had the major effect (25%–45%). However, when mean temperature during the growing season decreased from 25°C (high environments) to 18°C (low environments), the contribution of mean temperature to biomass variance increased from 7% to 34% (rainfed) and from 4% to 36% (irrigated). Varietal type had the larger interactions with other factors independently of the environment and irrigation. We demonstrated a need to quantify (i) the main G×E×M drivers of biomass variability based on environmental stress factors and (ii) the variance contribution of these drivers on sorghum biomass. Our regional‐scale estimations are key inputs for future robust biomass projections of energy sorghum genotypes integrating G×E×M under climate change scenarios.

Ojeda, Jonathan J.↗

Seasonal to Decadal Variability in Focused Groundwater and Contaminant Discharge along a Channelized Stream

Abstract Focused groundwater discharge to streams is problematic at contaminated sites because high fluxes can limit natural attenuation in the hyporheic zone. However, information on location, spatial evolution, and temporal persistence of springs in unlithified sediments over multiyear time scales is limited. We examine discharge at point (~1‐m) to reach (~300‐m) scales along a stream that intercepts trichloroethene and technetium‐99 plumes from a Superfund site. During 2011 to 2012, we seasonally monitored stream and spring flow and contaminant concentrations, along with probing streambed temperatures on a grid in winter and summer, building on prior monitoring during 1999 to 2002. Baseflow measured by both gauging and dye dilution generally increased with distance downstream, and stream and spring discharge varied seasonally, from minima in October to January to maxima in February to June. Thermal anomalies identified by probing occupied approximately 3% to 6% of the reach and typically coincided with visible springs or seeps. Locations of anomalies were similar to those identified in summer 2002, although some orifices disappeared and others emerged. Vertical groundwater fluxes calculated from probing tended to be less than net fluxes calculated from stream discharge, perhaps in part because the assumption of one‐dimensional, steady‐state flow in calculating point fluxes was simplistic. Maximum contaminant concentrations and fluxes decreased between 1999 to 2001 and 2011 to 2012 as a result of partial capture by an upgradient pump‐and‐treat system. Our findings confirm that springs in unlithified sediments can remain stationary within a few meters over decadal time scales, and seasonal variability in discharge can be greater than decadal variability.

Tripathi, Ganesh N.↗

Method for estimating the density of high-level nuclear waste glass

A database of over 1100 silicate glass compositions and densities was compiled and used to evaluate the efficacy of an algorithm for estimating the density of silicate glass compositions. We sought to develop a parsimonious algorithm based on the additivity of partial molar volumes of individual oxide components weighted by their mole fraction in a glass composition. Bound molar volumes were used for oxides in which the density of the oxide bound in a glass matrix was previously determined. The bound molar volumes were known for oxides covering 97.5 mole percent of the database compositional space. The measured glass densities were plotted against the estimated glass densities and a linear regression yielded an R 2 adj. = 0.95 and a slope and intercept of approximately one and zero, respectively. This regression suggests that glass densities estimated by the algorithm, within analysis uncertainty, are equal to the measured densities of the glasses. In addition to the development of the density estimation, we corroborated many of the referenced bound molar volume data used in the parameterization of the estimation algorithm via linear regression of the individual partial molar volumes versus the inverse measured densities (specific volumes) of the glasses in the database.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Spectroscopy outperforms leaf trait relationships for predicting photosynthetic capacity across different forest types

Leaf trait relationships are widely used to predict ecosystem function in Terrestrial Biosphere Models (TBMs), in which leaf maximum carboxylation capacity (V c,max ), an important trait for modelling photosynthesis, can be inferred from other easier- to-measure traits. However, whether trait-V c,max relationships are robust across different forest types remains unclear. Here we used measurements of leaf traits, including one morphological trait (leaf mass per area), three biochemical traits (leaf water content, area-based leaf nitrogen content, and leaf chlorophyll content), one physiological trait (V c,max ), as well as leaf reflectance spectra, and explored their relationships within and across three contrasting forest types in China. We found weak and forest type-specific relationships between V c,max and the four morphological and biochemical traits (R 2 ≤ 0.15), indicated by significantly changing slopes and intercepts across forest types. In contrast, reflectance spectroscopy effectively collapsed the differences in the trait-V c,max relationships across three forest biomes into a single robust model for V c,max (R 2 = 0.77), and also accurately estimated the four traits (R 2 = 0.75-0.94). Furthermore, these findings challenge the traditional use of empirical trait-V c,max relationships in TBMs for estimating terrestrial plant photosynthesis, but also highlight spectroscopy as an efficient alternative for characterizing V c,max and multi-trait variability, with critical insights into ecosystem modeling and functional trait ecology.

54 ENVIRONMENTAL SCIENCES↗

Short-term variation in leaf-level water use efficiency in a tropical forest

Here we report the representation of stomatal regulation of transpiration and CO 2 assimilation is key to forecasting terrestrial ecosystem responses to global change. Given its importance in determining the relationship between forest productivity and climate, accurate and mechanistic model representation of the relationship between stomatal conductance (g s ) and assimilation is crucial. We assess possible physiological and mechanistic controls on the estimation of the g 1 (stomatal slope, inversely proportional to water use efficiency) and g 0 (stomatal intercept) parameters, using diurnal gas exchange surveys and leaf level response curves of six tropical broadleaf evergreen tree species. g 1 estimated from ex-situ response curves averaged 50% less than g 1 estimated from survey data. While g 0 and g 1 varied between leaves of different phenological stages, the trend was not consistent among species. We identified a diurnal trend associated with g 1 and g 0 that significantly improved model projections of diurnal trends in transpiration. The accuracy of modelled g s can be improved by accounting for variation in stomatal behavior across diurnal periods, and between measurement approaches, rather than focusing on phenological variation in stomatal behavior. Additional investigation into the primary mechanisms responsible for diurnal variation in g 1 will be required to account for this phenomena in land surface models.

54 ENVIRONMENTAL SCIENCES↗

Learning an Eddy Viscosity Model Using Shrinkage and Bayesian Calibration: A Jet-in-Crossflow Case Study

In this paper, we demonstrate a statistical procedure for learning a high-order eddy viscosity model (EVM) from experimental data and using it to improve the predictive skill of a Reynolds-averaged Navier–Stokes (RANS) simulator. The method is tested in a three-dimensional (3D), transonic jet-in-crossflow (JIC) configuration. The process starts with a cubic eddy viscosity model (CEVM) developed for incompressible flows. It is fitted to limited experimental JIC data using shrinkage regression. The shrinkage process removes all the terms from the model, except an intercept, a linear term, and a quadratic one involving the square of the vorticity. The shrunk eddy viscosity model is implemented in an RANS simulator and calibrated, using vorticity measurements, to infer three parameters. The calibration is Bayesian and is solved using a Markov chain Monte Carlo (MCMC) method. A 3D probability density distribution for the inferred parameters is constructed, thus quantifying the uncertainty in the estimate. The phenomenal cost of using a 3D flow simulator inside an MCMC loop is mitigated by using surrogate models (“curve-fits”). A support vector machine classifier (SVMC) is used to impose our prior belief regarding parameter values, specifically to exclude nonphysical parameter combinations. The calibrated model is compared, in terms of its predictive skill, to simulations using uncalibrated linear and CEVMs. Finally, we find that the calibrated model, with one quadratic term, is more accurate than the uncalibrated simulator. The model is also checked at a flow condition at which the model was not calibrated.

42 ENGINEERING↗

A Forward Feedback Control Scheme for a Solar Thermochemical Moving Bed Counter-Current Flow Reactor

Abstract Pelletized thermochemical energy storage media has a potential for long-duration energy storage. Production of solid-state energy storage media can be done within a cavity chemical reactor that captures concentrated solar radiation from a solar thermal field. The temperature stability of a solar reactor is directly influenced by the solar flux intercepted. This paper presents a low-order physical model to simulate the dynamic response of temperature inside a tubular plug-flow reactor prototype. Solid granular particles are fed to the reactor from the top whereas a counter-current flowing gas enters the reactor from the bottom. An in-house code was developed to model transient heat transfer of the reactor wall, gas, and moving particles. The model was preliminarily validated with packed beds for different temperature ranges and two gas flowrates. Dynamic response of the reactor temperature is simulated for different input power and gas/particle flowrates. The results show that the system response can be controlled efficiently by utilizing input power (solar flux) as a control parameter. A conventional proportional integral (PI) controller is designed to control the temperature inside the reactor and to maintain it during the solar flux intermittency. The controller parameters are tuned using the Ziegler–Nichols method to ensure optimal system response. The results show that the feedback control model is successful in tracking different reference reactor temperatures within a reasonable settling time of 30 min and eliminated overshoot. This study can be extended to include a hybrid reactor with a multi-input, multi-output variable system.

Energy & Fuels↗

Is Radioactive Tank Waste a Bingham Fluid?

Slurries and sludges across the United States Department of Energy (DOE) complex rank among the most rheologically interesting. Their composition is heterogeneous, spanning a very broad range of particle sizes, densities, and interparticle forces. All exhibit shear thinning, some have yield stresses, and many are thixotropic. Despite the variety, these complex fluids are often represented using the historic Bingham fluid model, which fits higher shear rate data to a simple straight line. The intercept provides a yield stress, which has been a key design parameter in construction of large-scale waste processing facilities. However, many radioactive wastes are simply not Bingham fluids, which extrapolate but poorly across low to intermediate shear rates that are characteristic of typical processing conditions. Indeed, processing shear rates as high as 200 1/s, which has been a typical minimum shear rate used in fitting, are seldom encountered in nuclear waste processing. Therefore, more realistic rheological models are necessary to accurately predict waste processing performance. PNNL recently re-evaluated the rheology of reconstituted Hanford REDOX sludge waste against a wide variety of rheological models including the Bingham, Cross, Cross with yield stress, Carreau, biviscous, Herschel-Bulkley (which includes a power law dependence), Casson, and Gay models. They found that all of the models provided a closer fit than the Bingham model and that the biviscous model and Cross with yield stress model were convincing. However, reconstituted Hanford REDOX sludge waste is but one type of DOE waste and a direct contrast and comparison of these three models against undiluted, unmixed tank waste (actual not simulant) has not been performed previously. Therefore, the purpose of this paper is to evaluate the rheology of actual tank waste with these more accurate rheological models. In this paper, we evaluate select rheological data for slurry samples from Hanford’s AZ-101, AZ-102, and SY-101 waste tanks. In each of these cases, we find that Cross’ model with yield stress and the biviscous model significantly outperform the Bingham fluid model. Furthermore, the AZ-101 data also shows that the shear stress peak that startup significantly exceeds the Bingham yield stress, which is commonly observed in the incipient moments of rheological measurements on simulants. In contrast, Cross’ model may empirically accommodate an initial spike in shear stress at modest shear rates. These are important observations because computational and analytical fluid dynamics simulations rely on rheological constitutive models for accurately and conservatively predict waste processing performance. These findings suggest the need for better rheological modeling of and validation against radioactive waste.

Rheology, Radioactive Waste↗

Implications of Non-Bingham Rheology

One of the key challenges now facing the US Department of Energy (DOE) is the fate of radioactive waste remaining from World War II and the Cold War which is stored underground in tanks some 75 feet in diameter and over 30 feet tall. Over time, the waste has segmented into multiple layers with sludges and slurries at the bottom with salt crust layers often at the top and liquid in between. DOE’s current official baseline plans call for remaining sludges and slurries to be removed from the tanks and converted into a more stable glass waste form. Minimizing worker exposure to radiation drives DOE to use slurry processing techniques to suspend, mobilize, transport, mix, and process the waste. Therefore, a clear and quantitative understanding of Hanford waste rheology is essential for the success of the DOE mission. Historically much of the waste has been characterized using Eugene Bingham’s model more than a century old that provides a straight line fit to higher shear rate data with the intercept suggesting a yield stress and the slope providing the consistency. Yet, Bingham fits overestimate the shear stress at a given shear rate for low to intermediate shear rates, exactly the range of shear rates typically encountered in pipe flow, where shear rates peak along the pipe wall and vanish in the center. This model produces a fictitious yield stress for some of the wastes that do not exhibit yield phenomena. While over estimating the yield stress may be prudent, safe, and conservative for some applications (e.g., pump sizing to ensure that pumps can handle yield stresses), over estimating the rheology may be inaccurate and non-conservative for other applications (e.g., eroding settled particle beds). Therefore, this paper evaluates the slurry rheology of Hanford and Savannah River wastes using a more modern rheological model that fits the full range of experimental data. Although a bias has been recognized and alternative models proposed, the magnitude of this bias and the implications for tank waste have only been qualitatively suggested. The purpose of this paper is to evaluate quantitatively implications of the poor quality of fit between a Bingham model for rheology and experimental data at modest shear rates. We first demonstrate the magnitude of the bias between the data and the Bingham extrapolation. We then evaluate quantitatively the velocity profile under laminar conditions. This analysis shows that the bias may be large (hundreds of percent or more) at modest shear rates and that modest shear rates dominate pipe velocity profiles.

Rheology, Radioactive waste↗

Tunable UV upgrade on OMEGA EP

The OMEGA EP laser has been upgraded to provide a UV wavelength-tunable beam to support the study of wavelength detuning for the mitigation of cross-beam energy transfer in direct-drive inertial confinement fusion. The beamline delivers up to 0.5 TW in pulses up to 1-ns duration (0.1 TW up to 2.5 ns), to either the OMEGA or OMEGA EP target chambers with wavelength tunable from 350.2 to 353.4 nm. The upgrade leverages the existing optical parametric amplification (OPA) system in the short-pulse front end of OMEGA EP Beamline 1 for amplification of a new tunable, narrowband fiber front end over a broad spectral range. The tunable OPA output is spatially shaped to form a round OMEGA-like beam, which is amplified in the OMEGA EP beamline, then frequency tripled and characterized using the existing OMEGA EP long-pulse infrastructure. A new 3ω beam-transport system intercepts the tunable UV beam near the OMEGA EP target chamber and image relays it to the P9 port of the OMEGA target chamber for joint shots with the OMEGA 60-beam laser. Commissioning of the tunable UV capability has been completed, and four experimental campaigns have been supported with the tunable beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative ultrasound characterization of locally advanced breast cancer by estimation of its scatterer properties

Purpose: Tumor grading is an important part of breast cancer diagnosis and currently requires biopsy as its standard. Here, the authors investigate quantitative ultrasound parameters in locally advanced breast cancers that can potentially separate tumors from normal breast tissue and differentiate tumor grades. Methods: Ultrasound images and radiofrequency data from 42 locally advanced breast cancer patients were acquired and analyzed. Parameters related to the linear regression of the power spectrum—midband fit, slope, and 0-MHz-intercept—were determined from breast tumors and normal breast tissues. Mean scatterer spacing was estimated from the spectral autocorrelation, and the effective scatterer diameter and effective acoustic concentration were estimated from the Gaussian form factor. Parametric maps of each quantitative ultrasound parameter were constructed from the gated radiofrequency segments in tumor and normal tissue regions of interest. In addition to the mean values of the parametric maps, higher order statistical features, computed from gray-level co-occurrence matrices were also determined and used for characterization. Finally, linear and quadratic discriminant analyses were performed using combinations of quantitative ultrasound parameters to classify breast tissues. Results: Quantitative ultrasound parameters were found to be statistically different between tumor and normal tissue (p < 0.05). The combination of effective acoustic concentration and mean scatterer spacing could separate tumor from normal tissue with 82% accuracy, while the addition of effective scatterer diameter to the combination did not provide significant improvement (83% accuracy). Furthermore, the two advanced parameters, including effective scatterer diameter and mean scatterer spacing, were found to be statistically differentiating among grade I, II, and III tumors (p = 0.014 for scatterer spacing, p = 0.035 for effective scatterer diameter). The separation of the tumor grades further improved when the textural features of the effective scatterer diameter parametric map were combined with the mean value of the map (p = 0.004). Conclusions: Overall, the binary classification results (tumor versus normal tissue) were more promising than tumor grade assessment. Combinations of advanced parameters can further improve the separation of tumors from normal tissue compared to the use of linear regression parameters. While the linear regression parameters were sufficient for characterizing breast tumors and normal breast tissues, advanced parameters and their textural features were required to better characterize tumor subtypes.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-E-T-90: Concrete Forward-Scatter Fractions for Radiotherapy Shielding Applications

Purpose: There is little instruction within the primary shielding guidance document NCRP 151 for vault designs where the primary beam intercepts the maze. We have conducted a Monte-Carlo study to characterize forward-scattered radiation from concrete barriers with the intent of quantifying what amount of additional shielding outside the primary beam is needed in this situation. Methods: We reproduced our vault in MCNP 5 and simulated spectra obtained from the literature and from our treatment planning system for 10 and 18 MV beams. Neutron and gamma-capture contributions were not simulated. Energy deposited was scored at isocenter in a water phantom, within various cells that comprised the maze, and within cells that comprised the vault door. Tracks were flagged that scattered from within the maze to the door and their contributions were tallied separately. Three different concrete mixtures found in the literature were simulated. An empirically derived analytic equation was used for comparison, utilizing patient scatter fractions to approximate the scatter from concrete. Results: Our simulated data confirms that maze-scattered radiation is a significant contribution to total photon dose at the door. It contributes between 20-35% of the photon shielding workload. Forward-scatter fractions for concrete were somewhat dependent on concrete composition and the relative abundance of higher-Z elements. Scatter fractions were relatively insensitive to changes in the primary photon spectrum. Analytic results were of the same magnitude as simulated results. Conclusions: Forward-scattered radiation from the maze barrier needs to be included in the photon workload for shielding calculations in non-standard vault designs. Scatter fractions will vary with concrete composition, but should be insensitive to spectral changes between machine manufacturers. Further plans for investigation include refined scatter fractions for various concrete compositions, scatter fraction measurement, and simulation of photo-nuclear interaction within the maze barrier for high-energy beams.

07 ISOTOPE AND RADIATION SOURCES↗

SU-E-T-497: Semi-Automated in Vivo Radiochromic Film Dosimetry Using a Novel Image Processing Algorithm

Purpose: To validate an automated image processing algorithm designed to detect the center of radiochromic film used for in vivo film dosimetry against the current gold standard of manual selection. Methods: An image processing algorithm was developed to automatically select the region of interest (ROI) in *.tiff images that contain multiple pieces of radiochromic film (0.5x1.3cm{sup 2}). After a user has linked a calibration file to the processing algorithm and selected a *.tiff file for processing, an ROI is automatically detected for all films by a combination of thresholding and erosion, which removes edges and any additional markings for orientation. Calibration is applied to the mean pixel values from the ROIs and a *.tiff image is output displaying the original image with an overlay of the ROIs and the measured doses. Validation of the algorithm was determined by comparing in vivo dose determined using the current gold standard (manually drawn ROIs) versus automated ROIs for n=420 scanned films. Bland-Altman analysis, paired t-test, and linear regression were performed to demonstrate agreement between the processes. Results: The measured doses ranged from 0.2-886.6cGy. Bland-Altman analysis of the two techniques (automatic minus manual) revealed a bias of -0.28cGy and a 95% confidence interval of (5.5cGy,-6.1cGy). These values demonstrate excellent agreement between the two techniques. Paired t-test results showed no statistical differences between the two techniques, p=0.98. Linear regression with a forced zero intercept demonstrated that Automatic=0.997*Manual, with a Pearson correlation coefficient of 0.999. The minimal differences between the two techniques may be explained by the fact that the hand drawn ROIs were not identical to the automatically selected ones. The average processing time was 6.7seconds in Matlab on an IntelCore2Duo processor. Conclusion: An automated image processing algorithm has been developed and validated, which will help minimize user interaction and processing time of radiochromic film used for in vivo dosimetry.

60 APPLIED LIFE SCIENCES↗