Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A Multi-Center Space Data System Prototype Based on CCSDS Standards

Deep space missions beyond earth orbit will require new methods of data communications in order to compensate for increasing RF propagation delay. The Consultative Committee for Space Data Systems (CCSDS) standard protocols Spacecraft Monitor & Control (SM&C), Asynchronous Message Service (AMS), and Delay/Disruption Tolerant Networking (DTN) provide such a method. The maturity level of this protocol set is, however, insufficient for mission inclusion at this time. This prototype is intended to provide experience which will raise the Technical Readiness Level (TRL) of these protocols..

Rich, Thomas M.↗
DOE OSTI2017

The Thick Level-Set model for dynamic fragmentation

The Thick Level-Set (TLS) model is implemented to simulate brittle media undergoing dynamic fragmentation. This non-local model is discretized by the finite element method with damage represented as a continuous field over the domain. A level-set function defines the extent and severity of damage, and a length scale is introduced to limit the damage gradient. Numerical studies in one dimension demonstrate that the proposed method reproduces the rate-dependent energy dissipation and fragment length observations from analytical, numerical, and experimental approaches. In conclusion, additional studies emphasize the importance of appropriate bulk constitutive models and sufficient spatial resolution of the length scale.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Approximate l-fold cross-validation with Least Squares SVM and Kernel Ridge Regression

Kernel methods have difficulties scaling to large modern data sets. The scalability issues are based on computational and memory requirements for working with a large matrix. These requirements have been addressed over the years by using low-rank kernel approximations or by improving the solvers scalability. However, Least Squares Support VectorMachines (LS-SVM), a popular SVM variant, and Kernel Ridge Regression still have several scalability issues. In particular, the O(n^3) computational complexity for solving a single model, and the overall computational complexity associated with tuning hyperparameters are still major problems. We address these problems by introducing an O(n log n) approximate l-fold cross-validation method that uses a multi-level circulant matrix to approximate the kernel. In addition, we prove our algorithm s computational complexity and present empirical runtimes on data sets with approximately 1 million data points. We also validate our approximate method s effectiveness at selecting hyperparameters on real world and standard benchmark data sets. Lastly, we provide experimental results on using a multi-level circulant kernel approximation to solve LS-SVM problems with hyperparameters selected using our method.

Edwards, Richard E↗
DOE OSTI2014

Approximate l-Fold Cross-Validation with Least Squares SVM and Kernel Ridge Regression

Kernel methods have difficulties scaling to large modern data sets. The scalability issues are based on computational and memory requirements for working with a large matrix. These requirements have been addressed over the years by using low-rank kernel approximations or by improving the solvers' scalability. However, Least Squares Support Vector Machines (LS-SVM), a popular SVM variant, and Kernel Ridge Regression still have several scalability issues. In particular, the O(n^3) computational complexity for solving a single model, and the overall computational complexity associated with tuning hyper parameters are still major problems. We address these problems by introducing an O(nlog n) approximate l-fold cross-validation method that uses a multi-level circulant matrix to approximate the kernel. In addition, we prove our algorithm's computational complexity and present empirical runtimes on data sets with approximately one million data points. We also validate our approximate method's effectiveness at selecting hyper parameters on real world and standard benchmark data sets. Lastly, we provide experimental results on using a multi level circulant kernel approximation to solve LS-SVM problems with hyper parameters selected using our method.

Computer Science↗
NASA NTRS2004

The Bonding of NO2, NH3, and CH2NY to Models of a (10,0) Carbon Nanotube

We have studied the bonding of NO2, NH3, and CH2NH to a (10,O) carbon nanotube using the MP2 and ONIOM methods with extended basis sets. We find bond strengths of 3.5, 3.6. and 6.3 kcal/mol for NO2, NH3, and CH2NH, respectively, using the ONIOM method with the high accuracy part treated at the MP2/aug-CC-pVTZ level and the remainder of the CNT approximated at the UFF level and including an estimate of basis set superposition error using the counterpoise method.

Walch, Stephen P.↗
DOE OSTI2023

High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries

