SEARCH · Engineering Papers

Results for “projected gradient 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Reevaluating growing season length controls on net ecosystem production in evergreen conifer forests

Abstract Growing season length (GSL) is a key unifying concept in ecology that can be estimated from eddy covariance-derived estimates of net ecosystem production (NEP). Previous studies disagree on how increasing GSLs may affect NEP in evergreen coniferous forests, potentially due to the variety of methods used to quantify GSL from NEP. We calculated GSL and GSL-NEP regressions at eleven evergreen conifer sites across a broad climatic gradient in western North America using three common approaches: (1) variable length (3–7 days) regressions of day of year versus NEP, (2) a smoothed threshold approach, and (3) the carbon uptake period, followed by a new approach of a method-averaged ensemble. The GSL and the GSL-NEP relationship differed among methods, resulting in linear relationships with variable sign, slope, and statistical significance. For all combinations of sites and methods, the GSL explained between 6% and 82% of NEP with p -values ranging from 0.45 to < 0.01. These results demonstrate the variability among GSL methods and the importance of selecting an appropriate method to accurately project the ecosystem carbon cycling response to longer growing seasons in the future. To encourage this approach in future studies, we outline a series of best practices for GSL method selection depending on research goals and the annual NEP dynamics of the study site(s). These results contribute to understanding growing season dynamics at ecosystem and continental scales and underscore the potential for methodological variability to influence forecasts of the evergreen conifer forest response to climate variability.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Experimental and numerical validation of an ultrasonic visualization system for nuclear reactor application

This paper deals with the design of ultrasonic imaging systems for the next generation of nuclear reactors cooled by liquid metal. Indeed, a generation IV research nuclear reactor is being developed by the Belgian Nuclear Research Center (SCK-CEN) in the frame of the MYRRHA project (for Multipurpose hYbrid Research Reactor for High-tech Applications). This sub-critical/critical reactor is cooled by Lead-Bismuth Eutectic (LBE). The opacity of this liquid metal requires the development of an ultrasonic visualization system for internal inspection and object detection. But due to the peculiar conditions met in the core of the reactor, velocity as well as temperature gradients are expected and are likely to affect directly the ultrasonic propagation. The objective of this work is to validate the ultrasonic imaging strategy by tests performed in a dedicated test rig and by numerical simulations using a ray-tracing method. The experimental investigations have been performed on a specific water facility reproducing conditions similar to those encountered in the core of the MYRRHA reactor. These conditions include the propagation over large distance, and the presence of temperature and velocity gradients. In the MYRRHA reactor application, the distance to be travelled by the acoustic waves of the visualization system is about 5 m, including a reflection. The acoustic absorption, the scattering losses, the beam divergence and the transmitted energy during reflection have been determined as a function of the travelled distance. The experimental values are compared with the literature for validation. The presence of temperature and velocity gradients in the core of the reactor is due to the coolant circulation. These gradients are about 5 K over 0.1 m and 1 m/s over 0.2 m, respectively, and are reproduced in the facility for investigating their influence on the propagation of ultrasounds. The experimental data are used for improving and validating a ray-tracing algorithm. This algorithm is then applied to the full scale application for validating the ultrasonic visualization system for nuclear reactor application. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2014

Stellar kinematics and structural properties of virgo cluster dwarf early-type galaxies from the SMAKCED project. I. Kinematically decoupled cores and implications for infallen groups in clusters

