Search NASA⌕ Search

SEARCH · Search NASA

Results for “Level set method”

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,783 records · Page 99

Strict Constraint Feasibility in Analysis and Design of Uncertain Systems

This paper proposes a methodology for the analysis and design optimization of models subject to parametric uncertainty, where hard inequality constraints are present. Hard constraints are those that must be satisfied for all parameter realizations prescribed by the uncertainty model. Emphasis is given to uncertainty models prescribed by norm-bounded perturbations from a nominal parameter value, i.e., hyper-spheres, and by sets of independently bounded uncertain variables, i.e., hyper-rectangles. These models make it possible to consider sets of parameters having comparable as well as dissimilar levels of uncertainty. Two alternative formulations for hyper-rectangular sets are proposed, one based on a transformation of variables and another based on an infinity norm approach. The suite of tools developed enable us to determine if the satisfaction of hard constraints is feasible by identifying critical combinations of uncertain parameters. Since this practice is performed without sampling or partitioning the parameter space, the resulting assessments of robustness are analytically verifiable. Strategies that enable the comparison of the robustness of competing design alternatives, the approximation of the robust design space, and the systematic search for designs with improved robustness characteristics are also proposed. Since the problem formulation is generic and the solution methods only require standard optimization algorithms for their implementation, the tools developed are applicable to a broad range of problems in several disciplines.

Crespo, Luis G.↗

Simple Cell Balance Circuit

A method has been developed for continuous cell voltage balancing for rechargeable batteries (e.g. lithium ion batteries). A resistor divider chain is provided that generates a set of voltages representing the ideal cell voltage (the voltage of each cell should be as if the cells were perfectly balanced). An operational amplifier circuit with an added current buffer stage generates the ideal voltage with a very high degree of accuracy, using the concept of negative feedback. The ideal voltages are each connected to the corresponding cell through a current- limiting resistance. Over time, having the cell connected to the ideal voltage provides a balancing current that moves the cell voltage very close to that ideal level. In effect, it adjusts the current of each cell during charging, discharging, and standby periods to force the cell voltages to be equal to the ideal voltages generated by the resistor divider. The device also includes solid-state switches that disconnect the circuit from the battery so that it will not discharge the battery during storage. This solution requires relatively few parts and is, therefore, of lower cost and of increased reliability due to the fewer failure modes. Additionally, this design uses very little power. A preliminary model predicts a power usage of 0.18 W for an 8-cell battery. This approach is applicable to a wide range of battery capacities and voltages.

Johnson, Steven D.↗

Development of a Trip Energy Estimation Model Using Real-World Global Positioning System Driving Data: Preprint

A data-driven technique for estimation of energy requirements for a proposed vehicle trip has been developed. Based on over 700,000 miles of driving data, the technique has been applied to generate a model that estimates trip energy requirements. The model uses a novel binning approach to categorize driving by road type, traffic conditions, and driving profile. The trip-level energy estimations can easily be aggregated to any higher-level transportation system network desired. The model has been tested and validated on the Austin, Texas, data set used to build this model. Ground-truth energy consumption for the data set was obtained from Future Automotive Systems Technology Simulator (FASTSim) vehicle simulation results. The energy estimation model has demonstrated 12.1 percent normalized total absolute error. The energy estimation from the model can be used to inform control strategies in routing tools, such as change in departure time, alternate routing, and alternate destinations, to reduce energy consumption. The model can also be used to determine more accurate energy consumption of regional or national transportation networks if trip origin and destinations are known. Additionally, this method allows the estimation tool to be tuned to a specific driver or vehicle type.

33 ADVANCED PROPULSION SYSTEMS↗

Unbalanced Nested Random Effects Estimation of Variance Components

