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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Refined energetic ordering for sulphate–water ( n = 3–6) clusters using high-level electronic structure calculations

This work reports refinements of the energetic ordering of the known low-energy structures of sulfate-water clusters SO 4 2- (H2O)n (n = 3-6) using high-level electronic structure methods. Coupled cluster singles and doubles with perturbative triples (CCSD(T)) is used in combination with an estimate of basis set effects up to the complete basis set limit using second order Møller-Plesset theory. Harmonic zero point energy (ZPE), included at the B3LYP/6-311++G** level, was found to have a significant effect on the energetic ordering. Limitations of the ZPE calculations, both due to electronic structure errors, and use of the harmonic approximation, probably constitute the largest remaining errors. Due to the often small energy differences between cluster isomers, and the significant role of ZPE, deuteration can alter the relative energies of low-lying structures, and, when it is applied in conjunction with calculated harmonic ZPE’s, even alters the global minimum for n = 4.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Understanding density functional theory (DFT) and completing it in practice

We review some salient points in the derivation of density functional theory (DFT) and of the local density approximation (LDA) of it. We then articulate an understanding of DFT and LDA that seems to be ignored in the literature. We note the well-established failures of many DFT and LDA calculations to reproduce the measured energy gaps of finite systems and band gaps of semiconductors and insulators. We then illustrate significant differences between the results from self consistent calculations using single trial basis sets and those from computations following the Bagayoko, Zhao, and Williams (BZW) method, as enhanced by Ekuma and Franklin (BZW-EF). Unlike the former, the latter calculations verifiably attain the absolute minima of the occupied energies, as required by DFT. These minima are one of the reasons for the agreement between their results and corresponding, experimental ones for the band gap and a host of other properties. Further, we note predictions of DFT BZW-EF calculations that have been confirmed by experiment. Our subsequent description of the BZW-EF method ends with the application of the Rayleigh theorem in the selection, among the several calculations the method requires, of the one whose results have a full, physics content ascribed to DFT. This application of the Rayleigh theorem adds to or completes DFT, in practice, to preserve the physical content of unoccupied, low energy levels. Discussions, including implications of the method, and a short conclusion follow the description of the method. The successive augmentation of the basis set in the BZW-EF method, needed for the application of the Rayleigh theorem, is also necessary in the search for the absolute minima of the occupied energies, in practice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Use of 16S rRNA gene for identification of a broad range of clinically relevant bacterial pathogens

According to World Health Organization statistics of 2011, infectious diseases remain in the top five causes of mortality worldwide. However, despite sophisticated research tools for microbial detection, rapid and accurate molecular diagnostics for identification of infection in humans have not been extensively adopted. Time-consuming culture-based methods remain to the forefront of clinical microbial detection. The 16S rRNA gene, a molecular marker for identification of bacterial species, is ubiquitous to members of this domain and, thanks to ever-expanding databases of sequence information, a useful tool for bacterial identification. In this study, we assembled an extensive repository of clinical isolates (n = 617), representing 30 medically important pathogenic species and originally identified using traditional culture-based or non-16S molecular methods. This strain repository was used to systematically evaluate the ability of 16S rRNA for species level identification. To enable the most accurate species level classification based on the paucity of sequence data accumulated in public databases, we built a Naïve Bayes classifier representing a diverse set of high-quality sequences from medically important bacterial organisms. We show that for species identification, a model-based approach is superior to an alignment based method. Overall, between 16S gene based and clinical identities, our study shows a genus-level concordance rate of 96% and a species-level concordance rate of 87.5%. We point to multiple cases of probable clinical misidentification with traditional culture based identification across a wide range of gram-negative rods and gram-positive cocci as well as common gram-negative cocci.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Optimal Wind Power Uncertainty Intervals for Electricity Market Operation

