SEARCH · Search NASA

Results for “generator parameter estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Turbulence-induced bias in time-averaged laser absorption tomography of correlated concentration and temperature fields with a first-order correction

The influence of correlated scalar fluctuations on time-averaged laser absorption tomography measurements of temperature and species in a piloted turbulent premixed flame was examined using a coupled spectroscopic and fluid-dynamic analysis. To understand bias associated with turbulence, spatio-temporally resolved temperature and species mole fraction profiles predicted by large eddy simulations (LES) were used to synthetically generate time-resolved line-of-sight absorption measurements at short time scales (microsecond) to reflect the unsteady nature of a canonical jet burner across various transverse measurement planes. Inversion methods were employed on the time-averaged line-of-sight data to produce radially-resolved temperature and mole fraction profiles, analogous to those produced by laser absorption tomography performed on a time-averaged axisymmetric flowfield. It is shown that bias in the measurements compared to true time-averaged scalar fields is a function primarily of temperature dependence in absorptivity and non-zero correlation between temperature and species concentration scalars. Finally, a first-order correction to tomography measurements is proposed to account for the bias based on estimated correlations and the known spectroscopic parameters of the probed absorption transitions.

42 ENGINEERING↗

DOE OSTI2021

A nudged hybrid analysis and modeling approach for realtime wake-vortex transport and decay prediction

Here, we put forth a long short-term memory (LSTM) nudging framework for the enhancement of reduced order models (ROMs) of fluid flows utilizing noisy measurements for air traffic improvements. Toward emerging applications of digital twins in aviation, the proposed approach allows for constructing a realtime predictive tool for wake-vortex transport and decay systems. We build on the fact that in realistic application, there are uncertainties in initial and boundary conditions, model parameters, as well as measurements. Moreover, conventional nonlinear ROMs based on Galerkin projection (GROMs) suffer from imperfection and solution instabilities, especially for advection-dominated flows with slow decay in the Kolmogorov n-width. In the presented LSTM nudging (LSTM-N) approach, we fuse forecasts from a combination of imperfect GROM and uncertain state estimates, with sparse Eulerian sensor measurements to provide more reliable predictions in a dynamical data assimilation framework. We illustrate our concept by solving the two-dimensional vorticity transport equation. We investigate the effects of measurements noise and state estimate uncertainty on the performance of the LSTM-N behavior. We also demonstrate that it can sufficiently handle different levels of temporal and spatial measurement sparsity, and offer a huge potential in developing next-generation digital twin technologies for aerospace applications.

42 ENGINEERING↗

NASA NTRS2011

Method and System for Temporal Filtering in Video Compression Systems

Three related innovations combine improved non-linear motion estimation, video coding, and video compression. The first system comprises a method in which side information is generated using an adaptive, non-linear motion model. This method enables extrapolating and interpolating a visual signal, including determining the first motion vector between the first pixel position in a first image to a second pixel position in a second image; determining a second motion vector between the second pixel position in the second image and a third pixel position in a third image; determining a third motion vector between the first pixel position in the first image and the second pixel position in the second image, the second pixel position in the second image, and the third pixel position in the third image using a non-linear model; and determining a position of the fourth pixel in a fourth image based upon the third motion vector. For the video compression element, the video encoder has low computational complexity and high compression efficiency. The disclosed system comprises a video encoder and a decoder. The encoder converts the source frame into a space-frequency representation, estimates the conditional statistics of at least one vector of space-frequency coefficients with similar frequencies, and is conditioned on previously encoded data. It estimates an encoding rate based on the conditional statistics and applies a Slepian-Wolf code with the computed encoding rate. The method for decoding includes generating a side-information vector of frequency coefficients based on previously decoded source data and encoder statistics and previous reconstructions of the source frequency vector. It also performs Slepian-Wolf decoding of a source frequency vector based on the generated side-information and the Slepian-Wolf code bits. The video coding element includes receiving a first reference frame having a first pixel value at a first pixel position, a second reference frame having a second pixel value at a second pixel position, and a third reference frame having a third pixel value at a third pixel position. It determines a first motion vector between the first pixel position and the second pixel position, a second motion vector between the second pixel position and the third pixel position, and a fourth pixel value for a fourth frame based upon a linear or nonlinear combination of the first pixel value, the second pixel value, and the third pixel value. A stationary filtering process determines the estimated pixel values. The parameters of the filter may be predetermined constants.