To investigate the contributing factors of variance in the measurement of an iodine 127 (127I) sample, we implement a nested random effects analysis of variance (ANOVA). Historically, the reported uncertainty on a measurement of 127I has been obtained by methods of forward uncertainty propagation because there is typically only one replicate of a given sample for which to estimate the uncertainty. When samples are processed there are several types of quality control (QC) standards analyzed along-side the unknown samples with two to five replicates for each. Assuming the variance observed in these replicate QC standards is representative of that of the unknown samples, we use these data in a nested random effects ANOVA to estimate the total uncertainty of a measurement. We demonstrate this approach with two sets of measurements from Idaho National Laboratory and compare the results with the forward uncertainty propagation approach. Variance component estimates for the coarsest level of the nesting structure were most imprecise because replicates were most limited at these levels. We find that the results largely agree between forward propagation and ANOVA, and the greatest contributors of variance are due to instrument variation and chemical processing, with human processing being among the smallest contributors. This analysis provides reassurance that the reported uncertainties using forward propagation are reasonable and the process is well controlled. We propose a future designed experiment to increase replicates at the coarsest level of the hierarchy to improve estimates of these variance components.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scalable line and plane relaxation in a parallel structured multigrid solver

The efficient solution of sparse, linear systems that arise through the discretization of partial differential equations remains a key challenge for a range of high performance scientific simulations. One approach for reducing data movement and improving performance is by exposing and exploiting structure in a problem through the use of robust structured multilevel solvers. By choosing coarsening that preserves the structure of the problem, these methods maintain efficient structured computation and communication throughout the multigrid hierarchy. However, when coarsening is not permitted to be dependent on the operator, anisotropy must be addressed by the smoother — producing error compatible for coarse-grid correction with structured coarsening. Here, the components required in a scalable parallel structured solver are described with a focus on memory and communication efficiency of robust smoothers. While the implementation of communication and memory reduction techniques in smoothers integrated in a complete 3D solver present a significant engineering challenge, a novel approach is proposed that addresses these challenges systematically through a change to the solver’s execution model. Enabled by user-level threading paired with a set of data and communication abstractions, this approach permits seamless aggregation of communication in plane smoothers — directly reusing code for a 2D distributed multilevel cycle. Results show an effective reduction in communication costs for coarse-grid problems, and result in a speedup of 8.7x in smoothing routines shown in Fig. 12 using this approach. This produces a significant improvement to strong scalability while maintaining favorable weak scaling behavior. Finally, a parallel scaling study using a series of refined meshes is included that demonstrates the effectiveness of this approach in an application of interest.

97 MATHEMATICS AND COMPUTING↗

A Technology Roadmap for Advanced Geothermal Well Construction

Recent analysis has shown that by 2050, geothermal energy can supply up to 90 GW of power to the U.S. electrical grid. Reducing the lifecycle cost of geothermal wells is critical to achieving this level of adoption. While time-related costs of drilling remain an important issue, there is also a need to focus on the costs associated with casing and cementing of geothermal wells. This report provides a roadmap of R&D activities needed to reduce the life cycle cost of casing and cementing while improving life-of-well performance. The purpose of the roadmap is to develop priorities for R&D in the areas of (1) high-performance and cost-effective materials for target geothermal well conditions, (2) well construction methods and techniques to reduce construction costs, and (3) methods and techniques to decrease to decrease long-term operating costs. The roadmap sets targets for the next 10-year period and purposely excludes R&D efforts into improved rock reduction. The roadmap was developed by a joint working group of experts from Oak Ridge National Laboratory, Brookhaven National Laboratory, and Sandia National Laboratories as well as 35 experts from industry and academia. Input from a series of information gathering sessions was used to inform and guide the development of the roadmap.

15 GEOTHERMAL ENERGY↗

Evapotranspiration as a Regional Climate Priority: Results from a NASA/USDA Workshop