It is important to select an appropriate uncertainty level of the wind power forecast for power system scheduling and electricity market operation. Traditional methods hedge against a predefined level of wind power uncertainty, such as a specific confidence interval or uncertainty set, which leaves the questions of how to best select the appropriate uncertainty levels. To bridge this gap, this paper proposes a model to optimize the forecast uncertainty intervals of wind power for power system scheduling problems, with the aim of achieving the best trade-off between economics and reliability. Then we reformulate and linearize the models into a mixed integer linear programming (MILP) without strong assumptions on the shape of the probability distribution. In order to invest the impacts on cost, reliability, and prices in a electricity market, we apply the proposed model on a twosettlement electricity market based on a six-bus test system and on a power system representing the U.S. state of Illinois. The results show that the proposed method can not only help to balance the economics and reliability of the power system scheduling, but also help to stabilize the energy prices in electricity market operation.

electricity market↗
DOE OSTI2016

TH-AB-202-08: A Robust Real-Time Surface Reconstruction Method On Point Clouds Captured From a 3D Surface Photogrammetry System

Purpose: Surface photogrammetry (e.g. VisionRT, C-Rad) provides a noninvasive way to obtain high-frequency measurement for patient motion monitoring in radiotherapy. This work aims to develop a real-time surface reconstruction method on the acquired point clouds, whose acquisitions are subject to noise and missing measurements. In contrast to existing surface reconstruction methods that are usually computationally expensive, the proposed method reconstructs continuous surfaces with comparable accuracy in real-time. Methods: The key idea in our method is to solve and propagate a sparse linear relationship from the point cloud (measurement) manifold to the surface (reconstruction) manifold, taking advantage of the similarity in local geometric topology in both manifolds. With consistent point cloud acquisition, we propose a sparse regression (SR) model to directly approximate the target point cloud as a sparse linear combination from the training set, building the point correspondences by the iterative closest point (ICP) method. To accommodate changing noise levels and/or presence of inconsistent occlusions, we further propose a modified sparse regression (MSR) model to account for the large and sparse error built by ICP, with a Laplacian prior. We evaluated our method on both clinical acquired point clouds under consistent conditions and simulated point clouds with inconsistent occlusions. The reconstruction accuracy was evaluated w.r.t. root-mean-squared-error, by comparing the reconstructed surfaces against those from the variational reconstruction method. Results: On clinical point clouds, both the SR and MSR models achieved sub-millimeter accuracy, with mean reconstruction time reduced from 82.23 seconds to 0.52 seconds and 0.94 seconds, respectively. On simulated point cloud with inconsistent occlusions, the MSR model has demonstrated its advantage in achieving consistent performance despite the introduced occlusions. Conclusion: We have developed a real-time and robust surface reconstruction method on point clouds acquired by photogrammetry systems. It serves an important enabling step for real-time motion tracking in radiotherapy. This work is supported in part by NIH grant R01 CA169102-02.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1975

On the determination of haze levels from Landsat data

The paper describes two methods for determining haze levels (specified by haze optical depth at a wavelength of 0.5 microns) from Landsat multispectral scanner data. The channel correlation method relates the haze level to the y-intercept of the regression line through a plot of the data in the plane of the multispectral channels MSS 4 and MSS 5. The minimum value method relates haze level to the minimum value of individual lines in the MSS 4 data set.

Potter, J. F.↗
NASA NTRS1986

Transition and mixing in axisymmetric jets and vortex rings

A class of impulsively started, axisymmetric, laminar jets produced by a time dependent joint source of momentum are considered. These jets are different flows, each initially at rest in an unbounded fluid. The study is conducted at three levels of detail. First, a generalized set of analytic creeping flow solutions are derived with a method of flow classification. Second, from this set, three specific creeping flow solutions are studied in detail: the vortex ring, the round jet, and the ramp jet. This study involves derivation of vorticity, stream function, entrainment diagrams, and evolution of time lines through computer animation. From entrainment diagrams, critical points are derived and analyzed. The flow geometry is dictated by the properties and location of critical points which undergo bifurcation and topological transformation (a form of transition) with changing Reynolds number. Transition Reynolds numbers were calculated. A state space trajectory was derived describing the topological behavior of these critical points. This state space derivation yielded three states of motion which are universal for all axisymmetric jets. Third, the axisymmetric round jet is solved numerically using the unsteady laminar Navier Stokes equations. These equations were shown to be self similar for the round jet. Numerical calculations were performed up to a Reynolds number of 30 for a 60x60 point mesh. Animations generated from numerical solution showed each of the three states of motion for the round jet, including the Re = 30 case.