We present evidence for kinematically decoupled cores (KDCs) in two dwarf early-type (dE) galaxies in the Virgo cluster, VCC 1183 and VCC 1453, studied as part of the SMAKCED stellar absorption-line spectroscopy and imaging survey. These KDCs have radii of 1.''8 (0.14 kpc) and 4.''2 (0.33 kpc), respectively. Each of these KDCs is distinct from the main body of its host galaxy in two ways: (1) inverted sense of rotation and (2) younger (and possibly more metal-rich) stellar population. The observed stellar population differences are probably associated with the KDC, although we cannot rule out the possibility of intrinsic radial gradients in the host galaxy. We describe a statistical analysis method to detect, quantify the significance of, and characterize KDCs in long-slit rotation curve data. We apply this method to the two dE galaxies presented in this paper and to five other dEs for which KDCs have been reported in the literature. Among these seven dEs, there are four significant KDC detections, two marginal KDC detections, and one dE with an unusual central kinematic anomaly that may be an asymmetric KDC. The frequency of occurrence of KDCs and their properties provide important constraints on the formation history of their host galaxies. We discuss different formation scenarios for these KDCs in cluster environments and find that dwarf-dwarf wet mergers or gas accretion can explain the properties of these KDCs. Both of these mechanisms require that the progenitor had a close companion with a low relative velocity. This suggests that KDCs were formed in galaxy pairs residing in a poor group environment or in isolation whose subsequent infall into the cluster quenched star formation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2011

Progress Report for Diffusion Welding of the NGNP Process Application Heat Exchangers (Revision 1)

The U.S. Department of Energy selected the high temperature gas-cooled reactor as the basis for the Next Generation Nuclear Plant (NGNP). The NGNP will demonstrate the use of nuclear power for electricity, hydrogen production, and process heat applications. The NGNP Project is currently investigating the use of metallic, diffusion welded, compact heat exchangers to transfer heat from the primary (reactor side) heat transport system to the secondary heat transport system. An intermediate heat exchanger will transfer this heat to downstream applications such as hydrogen production, process heat, and electricity generation. The channeled plates that make up the heat transfer surfaces of the intermediate heat exchanger will have to be assembled into an array by diffusion welding. This report describes the preliminary results of a scoping study that evaluated the diffusion welding process parameters and the resultant mechanical properties of diffusion welded joints using Alloy 800H. The long-term goal of the program is to progress towards demonstration of small heat exchanger unit cells fabricated with diffusion welds. Demonstration through mechanical testing of the unit cells will support American Society of Mechanical Engineers rules and standards development, reduce technical risk, and provide proof of concept for heat exchanger fabrication methods needed to deploy heat exchangers in several potential NGNP configurations.1 Researchers also evaluated the usefulness of modern thermodynamic and diffusion computational tools (Thermo-Calc and Dictra) in optimizing the parameters for diffusion welding of Alloy 800H. The modeling efforts suggested a temperature of 1150°C for 1 hour with an applied pressure of 5 MPa using 15 μm nickel foil as joint filler to reduce chromium oxidation on the welded surfaces. Good agreement between modeled and experimentally determined concentration gradients was achieved

08 HYDROGEN↗

DOE OSTI2021

A new method for measuring the 3D turbulent velocity dispersion of molecular clouds

ABSTRACT The structure and star formation activity of a molecular cloud are fundamentally linked to its internal turbulence. However, accurately measuring the turbulent velocity dispersion is challenging due to projection effects and observational limitations, such as telescope resolution, particularly for clouds that include non-turbulent motions, such as large-scale rotation. Here, we develop a new method to recover the 3D turbulent velocity dispersion (σv,3D) from position–position–velocity (PPV) data. We simulate a rotating, turbulent, collapsing molecular cloud, and compare its intrinsic σv,3D with three different measures of the velocity dispersion accessible in PPV space: (1) the spatial mean of the 2nd-moment map, σi, (2) the standard deviation of the gradient/rotation-corrected 1st-moment map, σ(c − grad), and (3) a combination of (1) and (2), called the ‘gradient-corrected parent velocity dispersion’, $\sigma _{\mathrm{(p}-\mathrm{grad)}}=(\sigma _{\mathrm{i}}^2+\sigma _{(\mathrm{c}-\mathrm{grad)}}^2)^{1/2}$. We show that the gradient correction is crucial in order to recover purely turbulent motions of the cloud, independent of the orientation of the cloud with respect to the line of sight. We find that with a suitable correction factor and appropriate filters applied to the moment maps, all three statistics can be used to recover σv,3D, with method 3 being the most robust and reliable. We determine the correction factor as a function of the telescope beam size for different levels of cloud rotation, and find that for a beam full width at half-maximum f and cloud radius R, the 3D turbulent velocity dispersion can best be recovered from the gradient-corrected parent velocity dispersion via $\sigma _{v,\mathrm{3D}}= \left[(-0.29\pm 0.26)\, f/R + 1.93 \pm 0.15\right] \sigma _{\mathrm{(p}-\mathrm{grad)}}$ for f/R < 1, independent of the level of cloud rotation or LOS orientation.