On April 5 to 7, 2011, the National Aeronautics and Space Administration (NASA) and the United States Department of Agriculture-Agricultural Research Service (USDA-ARS) sponsored a Workshop on Evapotranspiration (ET) in Silver Spring Maryland. The workshop was a response to a recommendation in the 2009-2011 GEO (Group on Earth Observations) Work Plan that a workshop on ET should be held to discuss issues related to ET products and services and the potential for incorporating ET activities into the 2012-2015 GEO Work Plan. The workshop had a regional emphasis, although there were several excellent international and global presentations including one on the GEWEX LANDFLUX project. The different scales of these activities suggests that a framework is needed that can accommodate both regional and global ET activities. Despite limitations with the workshop's scheduling, it attracted 76 experts who contributed informative presentations and insightful discussions. The goals of the workshop involved the exchange of information and ideas and the development of plans for providing more visibility for ET issues. Specific objectives included 1) defining the needs and requirements for evapotranspiration data in weather and climate studies, in natural and agro-ecoystem monitoring, and in water resource management; 2) reviewing the methods used to measure and model evapotranspiration; 3) assessing surface and satellite observation systems required to support ET measurement, modeling and evaluation; 4) assessing the feasibility of developing a proposal for a task on evapotranspiration for the 2012-2015 GEO Work Plan, and 5) exploring the level of support and consensus for developing a strategy for establishing evapotranspiration as an Essential Climate Variable (ECV) within the Global Climate Observing System (GCOS) framework The workshop featured a combination of oral presentations and breakout group sessions focused on the above objectives. There were also poster presentations providing opportunities for one-on-one discussions of ET modeling and measurement techniques. Presentations by users of ET data set the tone for the workshop. In the USA at the national and regional levels water rights issues represent a major opportunity for ET applications. ET data play a major role in estimating water loss due to irrigation, the largest cause of consumptive water loss in the USA, particularly in the West. Irrigation requirements are relatively specific since the needs are clearly defined by the geometry and number of the irrigation systems and can be monitored with high resolution satellite data. There was a strong consensus that land surface temperature (LST) at high resolution is critical for monitoring irrigation. State governments have made commitments to more efficient water management in the western US, but they need full access to improved and more timely ET data and applications to implement this plan. Water managers also reported that in spite of the recent development of. new techniques, the procedures used in some of the water balance calculations in some states are out of date and do not take advantage of new observational and data assimilation systems. The development of ET forecasts for water management is also seen as a priority. Although ET forecasts are currently being produced on an experimental basis these predictions could be improved by considering ET as a dynamic prediction variable in models and by increasing the time resolution of these ET predictions.

Lawford, Richard↗

Initiating Molecular Growth in the Interstellar Medium via Dimeric Complexes of Observed Ions and Molecules

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

Bera, Partha P.↗

Probabilistic Analysis for Comparing Fatigue Data Based on Johnson-Weibull Parameters

Leonard Johnson published a methodology for establishing the confidence that two populations of data are different. Johnson's methodology is dependent on limited combinations of test parameters (Weibull slope, mean life ratio, and degrees of freedom) and a set of complex mathematical equations. In this report, a simplified algebraic equation for confidence numbers is derived based on the original work of Johnson. The confidence numbers calculated with this equation are compared to those obtained graphically by Johnson. Using the ratios of mean life, the resultant values of confidence numbers at the 99 percent level deviate less than 1 percent from those of Johnson. At a 90 percent confidence level, the calculated values differ between +2 and 4 percent. The simplified equation is used to rank the experimental lives of three aluminum alloys (AL 2024, AL 6061, and AL 7075), each tested at three stress levels in rotating beam fatigue, analyzed using the Johnson- Weibull method, and compared to the ASTM Standard (E739 91) method of comparison. The ASTM Standard did not statistically distinguish between AL 6061 and AL 7075. However, it is possible to rank the fatigue lives of different materials with a reasonable degree of statistical certainty based on combined confidence numbers using the Johnson- Weibull analysis. AL 2024 was found to have the longest fatigue life, followed by AL 7075, and then AL 6061. The ASTM Standard and the Johnson-Weibull analysis result in the same stress-life exponent p for each of the three aluminum alloys at the median, or L(sub 50), lives

Vlcek, Brian L.↗

Carbon Dioxide: Surprising Effects on Decision Making and Neurocognitive Performance

