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 487 records · Page 27

The Failure of Correlation to Describe Out-of-Plane Carbon=Carbon Bending

Carbon-carbon multiply bonded systems are improperly described with standard correlation methods and basis sets. For computations of vibrational modes, the out-of- plane bends can be reported as imaginary at worst or simply too low at best. Utilizing the simplest of aromatic structures (cyclopropenylidene) and various levels of theory, this work diagnoses this known behavior for the first time. A combined 1-particle and n-particle basis set effect conspire to produce these non-physical results. When moving from sp2 to sp3 hybridization in the carbon atoms, the larger number of basis functions overcorrects the energy. This is exacerbated by correlation methods. These allow for occupation of the π and π∗ orbitals in the expanded wave function that combine with the hydrogen s orbitals. As a result, the improperly described space can be further and non-physically stabilized by post-Hartree-Fock correlation. This represents a fundamental problem with at least Hartree-Fock based methods of all flavors in describing carbon. Beyond being a flaw in quantum chemical theory, other repercussions will be present in computations regarding spectroscopy as well as energy and environmental studies where highly-accurate hydrocabon vibrational transitions or thermochemical data are needed.

Fortenberry, Ryan C.↗

A Phase I/II Trial of Intensity Modulated Radiation (IMRT) Dose Escalation With Concurrent Fixed-dose Rate Gemcitabine (FDR-G) in Patients With Unresectable Pancreatic Cancer

Purpose: Local failure in unresectable pancreatic cancer may contribute to death. We hypothesized that intensification of local therapy would improve local control and survival. The objectives were to determine the maximum tolerated radiation dose delivered by intensity modulated radiation with fixed-dose rate gemcitabine (FDR-G), freedom from local progression (FFLP), and overall survival (OS). Methods and Materials: Eligibility included pathologic confirmation of adenocarcinoma, radiographically unresectable, performance status of 0-2, absolute neutrophil count of {>=}1500/mm{sup 3}, platelets {>=}100,000/mm{sup 3}, creatinine <2 mg/dL, bilirubin <3 mg/dL, and alanine aminotransferase/aspartate aminotransferase {<=}2.5 Multiplication-Sign upper limit of normal. FDR-G (1000 mg/m{sup 2}/100 min intravenously) was given on days -22 and -15, 1, 8, 22, and 29. Intensity modulated radiation started on day 1. Dose levels were escalated from 50-60 Gy in 25 fractions. Dose-limiting toxicity was defined as gastrointestinal toxicity grade (G) {>=}3, neutropenic fever, or deterioration in performance status to {>=}3 between day 1 and 126. Dose level was assigned using TITE-CRM (Time-to-Event Continual Reassessment Method) with the target dose-limiting toxicity (DLT) rate set to 0.25. Results: Fifty patients were accrued. DLTs were observed in 11 patients: G3/4 anorexia, nausea, vomiting, and/or dehydration (7); duodenal bleed (3); duodenal perforation (1). The recommended dose is 55 Gy, producing a probability of DLT of 0.24. The 2-year FFLP is 59% (95% confidence interval [CI]: 32-79). Median and 2-year overall survival are 14.8 months (95% CI: 12.6-22.2) and 30% (95% CI 17-45). Twelve patients underwent resection (10 R0, 2 R1) and survived a median of 32 months. Conclusions: High-dose radiation therapy with concurrent FDR-G can be delivered safely. The encouraging efficacy data suggest that outcome may be improved in unresectable patients through intensification of local therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Tapping the treasure trove of atypical phages

With advancements in genomics technologies, a vast diversity of ‘atypical’ phages, that is, with single-stranded DNA or RNA genomes, are being uncovered from different ecosystems. Though these efforts have revealed the existence and prevalence of these nonmodel phages, computational approaches often fail to associate these phages with their specific bacterial host(s), while the lack of methods to isolate these phages has limited our ability to characterize infectivity pathways and new gene function. In this review, we call for the development of generalizable experimental methods to better capture this understudied viral diversity via isolation and study them through gene-level characterization and engineering. Establishing a diverse set of new ‘atypical’ phage model systems has the potential to provide many new biotechnologies, including potential uses of these atypical phages in halting the spread of antibiotic resistance and engineering of microbial communities for beneficial outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

