Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set methods”

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

Backscattering and absorption coefficients for electrons: Solutions of invariant embedding transport equations using a method of convergence

The backscattering coefficient is a magnitude whose measurement is fundamental for the characterization of materials with techniques that make use of particle beams and particularly when performing microanalysis. In this work, we report the results of an analytic method to calculate the backscattering and absorption coefficients of electrons in similar conditions to those of electron probe microanalysis. Starting on a five level states ladder model in 3D, we deduced a set of integro-differential coupled equations of the coefficients with a method know as invariant embedding. By means of a procedure proposed by authors, called method of convergence, two types of approximate solutions for the set of equations, namely complete and simple solutions, can be obtained. Although the simple solutions were initially proposed as auxiliary forms to solve higher rank equations, they turned out to be also useful for the estimation of the aforementioned coefficients. In previous reports, we have presented results obtained with the complete solutions. In this paper, we present results obtained with the simple solutions of the coefficients, which exhibit a good degree of fit with the experimental data. Both the model and the calculation method presented here can be generalized to other techniques that make use of different sorts of particle beams.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Longitudinal, intermodality registration of quantitative breast PET and MRI data acquired before and during neoadjuvant chemotherapy: Preliminary results

Purpose: The authors propose a method whereby serially acquired DCE-MRI, DW-MRI, and FDG-PET breast data sets can be spatially and temporally coregistered to enable the comparison of changes in parameter maps at the voxel level. Methods: First, the authors aligned the PET and MR images at each time point rigidly and nonrigidly. To register the MR images longitudinally, the authors extended a nonrigid registration algorithm by including a tumor volume-preserving constraint in the cost function. After the PET images were aligned to the MR images at each time point, the authors then used the transformation obtained from the longitudinal registration of the MRI volumes to register the PET images longitudinally. The authors tested this approach on ten breast cancer patients by calculating a modified Dice similarity of tumor size between the PET and MR images as well as the bending energy and changes in the tumor volume after the application of the registration algorithm. Results: The median of the modified Dice in the registered PET and DCE-MRI data was 0.92. For the longitudinal registration, the median tumor volume change was −0.03% for the constrained algorithm, compared to −32.16% for the unconstrained registration algorithms (p = 8 × 10{sup −6}). The medians of the bending energy were 0.0092 and 0.0001 for the unconstrained and constrained algorithms, respectively (p = 2.84 × 10{sup −7}). Conclusions: The results indicate that the proposed method can accurately spatially align DCE-MRI, DW-MRI, and FDG-PET breast images acquired at different time points during therapy while preventing the tumor from being substantially distorted or compressed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2014

Explicit correlation treatment of the potential energy surface of CO{sub 2} dimer

We present an extensive study of the four-dimensional potential energy surface (4D-PES) of the carbon dioxide dimer, (CO{sub 2}){sub 2}. This PES is developed over the set of intermolecular coordinates. The electronic computations are carried out at the explicitly correlated coupled cluster method with single, double, and perturbative triple excitations [CCSD(T)-F12] level of theory in connection with the augmented correlation-consistent aug-cc-pVTZ basis set. An analytic representation of the 4D-PES is derived. Our extensive calculations confirm that “Slipped Parallel” is the most stable form and that the T-shaped structure corresponds to a transition state. Later on, this PES is employed for the calculations of the vibrational energy levels of the dimer. Moreover, the temperature dependence of the dimer second virial coefficient and of the first spectral moment of rototranslational collision-induced absorption spectrum is derived. For both quantities, a good agreement is found between our values and the experimental data for a wide range of temperatures. This attests to the high quality of our PES. Generally, our PES and results can be used for modeling CO{sub 2} supercritical fluidity and examination of its role in planetary atmospheres. It can be also incorporated into dynamical computations of CO{sub 2} capture and sequestration. This allows deep understanding, at the microscopic level, of these processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

A multilevel method for structural synthesis