Lu, Ligang↗

NASA NTRS2002

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗

DOE OSTI2016

Investigation of Heterogeneity Effects on Calculated Sensitivities of the ZPPR-15B Assemblies

In past years, sensitivity studies of ZPR and ZPPR assemblies were carried out with deterministic codes for several purposes, such as uncertainty estimation of calculated parameter values, cross section adjustment, representativity analysis. etc. However, these sensitivity studies used simplified R-Z models of ZPR/ZPPR assemblies where the multi-group cross section data was obtained from the homogenized material compositions of each reactor zone. The main reason for this simplification was the lack of a deterministic path able to compute sensitivity coefficients for complicated 3D geometries in transport theory. Also, each ZPR and ZPPR assembly consists of several different drawers each containing different plate types and plate arrangements and it was felt that performing sensitivity calculations with an equivalent R-Z reactor model would capture most of the feedback while the heterogeneity effects were believed to be of second order. Recently a new code, BuildZPRModel, was developed at ANL that automatically generates the inputs for cross section processing and core calculation of heterogeneous ZPR and ZPPR assemblies for 3D transport calculations. Additionally, the PERSENT code (recently developed at ANL as well) is able to carry out sensitivity calculations in transport theory on 3D geometries. With the benefit of the newly developed codes, a sensitivity study was performed for the heterogeneous model of ZPPR-15B assemblies 96 and 102 generated by BuildZPRModel. The results were compared with the sensitivities computed with the use of fully homogenized R-Z reactor models, in order to quantify the impact of heterogeneity effects neglected in the analysis performed in past years.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Machine-Learning Assisted Identification of Battery Life Models

Predictive battery life models are commonly utilized to extrapolate degradation trends observed during accelerated aging tests for simulation of degradation in real-world applications. Thus, fitting accelerated aging data as accurately as possible and with low uncertainty is crucial for making believable projections of battery lifetime, but it is challenging to identify algebraic expressions that accurately fit multivariate degradation trends. A review of models published in literature reveal some common expressions for fitting calendar aging data, which is only dependent on temperature and state-of-charge, but no consistency across many models for fitting cycle aging data, indicating the need for a statistically rigorous data driven approach for developing empirical models. This talk will describe a machine-learning assisted method for identification of predictive battery life models utilizing bilevel optimization and symbolic regression. Bilevel optimization with cross-validation is used to statistically determine cell- and stress-dependent model parameters, while symbolic regression identifies both linear and multiplicative candidate expressions to predict stress-dependent degradation rates by selecting low-order subsets of features from a generated feature library. Because model expressions are identified empirically, it is crucial to ensure resulting models behave according to physical expectations, so the stability of models for interpolation or extrapolation is interrogated qualitatively through simulation and quantitatively through cross-validation and uncertainty quantification via bootstrap resampling. This model identification approach substantially improves upon models identified purely using expert judgement in terms of both accuracy and uncertainty. Model simulation and validation is then conducted by deriving a state-equation form of the predictive model, enabling simulation of battery aging under dynamic stresses. This enables validation of the predictive battery model on lab-based tests with varying conditions or on drive-cycle or application-cycle testing protocols. Parameter uncertainty can be carried forward into model simulation, giving lifetime estimates and confidence windows for cell- or system-level lifetime. The financial impact of battery model uncertainty can be estimated by incorporating uncertainty into a technoeconomic model.

battery↗

DOE OSTI2016

Investigation of Heterogeneity Effects on Calculated Sensitivities of the ZPPR-15B Assemblies

In past years, sensitivity studies of ZPR and ZPPR assemblies were carried out with deterministic codes for several purposes, such as uncertainty estimation of calculated parameter values, cross section adjustment, representativity analysis. etc.. However, these sensitivity studies used simplified R-Z models of ZPR/ZPPR assemblies where the multi-group cross section data was obtained from the homogenized material compositions of each reactor zone. The main reason for this simplification was the lack of a deterministic path able to compute sensitivity coefficients for complicated 3D geometries in transport theory. Also, each ZPR and ZPPR assembly consists of several different drawers each containing different plate types and plate arrangements and it was felt that performing sensitivity calculations with an equivalent R-Z reactor model would capture most of the feedback while the heterogeneity effects were believed to be of second order. Recently a new code, BuildZPRModel, was developed at ANL that automatically generates the inputs for cross section processing and core calculation of heterogeneous ZPR and ZPPR assemblies for 3D transport calculations. Additionally, the PERSENT code (recently developed at ANL as well) is able to carry out sensitivity calculations in transport theory on 3D geometries. With the benefit of the newly developed codes, a sensitivity study was performed for the heterogeneous model of ZPPR-15B assemblies 96 and 102 generated by BuildZPRModel. The results were compared with the sensitivities computed with the use of fully homogenized R-Z reactor models, in order to quantify the impact of heterogeneity effects neglected in the analysis performed in past years. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1996

