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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Scale separation for multi-scale modeling of free-surface and two-phase flows with the conservative sharp interface method

In this paper we present a scale separation approach for multi-scale modeling of free-surface and two-phase flows with complex interface evolution. By performing a stimulus-response operation on the level-set function representing the interface, separation of resolvable and non-resolvable interface scales is achieved efficiently. Uniform positive and negative shifts of the level-set function are used to determine non-resolvable interface structures. Non-resolved interface structures are separated from the resolved ones and can be treated by a mixing model or a Lagrangian-particle model in order to preserve mass. Resolved interface structures are treated by the conservative sharp-interface model. Since the proposed scale separation approach does not rely on topological information, unlike in previous work, it can be implemented in a straightforward fashion into a given level set based interface model. A number of two- and three-dimensional numerical tests demonstrate that the proposed method is able to cope with complex interface variations accurately and significantly increases robustness against underresolved interface structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1965

Lower Bounds for the Energy Levels of Anharmonic Oscillators

The lower‐bounds method presented by Löwdin in the preceding paper has been applied to oscillators perturbed by third‐ and fourth‐power terms in the potential‐energy expression. For favorable cases agreement between upper and lower bounds is easily carried to many more figures than are likely to be physically significant. In many cases the lower bounds agreed more closely to the true eigenvalue than did the corresponding upper bounds. For a given basis set, this method gives closer bounds than that of Bazley and Fox, except for energy levels too high to be satisfactorily treated in the given basis. The only disadvantage found was that for close bounds double precision proved necessary, indicating more than ordinary loss of computational accuracy.

Reid, Charles E.↗
DOE OSTI2020

Learning-Based Load Control to Support Resilient Networked Microgrid Operations

Microgrids have proven to be an effective option for increasing the resiliency of critical end-use loads during extreme events. Building on past operational experiences, some microgrid operators are examining the potential to network microgrids to further improve resiliency. However, the frequency deviations experienced on isolated microgrids during transient events, such as switching operations, step increases in load, and loss of generation, are significantly larger than those typically seen on bulk transmission systems. The larger frequency deviations can cause a loss of inverter-connected assets, resulting in a loss of power to critical end-use loads. This paper presents a method of mitigating the impact of transient events by engaging end-use loads using Grid-Friendly Appliance TM (GFA) controllers. An online, i.e., real-time, device-level algorithm is presented, which adjusts individual GFA controller frequency set-points based on the operational characteristics of each end-use load, and on the changing grid dynamic characteristics. The presented method improves the dynamic stability of the networked microgrid operations while minimizing the interruptions to end-use loads. The presented work is validated with dynamic simulations using a modified version of the IEEE 123-node test system with three microgrids, using the GridLAB-D TM simulation environment.

Radhakrishnan, Nikitha↗
DOE OSTI2016

MO-AB-BRA-01: A Global Level Set Based Formulation for Volumetric Modulated Arc Therapy