SHαDE: survey description and mass–kinematics scaling relations for dwarf galaxies

ABSTRACT The Study of H α from Dwarf Emissions (SHαDE) is a high spectral resolution (R = 13 500) H α integral field survey of 69 dwarf galaxies with stellar masses 106 < M⋆ < 109 M⊙. The survey used FLAMES on the ESO Very Large Telescope. SHαDE is designed to study the kinematics and stellar populations of dwarf galaxies using consistent methods applied to massive galaxies and at matching level of detail, connecting these mass ranges in an unbiased way. In this paper, we set out the science goals of SHαDE, describe the sample properties, outline the data reduction, and analysis processes. We investigate the log M⋆–log S0.5 mass–kinematics scaling relation, which has previously shown potential for combining galaxies of all morphologies in a single scaling relation. We extend the scaling relation from massive galaxies to dwarf galaxies, demonstrating this relation is linear down to a stellar mass of M⋆ ∼ 108.6 M⊙. Below this limit, the kinematics of galaxies inside one effective radius appears to be dominated by the internal velocity dispersion limit of the H α-emitting gas, giving a bend in the log M⋆–log S0.5 relation. Replacing stellar mass with total baryonic mass using gas mass estimate reduces the severity but does not remove the linearity limit of the scaling relation. An extrapolation to estimate the galaxies’ dark matter halo masses, yields a log Mh–log S0.5 scaling relation that is free of any bend, has reduced curvature over the whole mass range, and brings galaxies of all masses and morphologies on to the virial relation.

Barat, Dilyar↗

Level-set Modeling Simulations of Chemical Vapor Infiltration for Ceramic Matrix Composites Manufacturing

Silicon-carbide (SiC) reinforced ceramic matrix composites (CMCs) are a key enabling technology to reduce fuel consumption and emissions of gas turbine engines. In one manufacturing approach, chemical vapor infiltration (CVI) is limited to only coating SiC fibers. The preform is then fabricated using a lay-up of basic plys or 2-D woven sheets composed of the precoated fibers. At the other extreme, CVI is used to completely densify a 3-D woven preform shaped almost like the gas turbine component itself. The latter approach is more suitable for highly engineered components which sit directly in the gas path of the engine, for example, a high pressure turbine blade. In this case, the geometry is necessarily complex for aerodynamic, stress, and lifing (multi-physics) requirements. Presently, optimizing the CVI-dominated manufacturing approach is largely by trial-and-error. In this work, a first-principles modeling of CVI is performed to realize optimization of SiC/SiC CMC manufacturing. The modeling is based on a level-set framework to describe the interface between the vapor and solid phases. A finite-difference numerical scheme using an immersed boundary method is developed for fixed, structured meshes. Massively parallel direct numerical simulations (DNS) of CVI through fiber-woven geometries are performed using one-step chemistry, and over a range of Thiele moduli. Illustrative applications of the resulting large DNS data sets are given, including the development of fiber-weave specific infiltration models and structure functions for mean-field (porous media) Computational Fluid Dynamics (CFD) simulations of CVI.

36 MATERIALS SCIENCE↗

Mariner Mars 9 stereophotogrammetry

Discussion of the equipment and orbital specifications of the Mariner Mars 9 spacecraft - a Martian orbiter which was placed in orbit on Nov. 14, 1971, for an expected 50-year stay and has made possible a practically 100% photographic coverage of the planet's surface. The orbit has a nominal 12-hr period, a nominal periapsis of 1300 km, and a nominal apoapsis of 18,000 km. Many of the outstanding topographic features of the planet have been covered by the convergent stereopairs carried by the orbiter. Theoretical considerations are set forth concerning the necessary instrument orientation and rotation requirements to achieve such performance levels in this extraterrestrial application of photogrammetry. The photogrammetric method used in this mission is assessed as useful in the evaluation of outstanding Martian features such as the mammoth volcano of the Nix Olympica region.

