Search NASA⌕ Search

SEARCH · Search NASA

Results for “image based modeling”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Extra-nuclear starbursts: young luminous Hinge clumps in interacting galaxies

Hinge clumps are luminous knots of star formation near the base of tidal features in some interacting galaxies. We use archival Hubble Space Telescope (HST) UV/optical/IR images and Chandra X-ray maps along with Galaxy Evolution Explorer UV, Spitzer IR, and ground-based optical/near-IR images to investigate the star forming properties in a sample of 12 hinge clumps in five interacting galaxies. The most extreme of these hinge clumps have star formation rates of 1-9 M {sub ☉} yr{sup –1}, comparable to or larger than the 'overlap' region of intense star formation between the two disks of the colliding galaxy system the Antennae. In the HST images, we have found remarkably large and luminous sources at the centers of these hinge clumps. These objects are much larger and more luminous than typical 'super star clusters' in interacting galaxies, and are sometimes embedded in a linear ridge of fainter star clusters, consistent with star formation along a narrow caustic. These central sources have FWHM diameters of ∼70 pc, compared to ∼3 pc in 'ordinary' super star clusters. Their absolute I magnitudes range from M{sub I} ∼ – 12.2 to –16.5; thus, if they are individual star clusters they would lie near the top of the 'super star cluster' luminosity function of star clusters. These sources may not be individual star clusters, but instead may be tightly packed groups of clusters that are blended together in the HST images. Comparison to population synthesis modeling indicates that the hinge clumps contain a range of stellar ages. This is consistent with expectations based on models of galaxy interactions, which suggest that star formation may be prolonged in these regions. In the Chandra images, we have found strong X-ray emission from several of these hinge clumps. In most cases, this emission is well-resolved with Chandra and has a thermal X-ray spectrum, thus it is likely due to hot gas associated with the star formation. The ratio of the extinction-corrected diffuse X-ray luminosity to the mechanical energy rate (the X-ray production efficiency) for the hinge clumps is similar to that in the Antennae galaxies, but higher than those for regions in the normal spiral galaxy NGC 2403. Two of the hinge clumps have point-like X-ray emission much brighter than expected for hot gas; these sources are likely 'ultra-luminous X-ray sources' due to accretion disks around black holes. The most extreme of these sources, in Arp 240, has a hard X-ray spectrum and an absorbed X-ray luminosity of ∼2 × 10{sup 41} erg s{sup –1}; this is above the luminosity expected by single high mass X-ray binaries (HMXBs), thus it may be either a collection of HMXBs or an intermediate mass black hole (≥80 M {sub ☉}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

WE-FG-202-09: Voxel-Level Analysis of Adverse Treatment Response in Pediatric Patients Treated for Ependymoma with Passive Scattering Proton Therapy

Purpose: We identified patients treated for ependymoma with passive scattering proton therapy who subsequently developed treatment-related imaging changes on MRI. We sought to determine if there is any spatial correlation between imaged response, dose, and LET. Methods: A group of 14 patients treated for ependymoma were identified as having post-treatment MR imaging changes observable as T2-FLAIR hyperintensity with or without enhancement on T1 post-contrast sequences. MR images were registered with treatment planning CT images and regions of treatment-related change contoured by a practicing radiation oncologist. The contoured regions were identified as response with voxels represented as 1 while voxels within the brain outside of the response region were represented as 0. An in-house Monte Carlo system was used to recalculate treatment plans to obtain dose and LET information. Voxels were binned according to LET values in 0.3 keV µm{sup −1} bins. Dose and corresponding response value (0 or 1) for each voxel for a given LET bin were then plotted and fit with the Lyman-Kutcher-Burman dose response model to determine TD{sub 50} and m parameters for each LET value. Response parameters from all patients were then collated, and linear fits of the data were performed. Results: The response parameters TD50 and m both show trends with LET. Outliers were observed due to low numbers of response voxels in some cases. TD{sub 50} values decreased with LET while m increased with LET. The former result would indicate that for higher LET values, the dose is more effective, which is consistent with relative biological effectiveness (RBE) models for proton therapy. Conclusion: A novel method of voxel-level analysis of image biomarker-based adverse patient treatment response in proton therapy according to dose and LET has been presented. Fitted TD{sub 50} values show a decreasing trend with LET supporting the typical models of proton RBE. Funding provided by NIH Program Project Grant 2U19CA021239-35.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

An automated, high-throughput method for standardizing image color profiles to improve image-based plant phenotyping

High-throughput phenotyping has emerged as a powerful method for studying plant biology. Large image-based datasets are generated and analyzed with automated image analysis pipelines. A major challenge associated with these analyses is variation in image quality that can inadvertently bias results. Images are made up of tuples of data called pixels, which consist of R, G, and B values, arranged in a grid. Many factors, for example image brightness, can influence the quality of the image that is captured. These factors alter the values of the pixels within images and consequently can bias the data and downstream analyses. Here, we provide an automated method to adjust an image-based dataset so that brightness, contrast, and color profile is standardized. The correction method is a collection of linear models that adjusts pixel tuples based on a reference panel of colors. We apply this technique to a set of images taken in a high-throughput imaging facility and successfully detect variance within the image dataset. In this case, variation resulted from temperature-dependent light intensity throughout the experiment. Using this correction method, we were able to standardize images throughout the dataset, and we show that this correction enhanced our ability to accurately quantify morphological measurements within each image. We implement this technique in a high-throughput pipeline available with this paper, and it is also implemented in PlantCV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2025

Boosting Noise2Inverse via enhanced model selection for denoising computed tomography data

Synchrotron-based x-ray tomographic imaging enables the examination of the internal structure of materials at high spatial and temporal resolution. Experimental constraints can impose dose and time limits on the measurements, introducing a higher level of noise and artifacts in the reconstructed images. Deep learning has emerged as a powerful tool to remove noise from reconstructed images. Recently, the Noise2Inverse method was designed specifically for denoising reconstructed images without requiring paired noisy and clean images. This method creates multiple statistically independent reconstructions used to pair the data in which training involves transforming one reconstruction into the other, and vice versa. Originally designed to be used after a fixed number of epochs, we see in practice that this approach may not produce the optimal model and may unnecessarily waste computational resources. Therefore, we propose an alternative method of identifying the best model during training that aligns with the Noise2Inverse method. During validation, we compare the model output of the multiple reconstructions among each other. We hypothesize that the best model is the one that produces images with the highest similarity, implying a convergence in the predicted material properties and absorption values. To compare model outputs, we consider the absolute error, square error, structural similarity index (SSIM), peak signal-to-noise ratio (PSNR), and cosine similarity. We evaluate our method on two simulated tomography datasets and two, real-world, low-contrast, high-energy x-ray tomography datasets. We show our approach is more effective at determining the best model, up to an increase of 12.50% and 12.53% in SSIM and PSNR, respectively, while only requiring a fifth of the training time compared to the original approach.

CT↗
DOE OSTI2022

Minimizing thickness variation in monolithic U-10Mo fuel foil and Zr interlayer during hot rolling: A microstructure-based finite element method analysis

Low-enriched uranium alloyed with 10 wt. % molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium fuel for the United States’ high performance research reactors. The monolithic U-10Mo fuel plate consists of a metallic U-10Mo fuel foil with a 25 µm Zr interlayer and a relatively thick cladding of aluminum alloy 6061. The Zr interlayer is typically applied during the hot co-rolling process, and this process dictates the uniformity of the Zr interlayer. Thickness variation observed in the U-10Mo and Zr interlayer has been attributed to several sources: the initial grain size of the U-10Mo castings, can materials, rolling temperature, inhomogeneous molybdenum content, and porosity in the cast U-10Mo. This thickness variation limits the ability to meet the dimensional specification; thus, a better understanding of the factors causing the nonuniform thickness is needed. In this work, we used a novel, microstructure-based finite element method to model the hot rolling process to address these concerns. Grain microstructures in U-10Mo were tessellated and explicitly considered in the finite element model. Each grain was assigned a random material property to mimic the grain strength variations induced by different grain orientations. Simulations were performed using six steel can thicknesses, four grain sizes, and with or without a Zr interlayer to investigate the influences of those variables on the thickness nonuniformity. The simulation results showed that a thinner steel can and finer U-10Mo grain size reduce thickness variations in both the U-10Mo fuel foil and Zr interlayer. The direct findings from the simulations and analysis can be used to optimize the hot rolling schedule, reduce fabrication defects, and meet the dimensional specifications. The proposed microstructure-based finite element model can be also coupled with experimental microstructure characterization data, images, and models to simulate multi-pass hot rolling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

New Orthophoto Generation Strategies from UAV and Ground Remote Sensing Platforms for High-Throughput Phenotyping

Remote sensing platforms have become an effective data acquisition tool for digital agriculture. Imaging sensors onboard unmanned aerial vehicles (UAVs) and tractors are providing unprecedented high-geometric-resolution data for several crop phenotyping activities (e.g., canopy cover estimation, plant localization, and flowering date identification). Among potential products, orthophotos play an important role in agricultural management. Traditional orthophoto generation strategies suffer from several artifacts (e.g., double mapping, excessive pixilation, and seamline distortions). The above problems are more pronounced when dealing with mid- to late-season imagery, which is often used for establishing flowering date (e.g., tassel and panicle detection for maize and sorghum crops, respectively). In response to these challenges, this paper introduces new strategies for generating orthophotos that are conducive to the straightforward detection of tassels and panicles. The orthophoto generation strategies are valid for both frame and push-broom imaging systems. The target function of these strategies is striking a balance between the improved visual appearance of tassels/panicles and their geolocation accuracy. The new strategies are based on generating a smooth digital surface model (DSM) that maintains the geolocation quality along the plant rows while reducing double mapping and pixilation artifacts. Moreover, seamline control strategies are applied to avoid having seamline distortions at locations where the tassels and panicles are expected. The quality of generated orthophotos is evaluated through visual inspection as well as quantitative assessment of the degree of similarity between the generated orthophotos and original images. Several experimental results from both UAV and ground platforms show that the proposed strategies do improve the visual quality of derived orthophotos while maintaining the geolocation accuracy at tassel/panicle locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Vector elastic deconvolution migration with dual wavefield decomposition

Here, we report elastic-wave imaging using multicomponent data can provide more useful subsurface information than acoustic-wave imaging, but it is usually algorithmically challenging. We have developed a vector elastic deconvolution migration method for high-resolution imaging of subsurface structures in isotropic and anisotropic elastic media. Our new method uses a vector deconvolution imaging condition based on dual wavefield decomposition, including an explicit directional wavefield separation using the Hilbert transform and a P/S vector wavefield decomposition using the low-rank decomposition method. Using three elastic models, we numerically determine that our new method produces notably higher resolution and more amplitude-balanced elastic images compared with a crosscorrelation-based vector elastic reverse time migration method.

58 GEOSCIENCES↗
DOE OSTI2022

Pulsed Thermal Tomography Nondestructive Examination of Additively Manufactured Reactor Materials and Components (Final Technical Report)

Metal Additive Manufacturing (AM) is a promising method for cost-efficient fabrication of complex shape structures for applications in harsh environment, such as in a nuclear reactor. However, internal defects (pores) occur in high-strength AM alloys, which are manufactured with Laser Powder Bed Fusion (LPBF) AM method. Pulsed Infrared Thermography (PIT) is an efficient nondestructive evaluation (NDE) method to examine actual structures, because this method offers one-sided non-contact measurements, and fast processing of large sample areas. However, imaging of material defects, particularly defects with sizes at microscopic level, is challenging. In this report, we benchmark the performance of several Unsupervised Learning (UL) algorithms designed to enhance imaging of microscopic defects in metals with PIT. UL aims to learn the latent principal patterns (dictionaries) in PIT data to detect defects with minimal human supervision. Performance of Independent Component Analysis (ICA), Sparse Coding (SC), Principal Component Analysis (PCA) and Exploratory Factor Analysis (EFA) was compared using F-score, UL model training time and defects reconstruction time. We obtained the average F-score of 0.75, and a highest F-score of 0.89 for the EFA algorithm. Overall, EFA outperforms other UL algorithms considered in this study. In another approach, we investigate Thermal Tomography (TT), which is a computational method for reconstruction of depth profile of internal material defects from PIT nondestructive evaluation (NDE). TT algorithm obtains depth reconstructions of thermal effusivity, which has been shown to provide visualization of subsurface internals defects in metals. In many applications, one needs to determine the defect shape and orientation from reconstructed effusivity images. Interpretation of TT images is non-trivial because of blurring, which increases with depth due to heat diffusion-based nature of image formation. We have developed a deep learning convolutional neural network (CNN) to classify size and orientation of subsurface material defects in TT images. CNN was trained with TT images produced with computer simulations of 2D metallic structures (thin plates) containing elliptical subsurface voids. Performance of CNN was investigated using test TT images developed with computer simulations of plates containing elliptical defects, and defects with shape imported from scanning electron microscopy (SEM) images. CNN demonstrated the ability to classify radii and angular orientation of elliptical defects in previously unseen test TT images. We have also demonstrated that CNN trained on TT images of elliptical defects is capable of classifying shape and orientation of irregular defects. Training the CNN on irregular defect shapes instead of on elliptical shapes would make the resulting classifications more descriptive of actual defect shapes. However, this requires a much higher volume of SEM images of material defects, which are difficult to obtain because of random occurrence of defects in LPBF. To address this challenge, we developed a generative adversarial network (GAN) to augment the existing dataset of SEM defect images. The GAN model is demonstrated to create novel yet realistic defect shapes that can be used as input for simulated PTT images to train CNN. We also investigate several approaches based on Gaussian Random Circle and Bezier Curves for constructing parametric models of irregular-shape defects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Physics-informed semantic inpainting: Application to geostatistical modeling

A fundamental problem in geostatistical modeling is to infer the heterogeneous geological field based on limited measurements and some prior spatial statistics. Semantic inpainting, a technique for image processing using deep generative models, has been recently applied for this purpose, demonstrating its effectiveness in dealing with complex spatial patterns. However, the original semantic inpainting framework incorporates only information from direct measurements, while in geostatistics indirect measurements are often plentiful. In this work, to overcome this limitation, we propose a physics-informed semantic inpainting framework, employing the Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) and jointly incorporating the direct and indirect measurements by exploiting the underlying physical laws. Our simulation results for a high-dimensional problem with 512 dimensions show that in the new method, the physical conservation laws are satisfied and contribute in enhancing the inpainting performance compared to using only the direct measurements.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2008