Stewart, Madeleine↗

DOE OSTI2025

Efficiently resolving the terrestrial-aquatic interface in E3SM with sub-grid methods to improve coastal simulations (Final Technical Report)

The broad goal of this project was to extend the capabilities of MPAS‐Ocean, the ocean component of DOE’s Energy Exascale Earth System Model (E3SM), to better resolve the hydrodynamics of the terrestrial‐aquatic interface without significantly increasing computing resource requirements. Specifically, our goal was to better resolve the small-scale features across the terrestrial‐aquatic interface which includes the low lying coastal floodplain and a dendritic network of hydraulically efficient channels connected to estuarine systems and back bays, barrier islands and high energy inlets. The processes driving flow across the terrestrial‐aquatic interface include upland and floodplain hydrology, tides, wind and atmospheric pressure (including coastally generated storm surge across continental shelves and within estuaries and bays), wind waves and water temperature and salinity gradients.

58 GEOSCIENCES↗

DOE OSTI2013

Energy-Saving Melting and Revert Reduction Technology (E-SMARRT), Task 5.3: Casting Porosity-Free Grain Refined Magnesium Alloys

The objective of this project was to identify the root causes for micro-porosity in magnesium alloy castings and recommend remedies that can be implemented in production. Three casting magnesium alloys were used to determine the causes of porosity that is a common problem in the magnesium casting industry. These alloys were AM50A, AM60A, and AZ91D, and all were grain refined with C 2 Cl 6 . Gas porosity was found to have little effect if proper degassing methods and foundry practices were used. It was found that sound gating practice, directional solidification towards the riser and generous risering should always dominate in order to produce sound castings. In sand cast magnesium, thermal gradients were found to be an important indicator of heat transfer in the mold. A minimum thermal gradient was found for each alloy that would produce sound castings. There were 2.02°C/cm for AM50A, 1.99°C/cm for AM60A, and 1.97°C/cm for AZ91D. If these gradients are met or exceeded when an area of the casting reaches 50°C below the liquidus for the alloy, the casting will be free of shrinkage porosity in that area. The cooling rates measured during casting magnesium alloys in a pre-heated step permanent mold ranged from 12°F/min in the three inches thick step to 200o F/min in the one inch step. The effect of the cooling rate on the grain size and the mechanical properties was noticeable but not substantial. Grain refining with C 2 Cl 6 had a marked effect on the grain size of the permanent mold step casting. It reduced it by 10-15% in AM50A and by a factor of six in AZ91E from 300 µm to around 55 µm. The mechanical properties of both permanent mold cast AZ91E and AM50A were improved by grain refinement. The most prominent increase was in the elongation of AM50A that more than doubled from 8% to 17%. Computer simulation of porosity was conducted for a well gated commercial wheel casting. The simulation did not predict significant porosity in the wheel. The mechanical testing of tensile specimens excised from the wheel indeed confirmed good properties, in line with the specifications for this alloy. Evidence of minimal micro-porosity was detected in sections cut from the wheel and investigated under optical microscopy. This finding confirms the key role played by utilizing optimal gating and risering practices in minimizing porosity in magnesium castings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1993

Measurements of atmospheric carbonyl sulfide during the NASA Chemical Instrumentation Test and Evaluation Project: Implications for the global COS budget