Here, we propose a method that morphs high-order meshes such that their boundaries and interfaces coincide/align with implicitly defined geometries. Our focus is particularly on the case when the target surface is prescribed as the zero isocontour of a smooth discrete function. Common examples of this scenario include using level set functions to represent material interfaces in multimaterial configurations, and evolving geometries in shape and topology optimization. The proposed method formulates the mesh optimization problem as a variational minimization of the sum of a chosen mesh-quality metric using the Target-Matrix Optimization Paradigm (TMOP) and a penalty term that weakly forces the selected faces of the mesh to align with the target surface. The distinct features of the method are use of a source mesh to represent the level set function with sufficient accuracy, and adaptive strategies for setting the penalization weight and selecting the faces of the mesh to be fit to the target isocontour of the level set field. We demonstrate that the proposed method is robust for generating boundary- and interface-fitted meshes for curvilinear domains using different element types in 2D and 3D.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Total hydrocarbon analysis by ion mobility spectrometry

Astronauts must be alerted quickly to chemical leaks that compromise their health and the success of their missions. An ideal leak detector would be equally sensitive to all compounds that might constitute a hazard and insensitive to nontoxic compounds. No ideal sensor exists; thus, selection of a methodology is a series of compromises. The commonly used methods are either insensitive at the low exposure levels set by OSHA, NASA, and other organizations or are selectively insensitive to important classes of chemicals such as Freons. After extensive study and experience, the Toxicology Group at JSC has selected ion mobility spectrometry (IMS) for development into a broad range, sensitive detector. In addition to the sensing method, signal processing is important leak detection because a background signal can be expected at all times. The leak-detecting instrument must be programmed to discriminate between authentic leaks and background fluctuations caused by routine operations. The results of an evaluation of the prototype THA is presented in terms related to spacecraft operations. The evaluation included determination of instrumental parameters such as stability and response times. We also included responses to some common components of spacecraft atmospheres in pure form and in binary and ternary mixtures. The output of the four algorithms to the mixtures was found to be noticeably different. These responses are compared on the basis of their utility for signaling a chemical leak. As a means of evaluating its resistance to a falsely positive response, the THA was challenged with carbon dioxide and methane, compounds whose concentrations normally increase in spacecraft air during human habitation. The instrument showed virtually no response to these interferences. Although the prototype THA is designed for space flight, this detector is expected to be useful for field screening at chemical waste dumps and other environmentally sensitive locations.

Cross, John H.↗
DOE OSTI2015

A comparative study of diffraction of shallow-water waves by high-level IGN and GN equations

This work is on the nonlinear diffraction analysis of shallow-water waves, impinging on submerged obstacles, by two related theories, namely the classical Green–Naghdi (GN) equations and the Irrotational Green–Naghdi (IGN) equations, both sets of equations being at high levels and derived for incompressible and inviscid flows. Recently, the high-level Green–Naghdi equations have been applied to some wave transformation problems. The high-level IGN equations have also been used in the last decade to study certain wave propagation problems. However, past works on these theories used different numerical methods to solve these nonlinear and unsteady sets of differential equations and at different levels. Moreover, different physical problems have been solved in the past. Therefore, it has not been possible to understand the differences produced by these two sets of theories and their range of applicability so far. We are thus motivated to make a direct comparison of the results produced by these theories by use of the same numerical method to solve physically the same wave diffraction problems. We focus on comparing these two theories by using similar codes; only the equations used are different but other parts of the codes, such as the wave-maker, damping zone, discretion method, matrix solver, etc., are exactly the same. This way, we eliminate many potential sources of differences that could be produced by the solution of different equations. The physical problems include the presence of various submerged obstacles that can be used for example as breakwaters or to represent the continental shelf. A numerical wave tank is created by placing a wavemaker on one end and a wave absorbing beach on the other. The nonlinear and unsteady sets of differential equations are solved by the finite-difference method. The results are compared with different equations as well as with the available experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Distributed renewable energy resource trust metrics and secure routing

To ensure reliable and predictable service in the electrical grid between distributed renewable distributed energy resources (DERs) it is important to gauge the level of trust present within critical components and DER aggregators (DERAs). Although trust throughout a smart grid is temporal and dynamically varies according to measured states, it is possible to accurately formulate communications and service level strategies based on such trust measurements. Utilizing an effective set of machine learning and statistical methods, it is shown that establishment of trust levels between DERAs using behavioral pattern analysis is possible. Further, it is also shown that the establishment of such trust can facilitate simple secure communications routing between DERAs. Providing secure routing between DERAs enables a grid operator to maintain service level agreements to its customers, reduce the attack surface and increase operational resiliency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Dual baseline search for muon antineutrino disappearance at 0.1 eV²<Δm²<100 eV²