Cammp Team

The International Space Station (ISS) Configuration Analysis Modeling and Mass Properties (CAMMP) Team is presenting a demo of certain CAMMP capabilities at a Booz Allen Hamilton conference in San Antonio. The team will be showing pictures of low fidelity, simplified ISS models, but no dimensions or technical data. The presentation will include a brief description of the contract and task, description and picture of the Topology, description of Generic Ground Rules and Constraints (GGR&C), description of Stage Analysis with constraints applied, and wrap up with description of other tasks such as Special Studies, Cable Routing, etc. The models include conceptual Crew Exploration Vehicle (CEV) and Lunar Lander images and animations created for promotional purposes, which are based entirely on public domain conceptual images from public NASA web sites and publicly available magazine articles and are not based on any actual designs, measurements, or 3D models. Conceptual Mars rover and lander are completely conceptual and are not based on any NASA designs or data. The demonstration includes High Fidelity Computer Aided Design (CAD) models of ISS provided by the ISS 3D CAD Team which will be used in a visual display to demonstrate the capabilities of the Teamcenter Visualization software. The demonstration will include 3D views of the CAD models including random measurements that will be taken to demonstrate the measurement tool. A 3D PDF file will be demonstrated of the Blue Book fidelity assembly complete model with no vehicles attached. The 3D zoom and rotation will be displayed as well as random measurements from the measurement tool. The External Configuration Analysis and Tracking Tool (ExCATT) Microsoft Access Database will be demonstrated to show its capabilities to organize and track hardware on ISS. The data included will be part numbers, serial numbers, historical, current, and future locations, of external hardware components on station. It includes dates of all external ISS events and flights and the associated hardware changes for each event. The hardware location information does not always reveal the exact location of the hardware, only the general location. In some cases the location is a module or carrier, in other cases it is a WIF socket, handrail, or attach point. Only small portions of the data will be displayed for demonstration purposes.