Atmospheric carbonyl sulfide COS concentrations were measured by three analytical systems during the Chemical Instrumentation Test and Evaluation (CITE 3) project. The three systems all used cryogenic sample preconcentration and gas chromatographic (GC) separation but differed in the method of detection. The FPD system used a flame photometric detector, the MS system used a mass selective detector, and the ECD-S system used a fluorinating catalyst followed by an electron capture detector. With the FPD system, we found a mean COS concentration of 510 ppt over the North Atlantic and 442 ppt over the Tropical Atlantic. With the ECD-S system, we found a mean COS concentration of 489 ppt over the North Atlantic and 419 ppt over the Tropical Atlantic. All three systems registered a latitudinal gradient in atmospheric COS of between 1.6 and 2.0 ppt per degree of latitude, with increasing COS concentrations northward which was similar to the gradient measured by Bingemer et al. (1990). It is difficult to reconcile the measured latitudinal concentration gradient with present theories of the global COS budget since the largest sink of COS is thought to be a flux to land plants, most of which are in the northern hemisphere.

Johnson, James E.↗

DOE OSTI2025

Development of Transformative Preparation Methods to Push up High Q&G Performance of FRIB Spare HWR Cryomodule Cavities

The FRIB accelerator project construction, a top priority of US nuclear science, was completed in January 2022, and is now moving to user operation. The stable and reliable operation of the accelerating cryomodules is essential in achieving/fulfilling DOE and user expectations. So far, FRIB cryomodules meet all FRIB specifications for cavity performance. However, during the lifetime of machine operation, degradation of cryomodule performance is possible, as reported in similar operating facilities (CEBAF, SNS). If cryomodule degradation is observed at FRIB, the under-performing cryomodule will require replacement/maintenance. In effort to manage operational reliability, FRIB plans to construct a 0.53 half-wave cryomodule to serve as an active spare. In a parallel effort, FRIB will also work toward increasing operational Q and gradient of spare cryomodule cavities to gain an overall performance margin to support future operational reliability. The current FRIB cavity designs have a potential to operate at gradients higher than 8 MV/m, but are currently limited by field emission (FE) and/or high field Q slope (HFQS); known issue in buffered chemical polished (BCP) treated cavities. The proposal looks to develop transformative surface preparation treatments to improve the operational gradient of spare cryomodules higher than 10 MV/m while maintaining high Q. Thus, increasing operational margin by 30 - 50%. With the goal to improve operational reliability set, the proposal will investigate multiple objectives as possible paths forward to achieve an overall increase in cavity performance and gain a better understanding of SRF limiting mechanisms. The proposal will study the application of different chemical surface treatments to 0.53 half-wave cavities, with the addition of low temperature bakes (LTB), and measure their effects on accelerating performance. Proposed chemical treatments to be explored in this proposal include conventional EP acid mixtures, as well as innovated EP and BCP acid mixtures designed to simplify processing paths in migrating FE and HFQS. The proposed transformative treatment wet N-doping also has the potential to replicate recent advancements in SRF technology relating to nitrogen doping and high Q operation without the requirement for an ultra-high vacuum annealing furnace; currently being developed at FNAL and JLAB. In parallel, high Q performance relating to flux trapping will be investigated with the installation of a second layer of magnetic shielding in the vertical test Dewar. The research objectives presented in the proposal, and their corresponding effects on cavity performance, will provide essential knowledge and future guidance to the SRF community and provide possible paths for future SRF based projects and applications.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Computational diagnostics for flame acceleration and transition to detonation in a hydrogen/air mixture