The occupants of modern submarines and the International Space Station (ISS) have much in common as far as their air quality is concerned. Air is polluted by materials offgassing, use of utility compounds, leaks of systems chemicals, and anthropogenic sources. The primary anthropogenic compound of concern to submariners and astronauts has been carbon dioxide (CO2). NASA and the US Navy rely on the National Research Council Committee on Toxicology (NRC-COT) to help formulate exposure levels to CO2 that are thought to be safe for exposures of 3-6 months. NASA calls its limits Spacecraft Maximum Allowable Concentrations (SMACs). Years of experience aboard the ISS and a recent publication on deficits in decision making in ground-based subjects exposed briefly to 0.25% CO2 suggest that exposure levels that have been presumed acceptable to preserve health and performance need to be reevaluated. The current CO2 exposure limits for 3-6 months set by NASA and the UK Navy are 0.7%, and the limit for US submariners is 0.5%, although the NRC-COT recommended a 90-day level of 0.8% as safe a few years ago. NASA has set a 1000-day SMAC at 0.5% for exploration-class missions. Anecdotal experience with ISS operations approaching the current 180-day SMAC of 0.7% suggest that this limit is too high. Temporarily, NASA has limited exposures to 0.5% until further peer-reviewed data become available. In the meantime, a study published last year in the journal Environmental Health Perspectives (Satish U, et al. 2012) demonstrated that complexdecision- making performance is somewhat affected at 0.1% CO2 and becomes "dysfunctional" for at least half of the 9 indices of performance at concentrations approaching 0.25% CO2. The investigators used the Strategic Management Simulation (SMS) method of testing for decisionmaking ability, and the results were so surprising to the investigators that they declared that their findings need to be independently confirmed. NASA has responded to the findings on decision making by developing a study to either confirm or refute the published results. In addition, other neurocognitive tests that have been and will be used aboard the ISS will be part of the ground-based study. Further evaluations during the exposures will include ocular and cardiovascular effects of CO2 exposure. In addition, the 1- carbon genetics of the test subjects will be evaluated to determine if any individual subjects are unusually susceptible to CO2 exposure because of genetic factors. Our plan is to extend the published study from Satish by employing a wider series of CO2 exposures to include concentrations of 0.06, 0.12, 0.25, and 0.50%. We believe that our findings will be of keen interest to submariners and to those designing "tight" buildings.

James, John T.↗

A Load Switching Group based Feeder-level Microgrid Energy Management Algorithm for Service Restoration in Power Distribution System

This paper presents a load switching group based energy management system (LSG-EMS) for operating microgrids on a distribution feeder powered by one or multiple grid-forming distributed energy resources. Loads on a distribution feeder are divided into load switching groups that can be remotely switched on and off. The LSG-EMS algorithm, formulated as a mixed-integer linear programming (MILP) problem, has an objective function of maximizing the served loads while minimizing the total number of switching actions. A new set of topology constraints are developed for allowing multiple microgrids to be formed on the feeder and selecting the optimal supply path. Customer comfort is accounted for by maximizing the supply duration in the customer preferred service period and enforcing a minimum service duration. The proposed method is demonstrated on a modified IEEE 33-bus system using actual customer data. Simulation results show that the LSG-EMS successfully coordinates multiple grid-forming sources by selecting an optimal supply topology that maximizes the supply period of both the critical and noncritical loads while minimizing customer service interruptions in the service restoration process.

Hu, Rongxing↗

Disentangling long and short distances in momentum-space TMDs

The extraction of nonperturbative TMD physics is made challenging by prescriptions that shield the Landau pole, which entangle long- and short-distance contributions in momentum space. The use of different prescriptions then makes the comparison of fit results for underlying nonperturbative contributions not meaningful on their own. We propose a model-independent method to restrict momentum-space observables to the perturbative domain. This method is based on a set of integral functionals that act linearly on terms in the conventional position-space operator product expansion (OPE). Artifacts from the truncation of the integral can be systematically pushed to higher powers in Λ QCD /k T . We demonstrate that this method can be used to compute the cumulative integral of TMD PDFs over k T ≤ k$^{cut}_{T}$ in terms of collinear PDFs, accounting for both radiative corrections and evolution effects. This yields a systematic way of correcting the naive picture where the TMD PDF integrates to a collinear PDF, and for unpolarized quark distributions we find that when renormalization scales are chosen near k$^{cut}_{T}$, such corrections are a percent-level effect. We also show that, when supplemented with experimental data and improved perturbative inputs, our integral functionals will enable model-independent limits to be put on the non-perturbative OPE contributions to the Collins-Soper kernel and intrinsic TMD distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An evaluation of point-insertion sequences for incremental Delaunay tessellations