A multilevel method for structural synthesis based on finite element modeling at both the system and the substructure level is set forth. At the system level, the design is described using basis vectors generated via substructure design modification. Whenever possible, substructure constraints are satisfied via local material redistribution. At both the system level and the substructure level, design modification is carried out by solving a sequence of approximate problems. The method generates upper bound approximations for the optimum design and it has potential for substantial gains in computational efficiency, because the uncoupled substructure design modification problems can be processed in parallel.

Schmit, L. A.↗
DOE OSTI2021

Characterization and Optimization of Building Blocks for Specialized Computing Platforms

As noise limits the performance of quantum processors, the ability to characterize this noise and develop methods to overcome it is essential for the future of quantum computing. In this report, we develop a complete set of tools for improving quantum processor performance at the application level, including low-level physical models of quantum gates, a numerically efficient method of producing process matrices that span a wide range of model parameters, and full-channel quantum simulations. We then provide a few examples of how to use these tools to study the effects of noise on quantum circuits.

42 ENGINEERING↗
DOE OSTI2024

Machine Learning Classification Strategy to Improve Streamflow Estimates in Diverse River Basins in the Colorado River Basin

Streamflow in the Colorado River Basin (CRB) is significantly altered by human activities including land use/cover alterations, reservoir operation, irrigation, and water exports. Climate is also highly varied across the CRB which contains snowpack-dominated watersheds and arid, precipitation-dominated basins. Recently, machine learning methods have improved the generalizability and accuracy of streamflow models. Previous successes with LSTM modeling have primarily focused on unimpacted basins, and few studies have included human impacted systems in either regional or single-basin modeling. We demonstrate that the diverse hydrological behavior of river basins in the CRB are too difficult to model with a single, regional model. We propose a method to delineate catchments into categories based on the level of predictability, hydrological characteristics, and the level of human influence. Lastly, we model streamflow in each category with climate and anthropogenic proxy data sets and use feature importance methods to assess whether model performance improves with additional relevant data. Overall, land use cover data at a low temporal resolution was not sufficient to capture the irregular patterns of reservoir releases, demonstrating the importance of having high-resolution reservoir release data sets at a global scale. On the other hand, the classification approach reduced the complexity of the data and has the potential to improve streamflow forecasts in human-altered regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

EPRI Research Results: International Nuclear Power Plant Waste Classification - 17444

Different countries around the globe use varying schemes for the classification of radioactive waste. Many countries develop country specific waste classification schemes based on the International Atomic Energy Agency (IAEA) waste classification framework. The United States of America use the Class A, B, and C waste classification. Additionally, some countries use concentration limits to classify wastes while others use dose rate limits or a combination of concentration and dose rate limits. General comparisons have been made between the IAEA and the United States classification schemes, but a clear and well researched comparison of the international schemes was not evident in the available research data sets. In 2013 Electric Power Research Institute (EPRI) initiated the Comparison of Global Low and Intermediate Level Waste Management Methods project. The objective of the project was to understand how six countries managed their low and intermediate level wastes from generation to disposal. A major task of this project was to conduct a comparison of the international waste classification schemes. EPRI developed a set of reference waste volumes that are representative of nuclear power plant low-level radioactive waste in form and radioactivity concentrations (e.g. dry active waste, high activity resin, low activity resin, high activity filters and low activity filters). These waste volumes were then classified using the waste classification systems in place in several countries around the world. In 2016, EPRI presented a paper providing a status update on this research. This 2017 paper provides the final results of that research, show where comparisons and correlations can be made, and provide a discussion on how this information might inform international radioactive research and development and benchmarking of radioactive waste management strategies, operational experiences, best practices, and technology implementation. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2019

New York City Panel on Climate Change 2019 Report Chapter 4: Coastal Flooding