A new computational diagnostic method for pressure-induced compressibility is proposed by projecting its local contribution to the chemical explosive mode (CEM) in the chemical explosive mode analysis (CEMA) framework. The new method is validated for the study of detonation development during the deflagration-to-detonation transition (DDT) process. The flame characteristics are identified through the quantification of individual CEM contributions of chemical reaction, diffusion, and pressure-induced compressibility. Numerical simulations are performed to investigate the DDT processes in a stoichiometric hydrogen-air mixture. A Godunov algorithm, fifth-order in space, and third-order in time are used to solve the fully compressible Navier-Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model (CDM) described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The new diagnostic method is first applied to onedimensional (1D) canonical flame configurations followed by two-dimensional (2D) simulations of DDT in an obstructed channel where different detonation initiation scenarios are examined using the new CEMA projection formulation. Detailed examinations of the idealized configuration of detonation initiation through shock focusing mechanism at a flame front are also studied using the new formulation. A comparison of the currently proposed CEMA projection and the original formulation by the authors suggests that including the pressure-induced compressibility is essential for the use of CEMA in DDT process. The results also show that the new formulation of CEMA projection can successively capture the detonation initiation through either a gradient mechanism or a direct initiation mechanism, and therefore can be used as an effective local analytical tool for the computational diagnostics of detonation initiation in a DDT process. It was found that detonation development is characterized by a strong contribution of chemistry role to the CEM which is pivotal to the initiation of detonation. The role of compressibility is found enhanced at the edge of the detonation front where diffusion was found to have minimal effects on detonation development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

MO-DE-BRA-06: MrRSCAL: A Radiological Simulation Tool for Resident Education

Purpose: The goal of this project was to create a readily accessible, comprehensive-yet-flexible interactive MRI simulation tool for use in training and education of radiology residents in particular. This tool was developed to take the place of an MR scanner in laboratory activities, as magnet time has become scarce while hospitals are optimizing clinical scheduling for improved throughput. Methods: MrRSCAL (Magnetic resonance Resident Simulation Console for Active Learning) was programmed and coded using Matlab on a Mac workstation utilizing OS X platform. MR-based brain images were obtained from one of the co-authors and processed to generate parametric maps. Scanner sounds are also generated via mp3 convolution of a single MR gradient slew with a time-profile of gradient waveforms. Results: MrRSCAL facilitates the simulation of multiple MR sequences with the ability to alter MR parameters via an intuitive GUI control panel. The application allows the user to gain real-time understanding of image transformation when varying these said parameters by examining the resulting images. Lab procedures can be loaded and displayed for more directed study. The panel is also configurable, providing a simple interface for elementary labs or a full array of controls for the expert user. Conclusion: Our introduction of MrRSCAL, which is readily available to users with a current laptop or workstation, allows for individual or group study of MR image acquisition with immediate educational feedback as the MR parameters are manipulated. MrRSCAL can be used at any time and any place once installed, offering a new tool for reviewing relaxometric and artifact principles when studying for boards or investigating properties of a pulse sequence. This tool promises to be extremely useful in conveying traditionally difficult and abstract concepts involved with MR to the radiology resident and other medical professionals at large.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Vickers hardness prediction from machine learning methods

Abstract The search for new superhard materials is of great interest for extreme industrial applications. However, the theoretical prediction of hardness is still a challenge for the scientific community, given the difficulty of modeling plastic behavior of solids. Different hardness models have been proposed over the years. Still, they are either too complicated to use, inaccurate when extrapolating to a wide variety of solids or require coding knowledge. In this investigation, we built a successful machine learning model that implements Gradient Boosting Regressor (GBR) to predict hardness and uses the mechanical properties of a solid (bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio) as input variables. The model was trained with an experimental Vickers hardness database of 143 materials, assuring various kinds of compounds. The input properties were calculated from the theoretical elastic tensor. The Materials Project’s database was explored to search for new superhard materials, and our results are in good agreement with the experimental data available. Other alternative models to compute hardness from mechanical properties are also discussed in this work. Our results are available in a free-access easy to use online application to be further used in future studies of new materials at www.hardnesscalculator.com .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

NASA NTRS2002

Low-Cost Evaluation of EO-1 Hyperion and ALI for Detection and Biophysical Characterization of Forest Logging in Amazonia (NCC5-481)