Benesh, M.↗

On instruction sets and their formats

Central to instruction set design is the issue of the instruction format. Some common format encoding techniques are discussed, and a method of representing high-level language parse trees by means of formats that cover successively larger portions of a tree is introduced. Variations are introduced on the method that represent directed acyclic graphs as well as simple parse trees, and that encode constants in a special fashion. For a particular representation, the number of times each format is executed to run a sample program to completion is measured.

Flynn, M. J.↗

The vibrational frequencies of ozone

The singles and doubles coupled-cluster (CCSD) method is used, together with the recently proposed CCSD(T) and CCSDT-1 methods, to determine the equilibrium structure and harmonic frequencies of ozone at various levels, on the basis of large one-particle basis sets. It is demonstrated that f-type functions are required in order to obtain quantitative accuracy for the equilibrium structure and harmonic frequencies. The CCSD(T) level of theory, in conjunction with a (4s 3p 2d 1f) atomic natural orbital set, yields an equilibrium geometry and harmonic frequencies that are in excellent agreement with experimental results; an effort has accordingly been made to determine the cubic and quartic force constants, and evaluate the anharmonic corrections to the fundamental vibrational frequencies.

Lee, Timothy J.↗

NREL Controllable Grid Interface for Testing MW-Scale Wind Turbine Generators (Poster)

In order to understand the behavior of wind turbines experiencing grid disturbances, it is necessary to perform a series of tests and accurate transient simulation studies. The latest edition of the IEC 61400-21 standard describes methods for such tests that include low voltage ride-through (LVRT), active power set-point control, ramp rate limitations, and reactive power capability tests. The IEC methods are being widely adopted on both national and international levels by wind turbine manufacturers, certification authorities, and utilities. On-site testing of wind turbines might be expensive and time consuming since it requires both test equipment transportation and personnel presence in sometimes remote locations for significant periods of time because such tests need to be conducted at certain wind speed and grid conditions. Changes in turbine control software or design modifications may require redoing of all tests. Significant cost and test-time reduction can be achieved if these tests are conducted in controlled laboratory environments that replicate grid disturbances and simulation of wind turbine interactions with power systems. Such testing capability does not exist in the United States today. An initiative by NREL to design and construct a 7-MVA grid simulator to operate with the existing 2.5 MW and new upcoming 5-MW dynamometer facilities will fulfill this role and bring many potential benefits to the U.S. wind industry with the ultimate goal of reducing wind energy integration costs.

17 WIND ENERGY↗

Propagating sample variance uncertainties in redshift calibration: simulations, theory, and application to the COSMOS2015 data

ABSTRACT Cosmological analyses of galaxy surveys rely on knowledge of the redshift distribution of their galaxy sample. This is usually derived from a spectroscopic and/or many-band photometric calibrator survey of a small patch of sky. The uncertainties in the redshift distribution of the calibrator sample include a contribution from shot noise, or Poisson sampling errors, but, given the small volume they probe, they are dominated by sample variance introduced by large-scale structures. Redshift uncertainties have been shown to constitute one of the leading contributions to systematic uncertainties in cosmological inferences from weak lensing and galaxy clustering, and hence they must be propagated through the analyses. In this work, we study the effects of sample variance on small-area redshift surveys, from theory to simulations to the COSMOS2015 data set. We present a three-step Dirichlet method of resampling a given survey-based redshift calibration distribution to enable the propagation of both shot noise and sample variance uncertainties. The method can accommodate different levels of prior confidence on different redshift sources. This method can be applied to any calibration sample with known redshifts and phenotypes (i.e. cells in a self-organizing map, or some other way of discretizing photometric space), and provides a simple way of propagating prior redshift uncertainties into cosmological analyses. As a worked example, we apply the full scheme to the COSMOS2015 data set, for which we also present a new, principled SOM algorithm designed to handle noisy photometric data. We make available a catalogue of the resulting resamplings of the COSMOS2015 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmarking the extraction of statistical neutron capture cross sections on short-lived nuclei for applications using the β -Oslo method