Evertt, Shonn F.↗
DOE OSTI2012

Predictive Modeling of Human Perception Subjectivity: Feasibility Study of Mammographic Lesion Similarity

The majority of clinical content-based image retrieval (CBIR) studies disregard human perception subjectivity, aiming to duplicate the consensus expert assessment of the visual similarity on example cases. The purpose of our study is twofold: (i) discern better the extent of human perception subjectivity when assessing the visual similarity of two images with similar semantic content, and (ii) explore the feasibility of personalized predictive modeling of visual similarity. We conducted a human observer study in which five observers of various expertise were shown ninety-nine triplets of mammographic masses with similar BI-RADS descriptors and were asked to select the two masses with the highest visual relevance. Pairwise agreement ranged between poor and fair among the five observers, as assessed by the kappa statistic. The observers' self-consistency rate was remarkably low, based on repeated questions where either the orientation or the presentation order of a mass was changed. Various machine learning algorithms were explored to determine whether they can predict each observer's personalized selection using textural features. Many algorithms performed with accuracy that exceeded each observer's self-consistency rate, as determined using a cross-validation scheme. This accuracy was statistically significantly higher than would be expected by chance alone (two-tailed p-value ranged between 0.001 and 0.01 for all five personalized models). The study confirmed that human perception subjectivity should be taken into account when developing CBIR-based medical applications.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2012