Major uncertainties exist regarding the rate and intensity of logging in tropical forests worldwide: these uncertainties severely limit economic, ecological, and biogeochemical analyses of these regions. Recent sawmill surveys in the Amazon region of Brazil show that the area logged is nearly equal to total area deforested annually, but conversion of survey data to forest area, forest structural damage, and biomass estimates requires multiple assumptions about logging practices. Remote sensing could provide an independent means to monitor logging activity and to estimate the biophysical consequences of this land use. Previous studies have demonstrated that the detection of logging in Amazon forests is difficult and no studies have developed either the quantitative physical basis or remote sensing approaches needed to estimate the effects of various logging regimes on forest structure. A major reason for these limitations has been a lack of sufficient, well-calibrated optical satellite data, which in turn, has impeded the development and use of physically-based, quantitative approaches for detection and structural characterization of forest logging regimes. We propose to use data from the EO-1 Hyperion imaging spectrometer to greatly increase our ability to estimate the presence and structural attributes of selective logging in the Amazon Basin. Our approach is based on four "biogeophysical indicators" not yet derived simultaneously from any satellite sensor: 1) green canopy leaf area index; 2) degree of shadowing; 3) presence of exposed soil and; 4) non-photosynthetic vegetation material. Airborne, field and modeling studies have shown that the optical reflectance continuum (400-2500 nm) contains sufficient information to derive estimates of each of these indicators. Our ongoing studies in the eastern Amazon basin also suggest that these four indicators are sensitive to logging intensity. Satellite-based estimates of these indicators should provide a means to quantify both the presence and degree of structural disturbance caused by various logging regimes. Our quantitative assessment of Hyperion hyperspectral and ALI multi-spectral data for the detection and structural characterization of selective logging in Amazonia will benefit from data collected through an ongoing project run by the Tropical Forest Foundation, within which we have developed a study of the canopy and landscape biophysics of conventional and reduced-impact logging. We will add to our base of forest structural information in concert with an EO-1 overpass. Using a photon transport model inversion technique that accounts for non-linear mixing of the four biogeophysical indicators, we will estimate these parameters across a gradient of selective logging intensity provided by conventional and reduced impact logging sites. We will also compare our physical ly-based approach to both conventional (e.g., NDVI) and novel (e.g., SWIR-channel) vegetation indices as well as to linear mixture modeling methods. We will cross-compare these approaches using Hyperion and ALI imagers to determine the strengths and limitations of these two sensors for applications of forest biophysics. This effort will yield the first physical ly-based, quantitative analysis of the detection and intensity of selective logging in Amazonia, comparing hyperspectral and improved multi-spectral approaches as well as inverse modeling, linear mixture modeling, and vegetation index techniques.

Asner, Gregory P.↗

NASA NTRS2005

Robotic Technologies for Surveying Habitats and Seeking Evidence of Life: Results from the 2004 Field Experiments of the "Life in the Atacama" Project

The Chilean Atacama Desert is the most arid region on Earth and in several ways analogous to Mars. Evidence suggests that the interior of the Atacama is lifeless, yet where the desert meets the Pacific coastal range dessication-tolerant microorganisms are known to exist. The gradient of biodiversity and habitats in the Atacama's subregions remain unexplored and are the focus of the Life in the Atacama project. Our field investigation attempts to bring further scientific understanding of the Atacama as a habitat for life through the creation of robotic astrobiology. This involves capabilities for autonomously traversing hundreds of kilometers while deploying sensors to survey the varying geologic and biologic properties of the environment, Fig. 1. Our goal is to make genuine discoveries about the limits of life on Earth and to generate knowledge about life in extreme environments that can be applied to future planetary missions. Through these experiments we also hope to develop and practice the methods by which a rover might best be employed to survey desert terrain in search of the habitats in which life can survive, or may have in the past.

Wettergreen, D.↗

DOE OSTI2016

SU-F-I-10: Spatially Local Statistics for Adaptive Image Filtering