Allen, G. A., Jr.↗
NASA NTRS1991

Compiling global name-space parallel loops for distributed execution

Distributed memory machines do not provide hardware support for a global address space. Thus programmers are forced to partition the data across the memories of the architecture and use explicit message passing to communicate data between processors. The compiler support required to allow programmers to express their algorithms using a global name-space is examined. A general method is presented for analysis of a high level source program and its translation into a set of independently executing tasks communicating via messages. If the compiler has enough information, this translation can be carried out at compile time. Otherwise, run-time code is generated to implement the required data movement. The analysis required in both situations is described and the performance of the generated code on the Intel iPSC/2 is presented.

Koelbel, Charles↗
NASA NTRS2012

Spacecraft Line-of-Sight Stabilization Using LWIR Earth Signature

The objective of this study is to investigate the potential of using the bright and near-uniform Earth infrared (or wavelength infrared, LWIR) signature as a stable reference for accurate (micro-rad or less) inertial pointing and tracking on-board an space vehicle, including the determination of the fundamental limits of applicability of the proposed method for space missions. We demonstrate sub-micro radian level pointing accuracy under a representative set of disturbances experienced by the spacecraft in orbit.

infrared detector↗
NASA NTRS2018

Improved Pyrotechnic Shock Scaling Method Based on Shock Response Spectrum and Statistical Energy Analysis

It has been more than forty years since the pyrotechnic shock scaling method was introduced. The scaling method estimates the attenuation of the Shock Response Spectrum (SRS) based on the distance from the source, structural configurations, types of structural joints and interfaces, and intervening structure [1]. The method has been successfully used in the spacecraft community countless times and still is frequently used to develop pyrotechnic shock requirements at various levels of assembly of a spacecraft. However, since the method was derived empirically from a limited set of shock test data [2], the aerospace community has been looking for an alternative approach. In this paper, a computational scaling method based on Statistical Energy Analysis (SEA) is re-visited. SEA is traditionally used as a method for investigating the diffusion of acoustic and vibratory energies of a system at steady-state [3]. Previously, an approach using the SRS as an acceleration constraint condition to the SEA model was introduced to estimate the attenuation from a shock source [4,5]. In the current investigation, additional examples are provided to further validate the approach. In addition, responses in the time domain are produced by the Local Modal Phase Reconstruction (LMPR) approach which was recently introduced to the community [6].

Fernandez, Juan P.↗
DOE OSTI2015

Multiscale simulations of hexa-peri-hexabenzocoronene and hexa-n-dodecyl-hexa-peri-hexabenzocoronene

This study concerns atomistic and coarse-grained molecular dynamics simulations of two disk-shaped molecules in the crystalline phase. The coarse-grained models were developed by applying the Iterative Boltzmann Inversion method which is a systematic coarse-graining method. In particular, a set of radial distribution functions and intramolecular distributions are reproduced at the coarse-grained level. Before applying coarse-graining, a reliable atomistic model was developed to reproduce the main experimental properties of these molecules. The crystalline phases are analyzed in terms of the Saupe ordering tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A Sensitivity-informed, Process Model Abstraction Implementation in GoldSim - 19276

System-level Performance Assessment (PA) models rely on specific process-level models for some inputs, such as infiltration. For probabilistic PA, the value of the process model output depends on the values of input parameters that also vary in the PA model. The simplest conceptual solution is to draw values of the input parameters from their distributions, run the process level model at these values, and then use the output where it is needed in the PA model. However, this solution is unsatisfactory for multiple reasons: (1) the PA Model should be set up such that new runs can be made and distributions can change without having to go back and run the process-level model; (2) calling the process model from the system model is possible, but technical difficulties remain in coupling the two; and (3) there is no uncertainty in the process-level model results incorporated into the PA model. Therefore, it is desirable to develop a method that maps the values of PA Model inputs to process level model outputs, in real time within the PA Model set up and using as simple of a map as possible without losing key features of the relationship. This has been termed 'model abstraction' because information about process model runs is 'abstracted' to be used with the PA model. A simple 'nearest neighbor' abstraction approach is presented to overcome these problems, in which a sensitivity-informed distance metric is used to map a system-level model realization to the most appropriate process-level model realization. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