Space-based Observational Constraints for 1-D Plume Rise Models

We use a space-based plume height climatology derived from observations made by the Multi-angle Imaging SpectroRadiometer (MISR) instrument aboard the NASA Terra satellite to evaluate the ability of a plume-rise model currently embedded in several atmospheric chemical transport models (CTMs) to produce accurate smoke injection heights. We initialize the plume-rise model with assimilated meteorological fields from the NASA Goddard Earth Observing System and estimated fuel moisture content at the location and time of the MISR measurements. Fire properties that drive the plume-rise model are difficult to estimate and we test the model with four estimates for active fire area and four for total heat flux, obtained using empirical data and Moderate Resolution Imaging Spectroradiometer (MODIS) re radiative power (FRP) thermal anomalies available for each MISR plume. We show that the model is not able to reproduce the plume heights observed by MISR over the range of conditions studied (maximum r2 obtained in all configurations is 0.3). The model also fails to determine which plumes are in the free troposphere (according to MISR), key information needed for atmospheric models to simulate properly smoke dispersion. We conclude that embedding a plume-rise model using currently available re constraints in large-scale atmospheric studies remains a difficult proposition. However, we demonstrate the degree to which the fire dynamical heat flux (related to active fire area and sensible heat flux), and atmospheric stability structure influence plume rise, although other factors less well constrained (e.g., entrainment) may also be significant. Using atmospheric stability conditions, MODIS FRP, and MISR plume heights, we offer some constraints on the main physical factors that drive smoke plume rise. We find that smoke plumes reaching high altitudes are characterized by higher FRP and weaker atmospheric stability conditions than those at low altitude, which tend to remain confined below the BL, consistent with earlier results. We propose two simplified parameterizations for computing injection heights for fires in CTMs and discuss current challenges to representing plume injection heights in large scale atmospheric models.