The MiniBooNE and SciBooNE collaborations report the results of a joint search for short baseline disappearance of ν¯ μ at Fermilab’s Booster Neutrino Beamline. The MiniBooNE Cherenkov detector and the SciBooNE tracking detector observe antineutrinos from the same beam, therefore the combined analysis of their data sets serves to partially constrain some of the flux and cross section uncertainties. Uncertainties in the ν μ background were constrained by neutrino flux and cross section measurements performed in both detectors. A likelihood ratio method was used to set a 90% confidence level upper limit on ν¯ μ disappearance that dramatically improves upon prior limits in the Δm²=0.1–100 eV² region.

Cheng, G.↗
DOE OSTI2025

A simple and highly efficient protocol for 13 C-labeling of plant cell wall for structural and quantitative analyses via solid-state nuclear magnetic resonance

Plant cell walls are made of a complex network of interacting polymers that play a critical role in plant development and responses to environmental changes. Thus, improving plant biomass and fitness requires the elucidation of the structural organization of plant cell walls in their native environment. The 13 C-based multi-dimensional solid-state nuclear magnetic resonance (ssNMR) has been instrumental in revealing the structural information of plant cell walls through 2D and 3D correlation spectral analyses. However, the requirement of enriching plants with 13 C limits the applicability of this method. To our knowledge, there is only a very limited set of methods currently available that achieve high levels of 13 C-labeling of plant materials using 13 CO 2 , and most of them require large amounts of 13 CO 2 in larger growth chambers. In this study, a simplified protocol for 13C-labeling of plant materials is introduced that allows ca 60% labeling of the cell walls, as quantified by comparison with commercially labeled samples. This level of 13 C-enrichment is sufficient for all conventional 2D and 3D correlation ssNMR experiments for detailed analysis of plant cell wall structure. The protocol is based on a convenient and easy setup to supply both 13 C-labeled glucose and 13 CO 2 using a vacuum-desiccator. The protocol does not require large amounts of 13 CO 2 . This study shows that our 13 C-labeling of plant materials can make the accessibility to ssNMR technique easy and affordable. The derived high-resolution 2D and 3D correlation spectra are used to extract structural information of plant cell walls. This helps to better understand the influence of polysaccharide-polysaccharide interaction on plant performance and allows for a more precise parametrization of plant cell wall models.

09 BIOMASS FUELS↗
DOE OSTI2013

Precision Search for Muon Antineutrino Disappearance Oscillations Using a Dual Baseline Technique

A search for short baseline muon antineutrino disappearance with the SciBooNE and MiniBooNE experiments at Fermi National Accelerator Laboratory in Batavia, Illinois is presented. Short baseline muon antineutrino disappearance measurements help constrain sterile neutrino models. The two detectors observe muon antineutrinos from the same beam, therefore the combined analysis of their data sets serves to partially constrain some of the flux and cross section uncertainties. A likelihood ratio method was used to set a 90% confidence level upper limit on muon antineutrino disappearance that dramatically improves upon prior sterile neutrino oscillation limits in the Δm 2 =0.1-100 eV 2 region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗
NASA NTRS2001

Left ventricular endocardial surface detection based on real-time 3D echocardiographic data

OBJECTIVE: A new computerized semi-automatic method for left ventricular (LV) chamber segmentation is presented. METHODS: The LV is imaged by real-time three-dimensional echocardiography (RT3DE). The surface detection model, based on level set techniques, is applied to RT3DE data for image analysis. The modified level set partial differential equation we use is solved by applying numerical methods for conservation laws. The initial conditions are manually established on some slices of the entire volume. The solution obtained for each slice is a contour line corresponding with the boundary between LV cavity and LV endocardium. RESULTS: The mathematical model has been applied to sequences of frames of human hearts (volume range: 34-109 ml) imaged by 2D and reconstructed off-line and RT3DE data. Volume estimation obtained by this new semi-automatic method shows an excellent correlation with those obtained by manual tracing (r = 0.992). Dynamic change of LV volume during the cardiac cycle is also obtained. CONCLUSION: The volume estimation method is accurate; edge based segmentation, image completion and volume reconstruction can be accomplished. The visualization technique also allows to navigate into the reconstructed volume and to display any section of the volume.

Non-NASA Center↗
DOE OSTI2018

Bayesian Static Parameter Estimation for Partially Observed Diffusions via Multilevel Monte Carlo