A new ab initio potential energy surface and infrared spectra for the Ar–CS{sub 2} complex

We report a new three-dimensional potential energy surface for Ar–CS{sub 2} involving the Q{sub 3} normal mode for the υ{sub 3} antisymmetric stretching vibration of the CS{sub 2} molecule. The potential energies were calculated using the supermolecular method at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples, using augmented correlation-consistent quadruple-zeta basis set plus midpoint bond functions. Two vibrationally averaged potentials with CS{sub 2} at both the ground (υ = 0) and the first excited (υ = 1)υ{sub 3} vibrational states were generated from the integration of the three-dimensional potential over the Q{sub 3} coordinate. Each potential was found to have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation /angular finite basis representation method and the Lanczos algorithm were applied to calculate the rovibrational energy levels. The calculated band origin shift of the complex (0.0622 cm{sup −1}) is very close to the observed one (0.0671 cm{sup −1}). The predicted infrared spectra and spectroscopic parameters based on the two averaged potentials are in excellent agreement with the available experimental data.

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

Controlling flexible structures: A survey of methods

Most of the presently available control system design techniques applicable to flexible structure problems were developed to design controllers for rigid body systems. Although many of these design methods can be applied to flexible dynamics problems, recently developed techniques may be more suitable for flexible structure controller design. The purpose of this presentation is to examine briefly the peculiarities of the dynamics of flexible structures and to stimulate discussion about top level controller design approaches when designing controllers for flexible structures. Presented here is a suggestion of a set of categories of design methods for designing controllers for flexible structures as well as a discussion of the advantages and disadvantages of each category. No attempt has been made herein to select one category of design techniques as the best for flexible structure controller design. Instead, it is hoped that the structure suggested by these categories will facilitate further discussion on the merits of particular methods that will eventually point to those design techniques suitable for further development.

Benson, Russell A.↗
NASA NTRS2004

Robust Requirements Tracing via Internet Search Technology: Improving an IV and V Technique

There are three major objectives to this phase of the work. (1) Improvement of Information Retrieval (IR) methods for Independent Verification and Validation (IV&V) requirements tracing. Information Retrieval methods are typically developed for very large (order of millions - tens of millions and more documents) document collections and therefore, most successfully used methods somewhat sacrifice precision and recall in order to achieve efficiency. At the same time typical IR systems treat all user queries as independent of each other and assume that relevance of documents to queries is subjective for each user. The IV&V requirements tracing problem has a much smaller data set to operate on, even for large software development projects; the set of queries is predetermined by the high-level specification document and individual requirements considered as query input to IR methods are not necessarily independent from each other. Namely, knowledge about the links for one requirement may be helpful in determining the links of another requirement. Finally, while the final decision on the exact form of the traceability matrix still belongs to the IV&V analyst, his/her decisions are much less arbitrary than those of an Internet search engine user. All this suggests that the information available to us in the framework of the IV&V tracing problem can be successfully leveraged to enhance standard IR techniques, which in turn would lead to increased recall and precision. We developed several new methods during Phase II; (2) IV&V requirements tracing IR toolkit. Based on the methods developed in Phase I and their improvements developed in Phase II, we built a toolkit of IR methods for IV&V requirements tracing. The toolkit has been integrated, at the data level, with SAIC's SuperTracePlus (STP) tool; (3) Toolkit testing. We tested the methods included in the IV&V requirements tracing IR toolkit on a number of projects.

Hayes, Jane↗
NASA NTRS1980

Prediction of noise constrained optimum takeoff procedures

An optimization method is used to predict safe, maximum-performance takeoff procedures which satisfy noise constraints at multiple observer locations. The takeoff flight is represented by two-degree-of-freedom dynamical equations with aircraft angle-of-attack and engine power setting as control functions. The engine thrust, mass flow and noise source parameters are assumed to be given functions of the engine power setting and aircraft Mach number. Effective Perceived Noise Levels at the observers are treated as functionals of the control functions. The method is demonstrated by applying it to an Advanced Supersonic Transport aircraft design. The results indicate that automated takeoff procedures (continuously varying controls) can be used to significantly reduce community and certification noise without jeopardizing safety or degrading performance.