Currently, incremental algorithms may be seen as the lowest-cost computational methods to generate Delaunay tessellations in several point distributions. In this work, eight point-insertion sequences in incremental algorithms for generating Delaunay tessellations are evaluated. More specifically, four point-insertion sequences in incremental algorithms for generating Delaunay tessellations are proposed: with orders given by the red–black tree with in-order and level-order traversals, spiral ordering, and H-indexing. These four incremental algorithms with such sequences are compared with four incremental algorithms with point-insertion orders given by the following sequences: the Hilbert and Lebesgue curves, cut-longest-edge kd-tree, and random order. Six 2-D and seven 3-D point distributions are tested, with sets ranging from 25,000 to 8,000,000 points. The results of computational and storage costs of these eight algorithms are analyzed. It follows that the incremental algorithm with a point-insertion sequence in the order given by the cut-longest-edge kd-tree shows the lowest computational and storage costs of the sequences tested.

97 MATHEMATICS AND COMPUTING↗

Exosomes and Exosomal MicroRNAs in Prostate Cancer Radiation Therapy

Despite current risk stratification systems using traditional clinicopathologic factors, many localized and locally advanced prostate cancers fail radical treatment (ie, radical prostatectomy, radiation therapy with or without androgen deprivation therapy). Therefore, a pressing need exists for enhanced methods of disease stratification through novel prognostic and predictive tools that can reliably be applied in clinical practice. Exosomes are 50- to 150-nm small vesicles released by cancer cells that reflect the genetic and nongenetic materials of parent cancer cells. Cancer cells can contain distinct sets of microRNA profiles, the expression of which can change owing to stress such as radiation therapy. These alterations or distinctions in contents allow exosomes to be used as prognostic and/or predictive biomarkers and to monitor the treatment response. Additionally, microRNAs have been shown to influence multiple processes in prostate tumorigenesis, including cell proliferation, induction of apoptosis, migration, oncogene inhibition, and radioresistance. Thus, comparative exosomal microRNA profiling at different levels could help portray tumor aggressiveness and response to radiation therapy. Although technical challenges persist in exosome isolation and characterization, recent improvements in microRNA profiling have evolved toward in-depth analyses of the exosomal cargo and its functions. We have reviewed the role of exosomes and exosomal microRNAs in biologic processes of prostate cancer progression and radiation therapy response, with a particular focus on the development of clinical assays for treatment personalization.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Limited-angle effect compensation for respiratory binned cardiac SPECT