Flow structure generated by perpendicular blade-vortex interaction and implications for helicopter noise prediction. Volume 1: Measurements

The perpendicular interaction of a streamwise vortex with an infinite span helicopter blade was modeled experimentally in incompressible flow. Three-component velocity and turbulence measurements were made using a sub-miniature four sensor hot-wire probe. Vortex core parameters (radius, peak tangential velocity, circulation, and centerline axial velocity deficit) were determined as functions of blade-vortex separation, streamwise position, blade angle of attack, vortex strength, and vortex size. The downstream development of the flow shows that the interaction of the vortex with the blade wake is the primary cause of the changes in the core parameters. The blade sheds negative vorticity into its wake as a result of the induced angle of attack generated by the passing vortex. Instability in the vortex core due to its interaction with this negative vorticity region appears to be the catalyst for the magnification of the size and intensity of the turbulent flowfield downstream of the interaction. In general, the core radius increases while peak tangential velocity decreases with the effect being greater for smaller separations. These effects are largely independent of blade angle of attack; and if these parameters are normalized on their undisturbed values, then the effects of the vortex strength appear much weaker. Two theoretical models were developed to aid in extending the results to other flow conditions. An empirical model was developed for core parameter prediction which has some rudimentary physical basis, implying usefulness beyond a simple curve fit. An inviscid flow model was also created to estimate the vorticity shed by the interaction blade, and to predict the early stages of its incorporation into the interacting vortex.

Wittmer, Kenneth S.↗

NASA NTRS2013

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗

NASA NTRS2020

Spectral Range Within Global aCDOM(440) Algorithms for Oceanic, Coastal, and Inland Waters with Application to Airborne Measurements

The optically active component of dissolved organic material in aquatic ecosystems, or colored dissolved organic matter (CDOM), is represented by the coefficient of absorption due to the dissolved aquatic constituents at 440 nm, aCDOM(440). Remote sensing of aCDOM(440) enables characterization of ecosystem processes and aids in retrieval of chlorophyll a, a proxy for phytoplankton biomass. Spectrally adjacent band-ratio domains, e.g., blue to green, have previously been applied for remote sensing of aCDOM(440) in coastal and oceanic waters with similar results compared to more complex semi-analytical algorithms. Estimation of aCDOM(440) from ratios of the most spectrally separated ocean color wavebands (end members), e.g., ultraviolet (UV) to near-infrared (NIR), termed end-member analysis (EMA), has previously been shown to increase the accuracy of global aCDOM(440) retrievals from in-water observations of diffuse attenuation and to enable a unified algorithmic perspective without requiring regional adjustment of internal bio-optical parameters. EMA of above-water observations is evaluated herein, with a focus on coastal and inland waters in which increasing optical complexity and likelihood of bottom reflectance challenge the oceanic algorithms developed for deep and optically simple (case-1) waters. Analysis herein of three independent, in situ, bio-optical datasets indicates significant correlation between aCDOM(440) and end-member band ratios (next-generation 320 and 780 nm or legacy 412 and 670 nm ratios) with a coefficient of determination, R2, of 0.87 (log-scale) or higher based on a dataset spanning the dynamic range of global, conservative water bodies. For applicable wavelengths, EMA algorithms are shown to agree with case-1 relationships and to produce consistent log-scale uncertainties across more than three orders of magnitude in aCDOM(440) values (0.001–2.305 m􀀀1). EMA using UV and NIR wavelengths (320 and 780 nm) is applied to low-altitude airborne observations and satisfies 25% uncertainty based on unbiased percent differences (UPDs) within each of three dissimilar match-up sites ranging in aCDOM(440) from 0.02–0.57 m􀀀1. Results demonstrate that EMA is a useful and robust approach for the remote sensing of aCDOM(440) in coastal and inland waters, which are generally shallower, contain more optically complex environments, and span a greater range in aCDOM(440) than oceanic waters.