Purpose: To facilitate adaptive image filtering operations, addressing spatial variations in both noise and signal. Such issues are prevalent in cone-beam projections, where physical effects such as X-ray scattering result in spatially variant noise, violating common assumptions of homogeneous noise and challenging conventional filtering approaches to signal extraction and noise suppression. Methods: We present a computational mechanism for probing into and quantifying the spatial variance of noise throughout an image. The mechanism builds a pyramid of local statistics at multiple spatial scales; local statistical information at each scale includes (weighted) mean, median, standard deviation, median absolute deviation, as well as histogram or dynamic range after local mean/median shifting. Based on inter-scale differences of local statistics, the spatial scope of distinguishable noise variation is detected in a semi- or un-supervised manner. Additionally, we propose and demonstrate the incorporation of such information in globally parametrized (i.e., non-adaptive) filters, effectively transforming the latter into spatially adaptive filters. The multi-scale mechanism is materialized by efficient algorithms and implemented in parallel CPU/GPU architectures. Results: We demonstrate the impact of local statistics for adaptive image processing and analysis using cone-beam projections of a Catphan phantom, fitted within an annulus to increase X-ray scattering. The effective spatial scope of local statistics calculations is shown to vary throughout the image domain, necessitating multi-scale noise and signal structure analysis. Filtering results with and without spatial filter adaptation are compared visually, illustrating improvements in imaging signal extraction and noise suppression, and in preserving information in low-contrast regions. Conclusion: Local image statistics can be incorporated in filtering operations to equip them with spatial adaptivity to spatial signal/noise variations. An efficient multi-scale computational mechanism is developed to curtail processing latency. Spatially adaptive filtering may impact subsequent processing tasks such as reconstruction and numerical gradient computations for deformable registration. NIH Grant No. R01-184173.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

SU-E-J-37: Combining Proton Radiography and X-Ray CT Information to Better Estimate Relative Proton Stopping Power in a Clinical Environment

Purpose: In standard proton therapy clinical practice, proton stopping power uncertainties are in the order of 3.5%, which affects the ability of placing the proton Bragg peak at the edge of the tumor. The innovating idea of this project is to approach the uncertainty problem in RSP by using combined information from X-ray CT and proton radiography along a few beam angles. In addition, this project aims to quantify the systematic error introduced by the theoretical models (Janni, ICRU49, Bischel) for proton stopping power in media. Methods: A 3D phantom of 36 cm3 composed of 9 materials randomly placed is created. Measured RSP values are obtained using a Gammex phantom with a proton beam. Theoretical RSP values are calculated with Beth-Block equation in combination with three databases (Janni, ICRU49 and Bischel). Clinical RSP errors are simulated by introducing a systematic (1.5%, 2.5%, 3.5%) and a random error (+/−0.5%) to the theoretical RSP. A ray-tracing algorithm uses each of these RSP tables to calculate energy loss for proton crossing the phantom through various directions. For each direction, gradient descent (GD) method is done on the clinical RSP table to minimize the residual energy difference between the simulation with clinical RSP and with theoretical RSP. The possibility of a systematic material dependent error is investigated by comparing measured RSP to theoretical RSP as calculated from the three models. Results: Using 10,000 iterations on GD algorithm, RSP differences between theoretical values and clinical RSP have converged (<1%) for each error introduced. Results produced with ICRU49 have the smallest average difference (0.021%) to the measured RSP. Janni (1.168%) and Bischel (−0.372%) database shows larger systematic errors. Conclusion: Based on these results, ray-tracing optimisation using information from proton radiography and X-ray CT demonstrates a potential to improve the proton range accuracy in a clinical environment.

07 ISOTOPE AND RADIATION SOURCES↗

NASA NTRS1983

Dual laser optical system and method for studying fluid flow

A dual laser optical system and method is disclosed for visualization of phenomena in transport substances which induce refractive index gradients such as fluid flow and pressure and temperature gradients in fluids and gases. Two images representing mutually perpendicular components of refractive index gradients may be viewed simultaneously on screen. Two lasers having wave lengths in the visible range but separated by about 1000 angstroms are utilized to provide beams which are collimated into a beam containing components of the different wave lengths. The collimated beam is passed through a test volume of the transparent substance. The collimated beam is then separated into components of the different wave lengths and focused onto a pair of knife edges arranged mutually perpendicular to produce and project images onto the screen.

Owen, R. B.↗