Padula, S. L.↗
DOE OSTI2015

Toward accurate tooth segmentation from computed tomography images using a hybrid level set model

Purpose: A three-dimensional (3D) model of the teeth provides important information for orthodontic diagnosis and treatment planning. Tooth segmentation is an essential step in generating the 3D digital model from computed tomography (CT) images. The aim of this study is to develop an accurate and efficient tooth segmentation method from CT images. Methods: The 3D dental CT volumetric images are segmented slice by slice in a two-dimensional (2D) transverse plane. The 2D segmentation is composed of a manual initialization step and an automatic slice by slice segmentation step. In the manual initialization step, the user manually picks a starting slice and selects a seed point for each tooth in this slice. In the automatic slice segmentation step, a developed hybrid level set model is applied to segment tooth contours from each slice. Tooth contour propagation strategy is employed to initialize the level set function automatically. Cone beam CT (CBCT) images of two subjects were used to tune the parameters. Images of 16 additional subjects were used to validate the performance of the method. Volume overlap metrics and surface distance metrics were adopted to assess the segmentation accuracy quantitatively. The volume overlap metrics were volume difference (VD, mm{sup 3}) and Dice similarity coefficient (DSC, %). The surface distance metrics were average symmetric surface distance (ASSD, mm), RMS (root mean square) symmetric surface distance (RMSSSD, mm), and maximum symmetric surface distance (MSSD, mm). Computation time was recorded to assess the efficiency. The performance of the proposed method has been compared with two state-of-the-art methods. Results: For the tested CBCT images, the VD, DSC, ASSD, RMSSSD, and MSSD for the incisor were 38.16 ± 12.94 mm{sup 3}, 88.82 ± 2.14%, 0.29 ± 0.03 mm, 0.32 ± 0.08 mm, and 1.25 ± 0.58 mm, respectively; the VD, DSC, ASSD, RMSSSD, and MSSD for the canine were 49.12 ± 9.33 mm{sup 3}, 91.57 ± 0.82%, 0.27 ± 0.02 mm, 0.28 ± 0.03 mm, and 1.06 ± 0.40 mm, respectively; the VD, DSC, ASSD, RMSSSD, and MSSD for the premolar were 37.95 ± 10.13 mm{sup 3}, 92.45 ± 2.29%, 0.29 ± 0.06 mm, 0.33 ± 0.10 mm, and 1.28 ± 0.72 mm, respectively; the VD, DSC, ASSD, RMSSSD, and MSSD for the molar were 52.38 ± 17.27 mm{sup 3}, 94.12 ± 1.38%, 0.30 ± 0.08 mm, 0.35 ± 0.17 mm, and 1.52 ± 0.75 mm, respectively. The computation time of the proposed method for segmenting CBCT images of one subject was 7.25 ± 0.73 min. Compared with two other methods, the proposed method achieves significant improvement in terms of accuracy. Conclusions: The presented tooth segmentation method can be used to segment tooth contours from CT images accurately and efficiently.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Vibrational corrections to the second hyperpolarizabilities of Al{sub n}P{sub n} clusters

In this work, we report results of vibrational corrections to the second hyperpolarizabilities of Al{sub 2}P{sub 2}, Al{sub 3}P{sub 3}, Al{sub 4}P{sub 4}, Al{sub 6}P{sub 6}, and Al{sub 9}P{sub 9} clusters. The vibrational corrections were calculated through the perturbation theoretic method of Bishop and Kirtman and also using a variational methodology at the second order Møller-Plesset perturbation theory level with the aug-cc-pVDZ basis set. Results show that the vibrational corrections are important, accounting for more than half of the corresponding electronic second hyperpolarizabilities at the static limit. Comparisons between results obtained through both methods show very good agreements for the terms [α{sup 2}] and [μβ] but significant differences for the term [μ{sup 2}α]. Dynamic vibrational corrections to the second hyperpolarizabilities related to the dc-second harmonic generation, intensity dependent refractive index, and dc-Kerr nonlinear optical processes are also reported.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