Martin, Maria Val↗
DOE OSTI2013

2D/3D registration algorithm for lung brachytherapy

Purpose: A 2D/3D registration algorithm is proposed for registering orthogonal x-ray images with a diagnostic CT volume for high dose rate (HDR) lung brachytherapy. Methods: The algorithm utilizes a rigid registration model based on a pixel/voxel intensity matching approach. To achieve accurate registration, a robust similarity measure combining normalized mutual information, image gradient, and intensity difference was developed. The algorithm was validated using a simple body and anthropomorphic phantoms. Transfer catheters were placed inside the phantoms to simulate the unique image features observed during treatment. The algorithm sensitivity to various degrees of initial misregistration and to the presence of foreign objects, such as ECG leads, was evaluated. Results: The mean registration error was 2.2 and 1.9 mm for the simple body and anthropomorphic phantoms, respectively. The error was comparable to the interoperator catheter digitization error of 1.6 mm. Preliminary analysis of data acquired from four patients indicated a mean registration error of 4.2 mm. Conclusions: Results obtained using the proposed algorithm are clinically acceptable especially considering the complications normally encountered when imaging during lung HDR brachytherapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Complete Performance Comparison Between the Optimized Image Construction Algorithm (U-MBIR) and the Existing Reconstruction Algorithm for Detecting Defects and Damage in Concrete