Purpose: In cardiac single photon emission computed tomography (SPECT), respiratory-binned study is used to combat the motion blur associated with respiratory motion. However, owing to the variability in respiratory patterns during data acquisition, the acquired data counts can vary significantly both among respiratory bins and among projection angles within individual bins. If not properly accounted for, such variation could lead to artifacts similar to limited-angle effect in image reconstruction. In this work, the authors aim to investigate several reconstruction strategies for compensating the limited-angle effect in respiratory binned data for the purpose of reducing the image artifacts. Methods: The authors first consider a model based correction approach, in which the variation in acquisition time is directly incorporated into the imaging model, such that the data statistics are accurately described among both the projection angles and respiratory bins. Afterward, the authors consider an approximation approach, in which the acquired data are rescaled to accommodate the variation in acquisition time among different projection angles while the imaging model is kept unchanged. In addition, the authors also consider the use of a smoothing prior in reconstruction for suppressing the artifacts associated with limited-angle effect. In our evaluation study, the authors first used Monte Carlo simulated imaging with 4D NCAT phantom wherein the ground truth is known for quantitative comparison. The authors evaluated the accuracy of the reconstructed myocardium using a number of metrics, including regional and overall accuracy of the myocardium, uniformity and spatial resolution of the left ventricle (LV) wall, and detectability of perfusion defect using a channelized Hotelling observer. As a preliminary demonstration, the authors also tested the different approaches on five sets of clinical acquisitions. Results: The quantitative evaluation results show that the three compensation methods could all, but to different extents, reduce the reconstruction artifacts over no compensation. In particular, the model based approach reduced the mean-squared-error of the reconstructed myocardium by as much as 40%. Compared to the approach of data rescaling, the model based approach further improved both the overall and regional accuracy of the myocardium; it also further improved the lesion detectability and the uniformity of the LV wall. When ML reconstruction was used, the model based approach was notably more effective for improving the LV wall; when MAP reconstruction was used, the smoothing prior could reduce the noise level and artifacts with little or no increase in bias, but at the cost of a slight resolution loss of the LV wall. The improvements in image quality by the different compensation methods were also observed in the clinical acquisitions. Conclusions: Compensating for the uneven distribution of acquisition time among both projection angles and respiratory bins can effectively reduce the limited-angle artifacts in respiratory-binned cardiac SPECT reconstruction. Direct incorporation of the time variation into the imaging model together with a smoothing prior in reconstruction can lead to the most improvement in the accuracy of the reconstructed myocardium.

60 APPLIED LIFE SCIENCES↗

Progress Update on Segmentation of Neutron Tomographic Images, 2018

The present report documents the techniques developed to segment fast-neutron tomographic images of objects consisting of assemblies of machined parts with distinct boundaries. This report satisfies the fiscal year 2018 technical deliverable for project OR16-3DTomography-PD3Jb, “3D Tomography and Image Processing Using Fast Neutrons,” to report on the extension of geometric shape-finding algorithms to three dimensions. The project has two overall goals. The first of these goals is to extend associated-particle fast-neutron transmission and, particularly, induced-reaction tomographic imaging algorithms to three dimensions. This aspect of the project is beyond the scope of this report. The second goal is to automatically segment the resultant tomographic images into constituent parts and then extract information about the parts, such as the class of shape and potentially the shape parameters. This report describes progress for the image segmentation part of the project. Imaging techniques have been developed as high-confidence methods for confirming the presence and configuration of special nuclear materials (SNM). This high confidence is achieved at the cost of revealing considerable information that may be undesirable to share with the operator of the equipment. One potential way to make use of the high-confidence of imaging methods without revealing imaging data is to employ automated analysis that can extract meaningful attributes of the SNM without showing imaging data to the operator. An essential step of this automated analysis is the segmentation of the image into its constituent parts. For example, a three-dimensional (3D) fast-neutron tomographic image of an assembly containing uranium could be reconstructed and then segmented into discrete parts. The properties of a constituent part identified as uranium would be of interest and could be extracted. For instance, the volume and density of the part would together determine its mass. Moreover, its shape might be expected to fall into a particular geometric class, such as a cube, cylinder, or sphere, and the parameters of the shape (such as side length, height, or diameter) would be expected to fall in a given range. In such an analysis, automated identification of the boundaries and properties of the uranium part is an essential step. A key goal of the present work is to infer the boundaries of objects in fast-neutron tomographic images via the general method of image segmentation. Development of this capability can be broken into more manageable steps, including 1. operator-guided segmentation of two-dimensional (2D) fast-neutron tomographic images, 2. operator-guided segmentation of 3D fast-neutron tomographic images, 3. development of automated segmentation algorithms, and 4. extraction of shape parameters from constituent volumes in 3D. In the first year, the project team reported progress on the first step. Last year, progress on the second step, the extension into 3D, was reported with application to simulated data. The present report describes additional progress on the second step, and progress toward the third step, automated segmentation algorithms. In particular, a multiphase level set approach for image segmentation in three dimensions has been applied to measured 3D data for the first time. Activities related to the fourth step will also be reported, but it will focus on estimation of 3D volumes rather than shape parameters. In this report, the steps before image segmentation will be described, and the multiphase level sets approach to image segmentation will be briefly reviewed. This will be followed by a discussion of the issues encountered when applying the level set technique to measured 3D transmission images, and then a description of the inclusion of induced-reaction data will be provided. Finally, segmentation results for various measured target assemblies will be presented and the image segmentation portion of this project will be summarized and future directions will be discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CarbStor: Development, Analysis and Modification of Carbon Storing Model Soil Communities