CDOM↗

DOE OSTI2018

SISTec: Mathematical Calibration of Large Clearance Monitors - 18265

One major path for clearance of large amount of waste arising from the decommission of nuclear power plants is based on measurement with large clearance monitors (also known as free release measurement facilities). These devices utilize multi-detector plastic scintillator arrays (large clearance monitors can consist of up to 24 plastic scintillators) to offer both, extremely low minimum detectable activities (MDA) and short measurement times for waste packages weighing up to one metric ton. However, a major problem with these devices is to obtain nuclide specific efficiencies, especially for geometrically complex large waste packages. Normal procedure for estimating the efficiency is to device dummy waste packages and attach point sources with well-known activities to unfavourable positions in these packages. In accordance with this procedure, calibration factors can be derived in a radiological conservative way. Devising such test dummies is not always possible and the procedure is prone to generating large measurement uncertainties. In this publication, we propose a new method to derive calibration factors by means of the Los Alamos National Laboratory (LANL) simulation tool Monte Carlo N-Particle (MCNP) and a statistical-phenomenological model for the detector characterization. In a multi-step approach, first the detector response for a distribution of activity in a waste package is calculated and subsequently the detectors are mathematically characterized. Applying state of the art statistical and optimizing procedures, parameters for the detector response are derived in a self-consistent way. This also allows to rigorously propagate the uncertainties to the final calculated calibration factors according to the Guide to the Expression of Uncertainty in Measurement (GUM supplement 1). To validate the derived parameters for the detector characterization a minimally invasive experimental setup was devised and installed into the large clearance monitor. It consisted of an aluminum swivel arm which could position point sources to the relevant coordinates. The so collected measurement data not only allowed to investigate the precision of the calculated calibration factors, but also to probe the location of concentrated activity sources in the waste package. Furthermore, data on measurements with dummy waste packages were collected to test the reproducibility of well-known calibration factors. We present the software tool SIStec, which will provide an intuitive user-interface to generate and simulate waste package models and execute the above-mentioned method. Using SIStec, different scenarios for the distribution of activity in the waste package can be probed. This will also allow to optimize the clearance process on various levels. For example, measurement time for pre-measurements can be reduced in exchange for unfavourable calibration factors, which would still allow to meet the regulator limits for clearance. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2018

Beyond the plane-parallel approximation for redshift surveys

Redshift -space distortions privilege the location of the observer in cosmological redshift surveys, breaking the translational symmetry of the underlying theory. This violation of statistical homogeneity has consequences for the modelling of clustering observables, leading to what are frequently called 'wide-angle effects'. We study these effects analytically, computing their signature in the clustering of the multipoles in configuration and Fourier space. We take into account both physical wide-angle contributions as well as the terms generated by the galaxy selection function. Similar considerations also affect the way power spectrum estimators are constructed. We quantify in an analytical way the biases that enter and clarify the relation between what we measure and the underlying theoretical modelling. The presence of an angular window function is also discussed. Motivated by this analysis, we present new estimators for the three dimensional Cartesian power spectrum and bispectrum multipoles written in terms of spherical Fourier-Bessel coefficients. In conclusion, we show how the latter have several interesting properties, allowing in particular a clear separation between angular and radial modes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Machine learning and algebraic approaches towards complete matter spectra in 4d F-theory

Motivated by engineering vector-like (Higgs) pairs in the spectrum of 4d F-theory compactifications, we combine machine learning and algebraic geometry techniques to analyze line bundle cohomologies on families of holomorphic curves. To quantify jumps of these cohomologies, we first generate 1.8 million pairs of line bundles and curves embedded in dP 3 , for which we compute the cohomologies. A white-box machine learning approach trained on this data provides intuition for jumps due to curve splittings, which we use to construct additional vector-like Higgs-pairs in an F-Theory toy model. We also find that, in order to explain quantitatively the full dataset, further tools from algebraic geometry, in particular Brill-Noether theory, are required. Using these ingredients, we introduce a diagrammatic way to express cohomology jumps across the parameter space of each family of matter curves, which reflects a stratification of the F-theory complex structure moduli space in terms of the vector-like spectrum. Furthermore, these insights provide an algorithmically efficient way to estimate the possible cohomology dimensions across the entire parameter space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Library based x-ray scatter correction for dedicated cone beam breast CT