Coastal flooding from storm surge is one of the most dangerous and damaging natural hazards that societies face. It was responsible for half of all hurricane related mortalities in the United States from 1963 to 2012, far more than any other factor. Coastal extreme water levels are increasing globally, mainly driven by rises in mean sea level. Sea level rise is also causing rapid increases in the annual number of shallow “nuisance floods” for low-lying neighborhoods. The objectives of this chapter are to review the latest knowledge on New York City flood risk from storms and tides, and to evaluate how climate change will affect this risk between now and the end of the century. Methods used by NPCC (2015) for assessing storm-driven extreme floods are generally repeated here, including the use of the Federal Emergency Management Agency (FEMA, 2013) baseline flood hazards (e.g., the 100-year flood(sup a)) and the methods for adding sea level rise and mapping the resulting hazard. New advancements include an innovative analysis of monthly tidal flooding based on a dynamic model, a broadened set of sea level rise scenarios supplemented with the Antarctic Rapid Ice Melt (ARIM) scenario (see Chapter 3), and sensitivity analyses that show how differing methods would affect our results. Wind is a primary factor for coastal storm surge, and a brief review is given in Appendix 4.A, with the latest scientific knowledge on what drives extreme wind events in the New York City area and how they may change in the future.

Orton, Philip↗
NASA NTRS2013

Uncertainties in Atomic Data and Their Propagation Through Spectral Models. I.

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data.We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of Oiii and Fe ii and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe ii]. Key words: atomic data - atomic processes - line: formation - methods: data analysis - molecular data - molecular processes - techniques: spectroscopic

linear sets↗
DOE OSTI2014

Charge transfer and atomic-level pressure in metallic glasses

This paper presents a systematic study on the charge transfer and ionicity in various metallic-glass forming systems, as well as its relationship with other atomic-level structure indicators, using the Bader analysis method and molecular dynamics simulation. It is shown that in a binary or multicomponent system, the chemical effects (when more than one elements present) appear to play a more important role in setting the absolute level of the atomic-level pressure, compared to the topological fluctuation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

A multi-level solution algorithm for steady-state Markov chains

A new iterative algorithm, the multi-level algorithm, for the numerical solution of steady state Markov chains is presented. The method utilizes a set of recursively coarsened representations of the original system to achieve accelerated convergence. It is motivated by multigrid methods, which are widely used for fast solution of partial differential equations. Initial results of numerical experiments are reported, showing significant reductions in computation time, often an order of magnitude or more, relative to the Gauss-Seidel and optimal SOR algorithms for a variety of test problems. The multi-level method is compared and contrasted with the iterative aggregation-disaggregation algorithm of Takahashi.

Horton, Graham↗
DOE OSTI2013

UNCERTAINTIES IN ATOMIC DATA AND THEIR PROPAGATION THROUGH SPECTRAL MODELS. I

We present a method for computing uncertainties in spectral models, i.e., level populations, line emissivities, and emission line ratios, based upon the propagation of uncertainties originating from atomic data. We provide analytic expressions, in the form of linear sets of algebraic equations, for the coupled uncertainties among all levels. These equations can be solved efficiently for any set of physical conditions and uncertainties in the atomic data. We illustrate our method applied to spectral models of O III and Fe II and discuss the impact of the uncertainties on atomic systems under different physical conditions. As to intrinsic uncertainties in theoretical atomic data, we propose that these uncertainties can be estimated from the dispersion in the results from various independent calculations. This technique provides excellent results for the uncertainties in A-values of forbidden transitions in [Fe II].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

TU-F-CAMPUS-I-05: Parameterization of the Noise Power Spectrum in X-Ray Computed Tomography