Soil microbial communities carry out a number of key processes including plant growth promotion, bioremediation and cycling of nutrients. Carbon cycling is among the most important of these nutrients that are metabolized and processed by the soil microbial community. Many of the carbon inputs are converted to alternative organic forms of carbon that can be used by plants or act as biomass for microbial growth. However, inorganic forms of carbon can also be produced by soil microbial communities including calcium carbonate (CaCO 3 ). Production of calcium carbonate is beneficial for the ecosystem in several ways: it can stabilize soils and improve soil health, especially denser soils with high clay content, it can act as a method of bioremediation, it can serve as an alternative carbon source for plants and it can be a way to store carbon in soil in a stable, inorganic manner for the long term. While the chemistry surrounding individual species carrying out this process is well known what is lacking is an understanding of how species interact in a community to drive carbonate production. As all microbial species in soil exist in a community setting gaining this knowledge is critical to our predicting and controlling this microbial phenotype to greatly improve soil health. The CarbStor project is focused on developing, analyzing and modifying defined microbial soil consortia that express phenotypes at both the species and community level to convert carbon into recalcitrant stable sources such as precipitated carbonate or microbial necromass. To take full advantage of the soil community for this process we will need to fill several key knowledge gaps (KG), three of which are the focus of CarbStor. KG1: Whether and to what degree microbial communities can be developed that produce precipitated carbon via microbial metabolism. KG2: What interspecies interactions drive the individual member phenotypes in defined communities that lead to carbon precipitation. KG3: How can these interactions be modified to enhance carbon sequestration beyond what native communities are capable of. We hypothesize that in a carbon sequestering community only a subset of species will express phenotypes related to carbon storage processes. We also hypothesize that these phenotypes are expressed as a result of interactions with other species in the community that are not involved in carbon storage processes and that these interactions can be harnessed to enhance community carbon sequestration.

54 ENVIRONMENTAL SCIENCES↗

Spectral Retrieval of Latent Heating Profiles from TRMM PR data: Moistening Estimates over Tropical Ocean Regions - Part 3

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of precipitation formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the tropics with the associated latent heating (LH) accounting for threefourths of the total heat energy available to the Earth's atmosphere. In the last decade, it has been established that standard products of LH from satellite measurements, particularly TRMM measurements, would be a valuable resource for scientific research and applications. Such products would enable new insights and investigations concerning the complexities of convection system life cycles, the diabatic heating controls and feedbacks related to rne-sosynoptic circulations and their forecasting, the relationship of tropical patterns of LH to the global circulation and climate, and strategies for improving cloud parameterizations In environmental prediction models. However, the LH and water vapor profile or budget (called the apparent moisture sink, or Q2) is closely related. This paper presented the development of an algorithm for retrieving Q2 using 'TRMM precipitation radar. Since there is no direct measurement of LH and Q2, the validation of algorithm usually applies a method called consistency check. Consistency checking involving Cloud Resolving Model (CRM)-generated LH and 42 profiles and algorithm-reconstructed is a useful step in evaluating the performance of a given algorithm. In this process, the CRM simulation of a time-dependent precipitation process (multiple-day time series) is used to obtain the required input parameters for a given algorithm. The algorithm is then used to "~econsti-LKt~h"e heating and moisture profiles that the CRM simulation originally produced, and finally both sets of conformal estimates (model and algorithm) are compared each other. The results indicate that discrepancies between the reconstructed and CM-simulated profiles for Q2, especially at low levels, are larger than those for latent heat. Larger discrepancies in Q2 at low levels are due to moistening for non-precipitating region that algorithm cannot reconstruct. Nevertheless, the algorithm-reconstructed total Q2 profiles are in good agreement with the CRM-simulated ones.

Shige, S.↗