Reinforced concrete (RC) is a composite material subjected to mechanical, thermal, and chemical loads throughout its service life. Because of these external stressors and the susceptibility of RC structural members to shrinkage and microcracking, the material degrades throughout its life cycle. This deterioration can lead to a decrease in member capacity and, ultimately, poses a threat to the structural integrity. Thus, it is crucial that the damage caused by aging and degradation be monitored and assessed at regular intervals throughout the material’s service life. Since coring of the material is typically not feasible for in-service structural systems, non-destructive evaluation (NDE) methodologies are used to assess remaining structural capacity. NDE methods enable surface and subsurface examination without damaging or degrading the medium. Moreover, RC is a critical component of nuclear power plants; thus, its safety and reliability must be thoroughly examined throughout the life cycle of the structural system. Ultrasonic measurements have been an industry standard for both surface and subsurface inspections. To this end, Oak Ridge National Laboratory (ORNL) has researched and developed advanced image reconstruction algorithms to capture internal damage. The results and discussion presented herein summarize the current state of the ultrasonic model–based iterative reconstruction (U-MBIR) algorithm developed at ORNL. More specifically, this report presents a comparison between reconstruction images produced via a widely employed ultrasonic NDE technique—the synthetic aperture focusing technique (SAFT)—and the ORNL-developed U-MBIR algorithm. These NDE methodologies are demonstrated using ultrasonic data collected from four concrete specimens. Overall, the U-MBIR algorithm eliminates artifacts and noise that are typically present within the SAFT reconstructions, and it shows defects and anomalies more clearly than the SAFT images. In conclusion, this algorithm is suitable for identifying concrete defects, although more improvements and optimization could be done to better define internal defects.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Distributed Assimilation of Satellite-based Snow Extent for Improving Simulated Streamflow in Mountainous, Dense Forests: An Example Over the DMIP2 Western Basins

Snow cover area affects snowmelt, soil moisture, evapotranspiration, and ultimately streamflow. For the Distributed Model Intercomparison Project - Phase 2 Western basins, we assimilate satellite-based fractional snow cover area (fSCA) from the Moderate Resolution Imaging Spectroradiometer, or MODIS, into the National Weather Service (NWS) SNOW-17 model. This model is coupled with the NWS Sacramento Heat Transfer (SAC-HT) model inside the National Aeronautics and Space Administration's (NASA) Land Information System. SNOW-17 computes fSCA from snow water equivalent (SWE) values using an areal depletion curve. Using a direct insertion, we assimilate fSCAs in two fully distributed ways: 1) we update the curve by attempting SWE preservation, and 2) we reconstruct SWEs using the curve. The preceding are refinements of an existing simple, conceptually-guided NWS algorithm. Satellite fSCA over dense forests inadequately accounts for below-canopy snow, degrading simulated streamflow upon assimilation during snowmelt. Accordingly, we implement a below-canopy allowance during assimilation. This simplistic allowance and direct insertion are found to be inadequate for improving calibrated results, still degrading them as mentioned above. However, for streamflow volume for the uncalibrated runs, we obtain: (1) substantial to major improvements (64-81 %) as a percentage of the control run residuals (or distance from observations), and (2) minor improvements (16-22 %) as a percentage of observed values. We highlight the need for detailed representations of canopy-snow optical radiative transfer processes in mountainous, dense forest regions if assimilation-based improvements are to be seen in calibrated runs over these areas.