Purpose: The current clinical Volumetric Modulated Arc Therapy (VMAT) optimization is formulated as a non-convex problem and various greedy heuristics have been employed for an empirical solution, jeopardizing plan consistency and quality. We introduce a novel global direct aperture optimization method for VMAT to overcome these limitations. Methods: The global VMAT (gVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term and an anisotropic total variation term. A level set function was used to describe the aperture shapes and adjacent aperture shapes were penalized to control MLC motion range. An alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc gVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme (GBM), lung (LNG), and 2 head and neck cases—one with 3 PTVs (H&N3PTV) and one with 4 PTVs (H&N4PTV). The plans were compared against the clinical VMAT (cVMAT) plans utilizing two overlapping coplanar arcs. Results: The optimization of the gVMAT plans had converged within 600 iterations. gVMAT reduced the average max and mean OAR dose by 6.59% and 7.45% of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N3PTV case. PTV coverages (D95, D98, D99) were within 0.25% of the prescription dose. By globally considering all beams, the gVMAT optimizer allowed some beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization. Conclusions: The novel VMAT approach allows for the search of an optimal plan in the global solution space and generates deliverable apertures directly. The single arc VMAT approach fully utilizes the digital linacs’ capability in dose rate and gantry rotation speed modulation. Varian Medical Systems, NIH grant R01CA188300, NIH grant R43CA183390.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1993

Predesign of diamond turned refractive/diffractive elements for IR objectives

Diffractive lenses have arrived. Literally hundreds of papers have been published and technology impact reports have been written about the exciting addition of a new tool for the lens designer. Sophisticated computer programs have been developed to aid in the optimization of these diffractive phase profiles for a wide variety of applications. Now, several fabrication methods are being pursued to produce these diffractive elements economically. The best known process is the etching of a multi-level relief grating, known as binary optics. This process uses sets of computer generated lithographic masks. Another, more recently developed method is Dry Photopolymer Embossing (DPE). This replication process uses master holograms. And now, diamond turning is being applied for the machining of these elements. Diamond turning is especially well suited for infrared optics. As any process has advantages and limitations, so has diamond turning. These advantages and limitations are discussed and general guidelines are presented to aid the designer and systems engineer in the project predesign stage.

Riedl, Max J.↗
DOE OSTI2018

A new analysis approach based on Haralick texture features for the characterization of microstructure on the example of low-alloy steels

Highlights: • Microstructures were characterized based on orientation-invariant image-texture features. • A new parameter called “amplitude”, was formulated and is based on image rotation and the Haralick parameters: texture contrast, correlation, homogeneity and energy. • The proposed method is insensitive to changing image resolution and image contrast/brightness settings. - Abstract: Microstructures were analyzed by an improved texture-based method using gray level co-occurrence matrices (GLCM). This method is based on a new parameter calculated from the stepwise rotation of images and thereby, calculating the values independent of the original texture orientation. The proposed method was applied on a database of etched and scanning electron microscopy (SEM)-imaged low-carbon steel microstructures that are currently extensively used for automated microstructure classification. The results on the microstructures consisting of pearlitic, lath martensitic and lower bainitic constituents revealed that the method allows a significant separation of various types of microstructures in the ideal case of square-shaped cutouts. For complete grains of the corresponding second phases, the results imply that the application of a classifier is advantageous to distinguish them with a sufficient accuracy. The robustness and workability of the method was further demonstrated by discussing the effect of varying the image resolution and contrast/brightness settings during image acquisition. It was shown that such user-dependent setting parameters do not impair the separability of the steel constituents by using the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2019

An empirical assessment of a single family‐wide hybrid capture locus set at multiple evolutionary timescales in Asteraceae

Premise Hybrid capture with high‐throughput sequencing (Hyb‐Seq) is a powerful tool for evolutionary studies. The applicability of an Asteraceae family‐specific Hyb‐Seq probe set and the outcomes of different phylogenetic analyses are investigated here. Methods Hyb‐Seq data from 112 Asteraceae samples were organized into groups at different taxonomic levels (tribe, genus, and species). For each group, data sets of non‐paralogous loci were built and proportions of parsimony informative characters estimated. The impacts of analyzing alternative data sets, removing long branches, and type of analysis on tree resolution and inferred topologies were investigated in tribe Cichorieae. Results Alignments of the Asteraceae family‐wide Hyb‐Seq locus set were parsimony informative at all taxonomic levels. Levels of resolution and topologies inferred at shallower nodes differed depending on the locus data set and the type of analysis, and were affected by the presence of long branches. Discussion The approach used to build a Hyb‐Seq locus data set influenced resolution and topologies inferred in phylogenetic analyses. Removal of long branches improved the reliability of topological inferences in maximum likelihood analyses. The Astereaceae Hyb‐Seq probe set is applicable at multiple taxonomic depths, which demonstrates that probe sets do not necessarily need to be lineage‐specific.

Jones, Katy E.↗
DOE OSTI2024

Exploring the Tl 2H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Abstract Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet TlH potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 MH (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Chemistry↗
DOE OSTI2016

High-Resolution Infrared Studies of Perdeutero-Spiropentane, C 5 D 8

Perdeutero-spiropentane (C 5 D 8 ) has been synthesized and infrared and Raman spectra are reported for the first time. Wavenumber assignments are made for most of the fundamental vibrational modes. Gas phase infrared spectra were recorded at a resolution (0.002 cm -1 ) sufficient to resolve individual rovibrational lines and show evidence of strong Coriolis and/or Fermi resonance interactions for most bands. However a detailed rovibrational analysis of the fundamental v 15 (b 2 ) parallel band proved possible and a fit of more than 1600 lines yielded a band origin of 1053.84465(10) cm -1 and ground state constants (in units of cm -1 ): B 0 = 0.1120700(9), D J = 1.51(3) x10 -8 , D JK = 3.42(15) x10 -8 . We note that the B 0 value is significantly less than a value of Ba = 0.1140 cm-1 calculated using structural parameters from an earlier electron diffraction (ED) study, whereas one expects B a to be lower than B 0 because of thermal averaging over higher vibrational levels. A similar discrepancy was noted in an earlier study of C 5 H 8 [1]. The structural and spectroscopic results are in good accord with values computed at the anharmonic level using the B3LYP density functional method with a cc-pVTZ basis set.

C5H8 and C5D8↗
DOE OSTI2014

Structure and dynamics of small van der Waals complexes

We illustrate computational aspects of the calculation of the potential energy surfaces of small (up to five atoms) van der Waals complexes with high-level quantum chemistry techniques such as the CCSD(T) method with extended basis sets. We discuss the compromise between the required accuracy and the computational time. Further, we show how these potential energy surfaces can be fitted and used in dynamical calculations such as non-reactive inelastic scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

An interface tracking model for droplet electrocoalescence.

This report describes an Early Career Laboratory Directed Research and Development (LDRD) project to develop an interface tracking model for droplet electrocoalescence. Many fluid-based technologies rely on electrical fields to control the motion of droplets, e.g. microfluidic devices for high-speed droplet sorting, solution separation for chemical detectors, and purification of biodiesel fuel. Precise control over droplets is crucial to these applications. However, electric fields can induce complex and unpredictable fluid dynamics. Recent experiments (Ristenpart et al. 2009) have demonstrated that oppositely charged droplets bounce rather than coalesce in the presence of strong electric fields. A transient aqueous bridge forms between approaching drops prior to pinch-off. This observation applies to many types of fluids, but neither theory nor experiments have been able to offer a satisfactory explanation. Analytic hydrodynamic approximations for interfaces become invalid near coalescence, and therefore detailed numerical simulations are necessary. This is a computationally challenging problem that involves tracking a moving interface and solving complex multi-physics and multi-scale dynamics, which are beyond the capabilities of most state-of-the-art simulations. An interface-tracking model for electro-coalescence can provide a new perspective to a variety of applications in which interfacial physics are coupled with electrodynamics, including electro-osmosis, fabrication of microelectronics, fuel atomization, oil dehydration, nuclear waste reprocessing and solution separation for chemical detectors. We present a conformal decomposition finite element (CDFEM) interface-tracking method for the electrohydrodynamics of two-phase flow to demonstrate electro-coalescence. CDFEM is a sharp interface method that decomposes elements along fluid-fluid boundaries and uses a level set function to represent the interface.

Erickson, Lindsay Crowl↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
NASA NTRS1994

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(CO)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule O4(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI,and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in O4(+) has also been studied using the Brueckner coupled-cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.↗
DOE OSTI2015

SU-E-J-131: Augmenting Atlas-Based Segmentation by Incorporating Image Features Proximal to the Atlas Contours

Purpose: For facilitating the current automatic segmentation, in this work we propose a narrow-shell strategy to enhance the information of each contour in the library and to improve the accuracy of the exiting atlas-based approach. Methods: In setting up an atlas-based library, we include not only the coordinates of contour points, but also the image features adjacent to the contour. 139 planning CT scans with normal appearing livers obtained during their radiotherapy treatment planning were used to construct the library. The CT images within the library were registered each other using affine registration. A nonlinear narrow shell with the regional thickness determined by the distance between two vertices alongside the contour. The narrow shell was automatically constructed both inside and outside of the liver contours. The common image features within narrow shell between a new case and a library case were first selected by a Speed-up Robust Features (SURF) strategy. A deformable registration was then performed using a thin plate splines (TPS) technique. The contour associated with the library case was propagated automatically onto the images of the new patient by exploiting the deformation field vectors. The liver contour was finally obtained by employing level set based energy function within the narrow shell. The performance of the proposed method was evaluated by comparing quantitatively the auto-segmentation results with that delineated by a physician. Results: Application of the technique to 30 liver cases suggested that the technique was capable of reliably segment organs such as the liver with little human intervention. Compared with the manual segmentation results by a physician, the average and discrepancies of the volumetric overlap percentage (VOP) was found to be 92.43%+2.14%. Conclusion: Incorporation of image features into the library contours improves the currently available atlas-based auto-contouring techniques and provides a clinically practical solution for auto-segmentation. This work is supported by NIH/NIBIB (1R01-EB016777), National Natural Science Foundation of China (No.61471226 and No.61201441), Research funding from Shandong Province (No.BS2012DX038 and No.J12LN23), and Research funding from Jinan City (No.201401221 and No.20120109)

60 APPLIED LIFE SCIENCES↗
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↗