Purpose: The purpose of this work was to develop a method so that the noise power spectrum (NPS) can be approximated for arbitrary levels of mAs, from a single determination in CT. Methods: The NPS is factorized into 2 components, 1) a parameterized function representing the 1D normalized spatial frequency distribution and 2) a function to scale the magnitude of 1) for arbitrary values of mAs. The 1D NPS, normalized by image variance (NNPS), was determined for 2 FBP reconstruction kernels (smoothing and edge enhancing) for 400 mAs. The NNPS were fit to the parameterized function and a scaling function was established to approximate the image variance at arbitrary values of mAs. Using the root mean square error normalized by the maximum value (NRMSE), the NPS approximated with the factorization method was compared to the NPS determined at 5 different mAs levels. Results: The factorization resulted in a set of 7 coefficients that can be used to approximate the 1D NPS, for arbitrary levels of mAs, for the convolution kernels studied in this work. The approximated NPS (factorization) agreed well with the determined NPS for all mAs levels. The greatest NRMSE was 0.02 and was observed for the edge enhancing kernel. Conclusion: The proposed factorization method has been demonstrated as applicable for FBP reconstruction. It can be used to approximate the 1D NPS for arbitrary levels of mAs, from a single NPS determination. Furthermore, approximations of the 1D NPS can conveniently be distributed since the factorization method only used 7 coefficients in the approximation.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1983

The special effort processing of FGGE data

The basic FGGE level IIb data set was enhanced. It focused on removing deficiencies in the objective methods of quality assurance, removing efficiencies in certain types of operationally produced satellite soundings, and removing deficiencies in certain types of operationally produced cloud tracked winds. The Special Effort was a joint NASA-NOAA-University of Wisconsin effort. The University of Wisconsin installed an interactive McIDAS capability on the Amdahl computer at the Goddard Laboratory of Atmospheric Sciences (GLAS) with one interactive video terminal at Goddard and the other at the World Weather Building. With this interactive capability a joint processing effort was undertaken to reprocess certain FGGE data sets. NOAA produced a specially edited data set for the special observing periods (SOPs) of FGGE. NASA produced an enhanced satellite sounding data set for the SOPs while the University of Wisconsin produced an enhanced cloud tracked wind set from the Japanese geostationary satellite images.

Source record↗
NASA NTRS1989

The determination of accurate dipole polarizabilities alpha and gamma for the noble gases

The static dipole polarizabilities alpha and gamma for the noble gases helium through xenon were determined using large flexible one-particle basis sets in conjunction with high-level treatments of electron correlation. The electron correlation methods include single and double excitation coupled-cluster theory (CCSD), an extension of CCSD that includes a perturbational estimate of connected triple excitations, CCSD(T), and second order perturbation theory (MP2). The computed alpha and gamma values are estimated to be accurate to within a few percent. Agreement with experimental data for the static hyperpolarizability gamma is good for neon and xenon, but for argon and krypton the differences are larger than the combined theoretical and experimental uncertainties. Based on our calculations, we suggest that the experimental value of gamma for argon is too low; adjusting this value would bring the experimental value of gamma for krypton into better agreement with our computed result. The MP2 values for the polarizabilities of neon, argon, krypton and zenon are in reasonabe agreement with the CCSD and CCSD(T) values, suggesting that this less expensive method may be useful in studies of polarizabilities for larger systems.

Rice, Julia E.↗
DOE OSTI2013

Texture analysis improves level set segmentation of the anterior abdominal wall

Purpose: The treatment of ventral hernias (VH) has been a challenging problem for medical care. Repair of these hernias is fraught with failure; recurrence rates ranging from 24% to 43% have been reported, even with the use of biocompatible mesh. Currently, computed tomography (CT) is used to guide intervention through expert, but qualitative, clinical judgments, notably, quantitative metrics based on image-processing are not used. The authors propose that image segmentation methods to capture the three-dimensional structure of the abdominal wall and its abnormalities will provide a foundation on which to measure geometric properties of hernias and surrounding tissues and, therefore, to optimize intervention.Methods: In this study with 20 clinically acquired CT scans on postoperative patients, the authors demonstrated a novel approach to geometric classification of the abdominal. The authors’ approach uses a texture analysis based on Gabor filters to extract feature vectors and follows a fuzzy c-means clustering method to estimate voxelwise probability memberships for eight clusters. The memberships estimated from the texture analysis are helpful to identify anatomical structures with inhomogeneous intensities. The membership was used to guide the level set evolution, as well as to derive an initial start close to the abdominal wall.Results: Segmentation results on abdominal walls were both quantitatively and qualitatively validated with surface errors based on manually labeled ground truth. Using texture, mean surface errors for the outer surface of the abdominal wall were less than 2 mm, with 91% of the outer surface less than 5 mm away from the manual tracings; errors were significantly greater (2–5 mm) for methods that did not use the texture.Conclusions: The authors’ approach establishes a baseline for characterizing the abdominal wall for improving VH care. Inherent texture patterns in CT scans are helpful to the tissue classification, and texture analysis can improve the level set segmentation around the abdominal region.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Accurate and Efficient ab Initio Calculations for Supramolecular Complexes: Symmetry-Adapted Perturbation Theory with Many-Body Dispersion