Yatheendradas, Soni↗
DOE OSTI2023

Application of unsupervised deep learning to image segmentation and in-situ contact angle measurements in a CO 2 -water-rock system

Rock surface wettability is a critical property that regulates multiphase flows in porous media, which can be quantified using the surface contact angle (CA). X-ray micro-computed tomography (μCT) provides an effective approach to in-situ measurements of surface CAs. However, the CA measurement accuracy depends significantly on the quality of CT image segmentation, which is the clustering of CT pixels into separate phases. Inspired by this, we developed a deep learning (DL)-based CA measurement workflow. Motivated by the recent tremendous progress in unsupervised learning techniques and aiming to avoid expensive manual data annotations, an unsupervised DL pipeline for CT image segmentation was proposed and implemented, which includes unsupervised model training and post-processing. The unsupervised model training was driven by a novel loss function constrained with feature similarity and spatial continuity and implemented by iterative forward and backward paths; the former clustered the pixel-wise feature vectors extracted by convolution neural networks, whereas the latter updated the parameters using gradient descent. An over-segmentation strategy was adopted for model training. The post-processing steps based on agglomerative hierarchical clustering (AHC) were implemented to further merge the over-segmented model output to the desired cluster number, which is intended to improve the efficiency of image segmentation. The developed unsupervised DL pipeline was compared with other commonly-used image segmentation methods using pixel-wise and physics-based evaluation metrics on a synthetic raw-image dataset, which had a known ground truth. The unsupervised DL pipeline showed the best performance. Next, the segmented images were input to an automatic CA measurement tool, and the results were validated by comparisons with manual measurements. The CA values from the manual and automatic measurements showed similar distributions and statistical properties. The automatic measurement demonstrated a wider spectrum because of the much larger number of measurement data points. The primary novelty of the unsupervised DL pipeline developed in this study lies in the novel loss function and the over-segmentation strategy associated with AHC post-processing. Finally, the workflow has been proven an efficient tool for pore-scale wettability characterization, which has a wide range of applications in fundamental studies of multiphase flows in natural porous media, which have critical implications to geological carbon sequestration, hydrocarbon energy recovery, and contaminant transport in groundwater.

42 ENGINEERING↗
NASA NTRS2021

JWST/MIRI Simulated Imaging: Insights into Obscured Star Formation and AGNs for Distant Galaxies in Deep Surveys

The James Webb Space Telescope MIRI instrument will revolutionize extragalactic astronomy with unprecedentedsensitivity and angular resolution in mid-IR. Here we assess the potential of MIRI photometry to constrain galaxyproperties in the Cosmic Evolution Early Release Science(CEERS)survey. We derive estimated MIRIfluxes fromthe spectral energy distributions(SEDs)of real sources that fall in a planned MIRI pointing. We also obtain MIRIfluxes for hypothetical active galactic nucleus(AGN)–galaxy mixed models varying the AGN fractionalcontribution to the total IR luminosity(fracAGN). Based on these modelfluxes, we simulate CEERS imaging(3.6 hrexposure)in six bands from F770W to F2100W usingMIRISIMand reduce these data usingJWST PIPELINE.Weperform point-spread-function-matched photometry withTPHOTandfit the source SEDs withX-CIGALE,simultaneously modeling photometric redshift and other physical properties. Adding the MIRI data, the accuracyof both redshift and fracAGNis generally improved by factors of2 for all sources atz3. Notably, for pure-galaxy inputs(fracAGN=0), the accuracy of fracAGNis improved by∼100 times thanks to MIRI. The simulatedCEERS MIRI data are slightly more sensitive to AGN detections than the deepest X-ray survey, based on theempiricalLX–L6μmrelation. Like X-ray observations, MIRI can also be used to constrain the AGN accretionpower(accuracy≈0.3 dex). Our work demonstrates that MIRI will be able to place strong constraints on the mid-IR luminosities from star formation and AGNs and thereby facilitate studies of the galaxy/AGN coevolution.

G Yang↗