Purpose: The image quality of dedicated cone beam breast CT (CBBCT) is limited by substantial scatter contamination, resulting in cupping artifacts and contrast-loss in reconstructed images. Such effects obscure the visibility of soft-tissue lesions and calcifications, which hinders breast cancer detection and diagnosis. In this work, we propose a library-based software approach to suppress scatter on CBBCT images with high efficiency, accuracy, and reliability. Methods: The authors precompute a scatter library on simplified breast models with different sizes using the GEANT4-based Monte Carlo (MC) toolkit. The breast is approximated as a semiellipsoid with homogeneous glandular/adipose tissue mixture. For scatter correction on real clinical data, the authors estimate the breast size from a first-pass breast CT reconstruction and then select the corresponding scatter distribution from the library. The selected scatter distribution from simplified breast models is spatially translated to match the projection data from the clinical scan and is subtracted from the measured projection for effective scatter correction. The method performance was evaluated using 15 sets of patient data, with a wide range of breast sizes representing about 95% of general population. Spatial nonuniformity (SNU) and contrast to signal deviation ratio (CDR) were used as metrics for evaluation. Results: Since the time-consuming MC simulation for library generation is precomputed, the authors’ method efficiently corrects for scatter with minimal processing time. Furthermore, the authors find that a scatter library on a simple breast model with only one input parameter, i.e., the breast diameter, sufficiently guarantees improvements in SNU and CDR. For the 15 clinical datasets, the authors’ method reduces the average SNU from 7.14% to 2.47% in coronal views and from 10.14% to 3.02% in sagittal views. On average, the CDR is improved by a factor of 1.49 in coronal views and 2.12 in sagittal views. Conclusions: The library-based scatter correction does not require increase in radiation dose or hardware modifications, and it improves over the existing methods on implementation simplicity and computational efficiency. As demonstrated through patient studies, the authors’ approach is effective and stable, and is therefore clinically attractive for CBBCT imaging.

60 APPLIED LIFE SCIENCES↗

NASA NTRSDate not supplied

Hourly Carbon Fluxes Estimation Using the GOES Advanced Baseline Imager (ABI) Data Over the Conterminous USA

Tremendous efforts by Fluxnet scientists over the past few decades have made thousands of site-years of carbon flux observations available for advancing our understanding of carbon cycling in terrestrial ecosystems. One of key Fluxnet measurements is net ecosystem exchange (NEE) as it is directly related to carbon budget of terrestrial ecosystems. However, carbon flux estimation studies using satellite remote sensing have focused mainly on daily Gross Primary Production (GPP). The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. Because daily NEE is close to zero value, the carbon flux models using the polar orbiting satellite data have not been well used for NEE estimation. The new generation of geostationary satellite sensors (e.g., GOES Advanced Baesline Imager (ABI) and Himawari Advanced Himawari Imager (AHI)) provide frequent observations, often less than every 10 minutes. Here, we use GOES ABI data to estimate hourly NEE over the conterminous USA. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and, in particular, solar radiation that is directly derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Fluxnet data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

Geostationary satellite↗

NASA NTRS2024

Real-Time Temperature Field State Observer and Digital Replica to Support Cyber-Physical Testing

Cyber-physical testing is a class of experimental methods in which a system is partitioned into physical and numerical components to study its overall behavior. Transfer systems are typically needed to capture relevant interactions by enforcing appropriate interface conditions between the physical and numerical portions. While cyber-physical testing is often performed using lumped-parameter systems, thermomechanical cyber-physical testing requires the use of thermal actuators as transfer systems to apply a thermal condition over a spatially distributed region or surface. Thus, the interactions between the numerical and physical components are of a distributed nature, i.e., a field rather than point values. In this paper, we develop and experimentally validate a state observer for temperature fields to enable thermomechanical cyber-physical testing. The state observer provides continuous estimates of temperature over a surface of the thermal actuator using only temperature measurements at discrete locations. The estimated temperature field is then leveraged to account for localized mechanical defects in the interface condition. The method consists of building a digital replica of one portion of the control plant. By imposing defects in the digital replica it is used to generate and output the response of a thermal load pattern that represents that of a physical system with mechanical defects.

Manuel Salmeron↗