Numerous scientific fields including astrophysics, nuclear power, and nuclear forensics require a knowledge of basic nuclear properties for large numbers of short-lived, radioactive isotopes far removed from stable nuclei. Neutron-capture cross sections are one such piece of nuclear data where direct measurements are not possible and theoretical predictions can vary by orders of magnitude. A recently developed indirect technique for inferring neutron capture rates, the β-Oslo method, has been introduced but not compared against a known, directly measured neutron capture cross section. To provide this benchmark, two indirect methods based on β decay and charged-particle reactions were used to extract the nuclear level density and γ-ray strength function of 51 Ti. The nuclear level density and γ-ray strength function from the two data sets were found to be equivalent and were used to extract the neutron capture cross section of 50 Ti which agrees with previous direct measurements at high neutron energies. Here, the results demonstrate the validity of the β-Oslo method for extracting neutron capture cross sections of short-lived nuclei and provide a sufficiently small uncertainty to be used in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theoretical studies for the N{sub 2}–N{sub 2}O van der Waals complex: The potential energy surface, intermolecular vibrations, and rotational transition frequencies

Theoretical studies of the potential energy surface (PES) and bound states are performed for the N{sub 2}–N{sub 2}O van der Waals (vdW) complex. A four-dimensional intermolecular PES is constructed at the level of single and double excitation coupled-cluster method with a non-iterative perturbation treatment of triple excitations [CCSD(T)] with aug-cc-pVTZ basis set supplemented with bond functions. Two equivalent T-shaped global minima are located, in which the O atom of N{sub 2}O monomer is near the N{sub 2} monomer. The intermolecular fundamental vibrational states are assigned by inspecting the orientation of the nodal surface of the wavefunctions. The calculated frequency for intermolecular disrotation mode is 23.086 cm{sup −1}, which is in good agreement with the available experimental data of 22.334 cm{sup −1}. A negligible tunneling splitting with the value of 4.2 MHz is determined for the ground vibrational state and the tunneling splitting increases as the increment of the vibrational frequencies. Rotational levels and transition frequencies are calculated for both isotopomers {sup 14}N{sub 2}–N{sub 2}O and {sup 15}N{sub 2}–N{sub 2}O. The accuracy of the PES is validated by the good agreement between theoretical and experimental results for the transition frequencies and spectroscopic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Synthesis and Control of Flexible Systems with Component-Level Uncertainties

An efficient and computationally robust method for synthesis of component dynamics is developed. The method defines the interface forces/moments as feasible vectors in transformed coordinates to ensure that connectivity requirements of the combined structure are met. The synthesized system is then defined in a transformed set of feasible coordinates. The simplicity of form is exploited to effectively deal with modeling parametric and non-parametric uncertainties at the substructure level. Uncertainty models of reasonable size and complexity are synthesized for the combined structure from those in the substructure models. In particular, we address frequency and damping uncertainties at the component level. The approach first considers the robustness of synthesized flexible systems. It is then extended to deal with non-synthesized dynamic models with component-level uncertainties by projecting uncertainties to the system level. A numerical example is given to demonstrate the feasibility of the proposed approach.

Maghami, Peiman G.↗

Comparison of different methods for estimating snowcover in forested, mountainous basins using LANDSAT (ERTS) images

Snow-covered areas on LANDSAT (ERTS) images of the Santiam River basin, Oregon, and other basins in Washington were measured using several operators and methods. Seven methods were used: (1) Snowline tracing followed by measurement with planimeter, (2) mean snowline altitudes determined from many locations, (3) estimates in 2.5 x 2.5 km boxes of snow-covered area with reference to snow-free images, (4) single radiance-threshold level for entire basin, (5) radiance-threshold setting locally edited by reference to altitude contours and other images, (6) two-band color-sensitive extraction locally edited as in (5), and (7) digital (spectral) pattern recognition techniques. The seven methods are compared in regard to speed of measurement, precision, the ability to recognize snow in deep shadow or in trees, relative cost, and whether useful supplemental data are produced.