In this article we consider static Bayesian parameter estimation for partially observed diffusions that are discretely observed. We work under the assumption that one must resort to discretizing the underlying diffusion process, for instance, using the Euler--Maruyama method. Given this assumption, we show how one can use Markov chain Monte Carlo (MCMC) and particularly particle MCMC [C. Andrieu, A. Doucet, and R. Holenstein, J. R. Stat. Soc. Ser. B Stat. Methodol., 72 (2010), 269--342] to implement a new approximation of the multilevel (ML) Monte Carlo (MC) collapsing sum identity. Our approach comprises constructing an approximate coupling of the posterior density of the joint distribution over parameter and hidden variables at two different discretization levels and then correcting by an importance sampling method. The variance of the weights are independent of the length of the observed data set. The utility of such a method is that, for a prescribed level of mean square error, the cost of this MLMC method is provably less than i.i.d. sampling from the posterior associated to the most precise discretization. However the method here comprises using only known and efficient simulation methodologies. Finally, the theoretical results are illustrated by inference of the parameters of two prototypical processes given noisy partial observations of the process: the first is an Ornstein--Uhlenbeck process and the second is a more general Langevin equation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2018

Lessons and Recommendations for Board-Level Testing with Protons

Protons with sufficiently high energy, provided in a broad field covering on the order of 0.1m2 can be used to perform board-level testing for single event effects (SEE). NASA has used this approach for board-level testing over the last 20 years. Although many difficulties inherent in SEE testing are simplified when using a board-level test, including reduced cost, the method is inherently risky because of the limited value of the collected data and the potential to make critical mistakes when performing SEE testing this way, leading to data of less value. Historically, NASA’s approach to proton board-level testing has been limited to lower criticality applications. However, with users both inside and outside NASA using this method for higher levels of mission assurance, we have put together a set of lessons and recommendations to improve the value of data collected using this method. Focus areas covered include test preparation, test execution, and interpretation of results.

Guertin, Steven M.↗
DOE OSTI2018

Relativistic MR-MP Energy Levels for L-shell Ions of Sulfur and Argon

Calculated level energies for valence and K-vacancy states are given for the ion series S vii–S xiv and Ar ix–Ar xvi. The calculations were performed with the relativistic Multi-Reference Møller–Plesset Perturbation Theory method (MR-MP). The data set includes all the level energies with configurations $1{s}^{2}2{(s,p)}^{q}$, $1{s}^{2}2{(s,p)}^{q-1}{nl}$, $1{s}^{1}2{(s,p)}^{q+1}$, $1{s}^{1}2{(s,p)}^{q}{nl}$, $2{(s,p)}^{q+2}$, and $2{(s,p)}^{q+1}{nl}$, where 1 ≤ q ≤ 8, n ≤ 5, and l ≤ 3. We have compared our results with data from the National Institute of Standards and Technology (NIST) online database and with previous calculations. The average deviation of valence level energies ranges from 0.16 eV in Ne-like ions to 0.01 eV in Li-like ions, demonstrating that the present MR-MP valence level energies are highly accurate. In the case of K-vacancy states, the deviation is generally below 0.3 eV for Li-like S xiv and Ar xvi. The deviation for K-vacancy energies in other L-shell ions (Be-, B-, C-, N-, and O-like Ar ions) is higher but likely because the NIST-recommended values have a higher uncertainty. The data set includes many n = 4 and n = 5 valence and K-vacancy levels in L-shell ions of S and Ar that have not been previously reported. The data can be used for line identification and modeling of L-shell ions of S and Ar in astrophysical and laboratory-generated plasmas, and as energy references in the absence of more accurate laboratory calculations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Anharmonic Bloch oscillations of electrons in electrically biased superlattices

The phenomenon of anharmonic Bloch oscillations (i.e., oscillations with a frequency which is a multiple of the Bloch frequency) is considered. The energy-band structure of silicon-carbide polytypes where these oscillations are observed is calculated ab initio. A one-dimensional model potential making it possible to calculate the Stark wave functions in a biased superlattice is constructed for these polytypes. The transfer matrix method is used to find the set of energy levels (the so-called Stark ladder) and calculate the electron wave functions. It is shown that both radiative transitions between neighboring Stark levels (at the Bloch frequency) and transitions in the form of jumps via several levels, accompanied by emission at frequencies that are multiples of the Bloch frequency, may occur. The calculated probabilities of transitions between Stark-ladder levels increase with increasing applied field strength.

36 MATERIALS SCIENCE↗