Symmetry-adapted perturbation theory (SAPT) provides a chemically meaningful energy decomposition scheme for nonbonded interactions that is useful for interpretive purposes. Although formally a dimer theory, we have previously introduced an “extended” version (XSAPT) that incorporates many-body polarization via self-consistent charge embedding. Here, we extend the XSAPT methodology to include nonadditive dispersion, using a modified form of the many-body dispersion (MBD) method of Tkatchenko and co-workers. Dispersion interactions beyond the pairwise atom–atom approximation improve total interaction energies even in small systems, and for large π-stacked complexes these corrections can amount to several kilocalories per mole. The XSAPT+MBD method introduced here achieves errors of ≲1 kcal/mol (as compared to high-level ab initio benchmarks) for the L7 data set of large dispersion-bound complexes and ≲4 kcal/mol (as compared to experiment) for the S30L data set of host–guest complexes. This is superior to the best contemporary density functional methods for noncovalent interactions, at comparable or lower cost. XSAPT+MBD represents a promising method for application to supramolecular assemblies, including protein–ligand binding.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2012

Direct dynamics simulation of dioxetane formation and decomposition via the singlet ·O–O–CH2–CH2 · biradical: Non-RRKM dynamics

Electronic structure calculations and direct chemical dynamics simulations are used to study the formation and decomposition of dioxetane on its ground state singlet potential energy surface. The stationary points for 1O2 + C2H4, the singlet {center_dot}O-O-CH2-CH2{center_dot} biradical, the transition state (TS) connecting this biradical with dioxetane, and the two transition states and gauche {center_dot}O-CH2-CH2-O{center_dot} biradical connecting dioxetane with the formaldehyde product molecules are investigated at different levels of electronic structure theory including UB3LYP, UMP2, MRMP2, and CASSCF and a range of basis sets. The UB3LYP/6-31G* method was found to give representative energies for the reactive system and was used as a model for the simulations. UB3LYP/6-31G* direct dynamics trajectories were initiated at the TS connecting the {center_dot}O-O-CH2-CH2{center_dot} biradical and dioxetane by sampling the TS's vibrational energy levels, and rotational and reaction coordinate energies, with Boltzmann distributions at 300, 1000, and 1500 K. This corresponds to the transition state theory model for trajectories that pass the TS. The trajectories were directed randomly towards both the biradical and dioxetane. A small fraction of the trajectories directed towards the biradical recrossed the TS and formed dioxetane. The remainder formed 1O2 + C2H4 and of these {approx} 40% went directly from the TS to 1O2 + C2H4 without getting trapped and forming an intermediate in the {center_dot}O-O- CH2-CH2{center_dot} biradical potential energy minimum, a non-statistical result. The dioxetane molecules which are formed dissociate to two formaldehyde molecules with a rate constant two orders of magnitude smaller than that predicted by Rice-Ramsperger-Kassel-Marcus theory. The reaction dynamics from dioxetane to the formaldehyde molecules do not follow the intrinsic reaction coordinate or involve trapping in the gauche {center_dot}O-CH2-CH2-O{center_dot} biradical potential energy minimum. Important non-statistical dynamics are exhibited for this reactive system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