Meier, M. J.↗

Review and evaluation of recent developments in melic inlet dynamic flow distortion prediction and computer program documentation and user's manual estimating maximum instantaneous inlet flow distortion from steady-state total pressure measurements with full, limited, or no dynamic data

A review of the Melick method of inlet flow dynamic distortion prediction by statistical means is provided. These developments include the general Melick approach with full dynamic measurements, a limited dynamic measurement approach, and a turbulence modelling approach which requires no dynamic rms pressure fluctuation measurements. These modifications are evaluated by comparing predicted and measured peak instantaneous distortion levels from provisional inlet data sets. A nonlinear mean-line following vortex model is proposed and evaluated as a potential criterion for improving the peak instantaneous distortion map generated from the conventional linear vortex of the Melick method. The model is simplified to a series of linear vortex segments which lay along the mean line. Maps generated with this new approach are compared with conventionally generated maps, as well as measured peak instantaneous maps. Inlet data sets include subsonic, transonic, and supersonic inlets under various flight conditions.

Schweikhard, W. G.↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-398: Evaluation of Radiobiological Parameters Using Serial Tumor Imaging During Radiotherapy as An Inverse Ill-Posed Problem

Purpose: Combination of serial tumor imaging with radiobiological modeling can provide more accurate information on the nature of treatment response and what underlies resistance. The purpose of this article is to improve the algorithms related to imaging-based radiobilogical modeling of tumor response. Methods: Serial imaging of tumor response to radiation therapy represents a sum of tumor cell sensitivity, tumor growth rates, and the rate of cell loss which are not separated explicitly. Accurate treatment response assessment would require separation of these radiobiological determinants of treatment response because they define tumor control probability. We show that the problem of reconstruction of radiobiological parameters from serial imaging data can be considered as inverse ill-posed problem described by the Fredholm integral equation of the first kind because it is governed by a sum of several exponential processes. Therefore, the parameter reconstruction can be solved using regularization methods. Results: To study the reconstruction problem, we used a set of serial CT imaging data for the head and neck cancer and a two-level cell population model of tumor response which separates the entire tumor cell population in two subpopulations of viable and lethally damage cells. The reconstruction was done using a least squared objective function and a simulated annealing algorithm. Using in vitro data for radiobiological parameters as reference data, we shown that the reconstructed values of cell surviving fractions and potential doubling time exhibit non-physical fluctuations if no stabilization algorithms are applied. The variational regularization allowed us to obtain statistical distribution for cell surviving fractions and cell number doubling times comparable to in vitro data. Conclusion: Our results indicate that using variational regularization can increase the number of free parameters in the model and open the way to development of more advanced algorithms which take into account tumor heterogeneity, for example, related to hypoxia.

60 APPLIED LIFE SCIENCES↗

A method for calculating transient thrust and flow-rate levels for Mariner type attitude control nitrogen gas jets

The purpose of this report is to define and program the transient pneumatic flow equations necessary to determine, for a given set of conditions (geometry, pressures, temperatures, valve on time, etc.), the total nitrogen impulse and mass flow per pulse for the single pulsing of a Mariner type reaction control assembly valve. The rates of opening and closing of the valves are modeled, and electrical pulse durations from 20 to 100 ms are investigated. In developing the transient flow analysis, maximum use was made of the steady-state analysis. The impulse results are also compared to an equivalent square-wave impulse for both the Mariner Mars 1971 (MM'71) and Mariner Mars 1964 (MM'64) systems. It is demonstrated that, whereas in the MM'64 system, the actual impulse was as much as 56 percent higher than an assumed impulse (which is the product of the steady-state thrust and value on time i.e., the square wave), in the MM'71 system, these two values were in error in the same direction by only approximately 4 percent because of the larger nozzle areas and shorter valve stroke used.

Ferrera, J. D.↗