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 1,513 records · Page 84
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Waste Determination Equivalency - 12172

The Savannah River Site (SRS) is a Department of Energy (DOE) facility encompassing approximately 800 square kilometers near Aiken, South Carolina which began operations in the 1950's with the mission to produce nuclear materials. The SRS contains fifty-one tanks (2 stabilized, 49 yet to be closed) distributed between two liquid radioactive waste storage facilities at SRS containing carbon steel underground tanks with storage capacities ranging from 2,800,000 to 4,900,000 liters. Treatment of the liquid waste from these tanks is essential both to closing older tanks and to maintaining space needed to treat the waste that is eventually vitrified or disposed of onsite. Section 3116 of the Ronald W. Reagan National Defense Authorization Act of Fiscal Year 2005 (NDAA) provides the Secretary of Energy, in consultation with the Nuclear Regulatory Commission (NRC), a methodology to determine that certain waste resulting from prior reprocessing of spent nuclear fuel are not high-level radioactive waste if it can be demonstrated that the waste meets the criteria set forth in Section 3116(a) of the NDAA. The Secretary of Energy, in consultation with the NRC, signed a determination in January 2006, pursuant to Section 3116(a) of the NDAA, for salt waste disposal at the SRS Saltstone Disposal Facility. This determination is based, in part, on the Basis for Section 3116 Determination for Salt Waste Disposal at the Savannah River Site and supporting references, a document that describes the planned methods of liquid waste treatment and the resulting waste streams. The document provides descriptions of the proposed methods for processing salt waste, dividing them into 'Interim Salt Processing' and later processing through the Salt Waste Processing Facility (SWPF). Interim Salt Processing is separated into Deliquification, Dissolution, and Adjustment (DDA) and Actinide Removal Process/Caustic Side Solvent Extraction Unit (ARP/MCU). The Waste Determination was signed by the Secretary of Energy in January of 2006 based on proposed processing techniques with the expectation that it could be revised as new processing capabilities became viable. Once signed, however, it became evident that any changes would require lengthy review and another determination signed by the Secretary of Energy. With the maturation of additional salt removal technologies and the extension of the SWPF start-up date, it becomes necessary to define 'equivalency' to the processes laid out in the original determination. For the purposes of SRS, any waste not processed through Interim Salt Processing must be processed through SWPF or an equivalent process, and therefore a clear statement of the requirements for a process to be equivalent to SWPF becomes necessary. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Exact methods for self interacting neutrinos

The effective many-body Hamiltonian which describes vacuum oscillations and self interactions of neutrinos in a two flavor mixing scheme under the single angle approximation has the same dynamical symmetries as the well known BCS pairing Hamiltonian. These dynamical symmetries manifest themselves in terms of a set of constants of motion and can be useful in formulating the collective oscillation modes in an intuitive way. In particular, we show that a neutrino spectral split can be simply viewed as an avoided level crossing between the eigenstates of a mean field Hamiltonian which includes a Lagrange multiplier in order to fix the value of an exact many-body constant of motion. We show that the same dynamical symmetries also exist in the three neutrino mixing scheme by explicitly writing down the corresponding constants of motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Groundwater recharge estimation in arid hardrock-alluvium aquifers using combined water-table fluctuation and groundwater balance approaches

This article proposes an approach to estimate groundwater recharge using an optimization-based water-table fluctuation method combined with a groundwater balance model in an arid hardrock-alluvium region, located at the Oman–United Arab Emirates border. Here, we introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long-term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm 3 , respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12-year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.

58 GEOSCIENCES↗
DOE OSTI2025

A fully ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 system

A ro-vibrationally resolved corona model for the molecular hydrogen Fulcher-𝛼 transition based on the Yacora solver is introduced. The model couples 1365 ro-vibrational levels of the 𝑋 1 ⁢𝛴$^+_𝑔$, 𝑑 3⁢ 𝛱 𝑢 and 𝑎 3 ⁢𝛴$^+_𝑔$ states via electron impact excitation from the ground state into the 𝑑 3⁢ 𝛱 𝑢 state and subsequent spontaneous emission into the 𝑎 3 ⁢𝛴$+_𝑔$ state. For the process of electron impact excitation a set of 45260 fully ro-vibrationally resolved cross sections calculated with the molecular convergent close-coupling (MCCC) method in the adiabatic-nuclei formulation is applied. The MCCC cross sections are compared with results of a scaling method used in other works, demonstrating the need for dedicated ro-vibrationally resolved cross sections. By post-processing the model output, entire Fulcher spectra can be simulated. These spectra are benchmarked with measurements from an inductively coupled plasma discharge at a pressure between 1.1–10 Pa and a RF power of 700 W. The model results agree very well with the experiment both in relative shape and in absolute value. A non-ro-vibrationally resolved, purely electronically resolved collisional radiative model for molecular hydrogen (including further states and reaction channels) is applied to investigate the validity of the corona approximation for the benchmark plasmas and the relevance of possible process extensions of the corona model. Furthermore, the influence of collisional quenching, autoionization and predissociation is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗
DOE OSTI2024

IsoForma: An R Package for Quantifying and Visualizing Positional Isomers in Top-Down LC-MS/MS Data

Proteoforms, the different forms of a protein with sequence variations including post-translational modifications (PTMs), execute vital functions in biological systems such as cell signaling and epigenetic regulation. Precisely defining the stoichiometry of PTMs has been challenging because, in the widely used bottom-up proteomics methods, the detection occurs at the peptide level and thus the link between peptides and their specific modification site is lost, resulting in proteoform ambiguity. Advances in top-down mass spectrometry (MS) technology have permitted the direct characterization of intact proteoforms and their exact number of modification sites, allowing for the relative quantification of positional isomers (PI). Proteins with positional isomers refers to proteoforms with identical total mass and set of modifications but varying PTM site combinations. The relative abundance of PI can be estimated by matching proteoform-specific fragment ions to top-down tandem MS (MS2) data to localize and quantify modifications. However, current approaches heavily rely on manual annotation. Here, we present IsoForma, an open-source R package for relative quantification of PI within a single tool. We benchmarked IsoForma’s performance against two existing workflows and highlight the similarity of the results and improvements in speed. Overall, IsoForma provides a streamlined process, reduces the time of conducting isoform-based analyses, and offers an essential framework for developing customized proteoform analysis workflows. Finally, the software is open source and available at https://github.com/EMSL-Computing/isoforma-lib.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Coupled cluster evaluation of the frequency dispersion of the first and second hyperpolarizabilities of water, methanol, and dimethyl ether

The static and dynamic first (β{sub ‖}) and second (γ{sub ‖}) hyperpolarizabilities of water, methanol, and dimethyl ether have been evaluated within the response function approach using a hierarchy of coupled cluster levels of approximation and doubly augmented correlation consistent atomic basis sets. For the three compounds, the electronic β{sub ‖} and γ{sub ‖} values calculated at the CCSD and CC3 levels are in good agreement with gas phase electric field-induced second harmonic generation (EFISHG) measurements. In addition, for dimethyl ether, the frequency dispersion of both properties follows closely recent experimental values [V. W. Couling and D. P. Shelton, J. Chem. Phys. 143, 224307 (2015)] demonstrating the reliability of these methods and levels of approximation. This also suggests that the vibrational contributions to the EFISHG responses of these molecules are small.

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

Statistical anisotropy of CMB as a probe of conformal rolling scenario

Search for the statistical anisotropy in the CMB data is a powerful tool for constraining models of the early Universe. In this paper we focus on the recently proposed cosmological scenario with conformal rolling. We consider two sub-scenarios, one of which involves a long intermediate stage between conformal rolling and conventional hot epoch. Primordial scalar perturbations generated within these sub-scenarios have different direction-dependent power spectra, both characterized by a single parameter h{sup 2}. We search for the signatures of this anisotropy in the seven-year WMAP data using quadratic maximum likelihood method, first applied for similar purposes by Hanson and Lewis. We confirm the large quadrupole anisotropy detected in V and W bands, which has been argued to originate from systematic effects rather than from cosmology. We construct an estimator for the parameter h{sup 2}. In the case of the sub-scenario with the intermediate stage we set an upper limit h{sup 2} < 0.045 at the 95% confidence level. The constraint on h{sup 2} is much weaker in the case of another sub-scenario, where the intermediate stage is absent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Model Residuals as Shields: A Two-Level Formulation to Defend Smart Grids From Poisoning Attacks

The advancement of smart grids presents both vast opportunities and heightened cybersecurity risks. Data-driven defense mechanisms, though designed as a shield against these threats, can fall prey to poisoning attacks. We delve into regression settings, underscoring the imperative to fortify defenses against a spectrum of poison ratios, notably those above 0.5—an issue scarcely addressed in prior studies. Recognizing the susceptibilities of smart grids and their manipulable sensors, we exploit the very intent of poisoning attacks, compromising model accuracy, as our defense mechanism. Our proposed two-level optimization framework discerns between poisoned and authentic data based on model residuals, outperforming or matching existing methods in 72% to 77% of precision and 75% to 80% of recalls across various poisoning attacks, poison ratios, and datasets. Once the authentic data are identified, the trained model is adaptable for a variety of applications. Comprehensive evaluations on different smart grid datasets, pitted against myriad poisoning schemes, validate our methodology’s edge over existing methods. Here, we also shed light on the implications of model misspecification originating from temporal auto-correlation, a common feature in Internet of Things and smart grid data.

Adversarial machine learning (ML)↗
DOE OSTI2018

A Solar Time-based Analog Ensemble Method for Regional Solar Power Forecasting

This paper presents a new analog ensemble method for day-ahead regional photovoltaic (PV) power forecasting with hourly resolution. By utilizing open weather forecast and power measurement data, this prediction method is processed within a set of historical data with similar meteorological data (temperature and irradiance), and astronomical date (solar time and earth declination angle). Further, clustering and blending strategies are applied to improve its accuracy in regional PV forecasting. The robustness of the proposed method is demonstrated with three different numerical weather prediction models, the North American Mesoscale Forecast System, the Global Forecast System, and the Short-Range Ensemble Forecast, for both region level and single site level PV forecasts. Using real measured data, the new forecasting approach is applied to the load zone in Southeastern Massachusetts as a case study. In conclusion, the normalized root mean square error (NRMSE) has been reduced by 13.80%-61.21% when compared with three tested baselines.

14 SOLAR ENERGY↗
DOE OSTI2018

Breeding progress and preparedness for mass‐scale deployment of perennial lignocellulosic biomass crops switchgrass, miscanthus, willow and poplar

Abstract Genetic improvement through breeding is one of the key approaches to increasing biomass supply. This paper documents the breeding progress to date for four perennial biomass crops (PBCs) that have high output–input energy ratios: namely Panicum virgatum (switchgrass), species of the genera Miscanthus (miscanthus), Salix (willow) and Populus (poplar). For each crop, we report on the size of germplasm collections, the efforts to date to phenotype and genotype, the diversity available for breeding and on the scale of breeding work as indicated by number of attempted crosses. We also report on the development of faster and more precise breeding using molecular breeding techniques. Poplar is the model tree for genetic studies and is furthest ahead in terms of biological knowledge and genetic resources. Linkage maps, transgenesis and genome editing methods are now being used in commercially focused poplar breeding. These are in development in switchgrass, miscanthus and willow generating large genetic and phenotypic data sets requiring concomitant efforts in informatics to create summaries that can be accessed and used by practical breeders. Cultivars of switchgrass and miscanthus can be seed‐based synthetic populations, semihybrids or clones. Willow and poplar cultivars are commercially deployed as clones. At local and regional level, the most advanced cultivars in each crop are at technology readiness levels which could be scaled to planting rates of thousands of hectares per year in about 5 years with existing commercial developers. Investment in further development of better cultivars is subject to current market failure and the long breeding cycles. We conclude that sustained public investment in breeding plays a key role in delivering future mass‐scale deployment of PBCs.

09 BIOMASS FUELS↗
DOE OSTI2016

SU-G-JeP2-12: Quantification of 3D Geometric Distortion for 1.5T and 3T MRI Scanners Used for Radiation Therapy

Purpose: To quantify the magnitude of geometric distortion for MRI scanners and provide recommendations for MRI imaging for radiation therapy Methods: A novel phantom, QUASAR MRID3D [Modus Medical Devices Inc.], was scanned to evaluate the level of 3D geometric distortion present in five MRI scanners used for radiation therapy in our department. The phantom was scanned using the body coil with 1mm image slice thickness to acquire 3D images of the phantom body. The phantom was aligned to its geometric center for each scan, and the field of view was set to visualize the entire phantom. The dependence of distortion magnitude with distance from imaging isocenter and with magnetic field strength (1.5T and 3T) was investigated. Additionally, the characteristics of distortion for Siemens and GE machines were compared. The image distortion for each scanner was quantified in terms of mean, standard deviation (STD), maximum distortion, and skewness. Results: The 3T and 1.5T scans show a similar absolute distortion with a mean of 1.38mm (0.33mm STD) for 3T and 1.39mm (0.34mm STD) for 1.5T for a 100mm radius distance from isocenter. Some machines can have a distortion larger than 10mm at a distance of 200mm from the isocenter. The distortions are presented with plots of the x, y, and z directional components. Conclusion: The results indicate that quantification of MRI image distortion is crucial in radiation oncology for target and organ delineation and treatment planning. The magnitude of geometric distortion determines the margin needed for target contouring which is usually neglected in treatment planning process, especially for SRS/SBRT treatments. Understanding the 3D distribution of the MRI image distortion will improve the accuracy of target delineation and, hence, treatment efficacy. MRI imaging with proper patient alignment to the isocenter is vital to reducing the effects of MRI distortion in treatment planning.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

White Box Access to Quantum Testbeds for Co-Design

At Lawrence Livermore National Laboratory (LLNL), we operate and maintain the Quantum Device and Integration Testbed (QuDIT) facility, a small quantum testbed that supports about 10 active research teams (including our own) and over 50 internal and external collaborators. This testbed is designed to give remote white box access to users for research, training, and outreach. A guiding principle behind the development of our testbed infrastructure, software and user interfaces is to empower users to perform experiments at the cutting edge of quantum information science at any level of abstraction, from materials studies, device physics and control and characterization techniques to algorithm development and quantum operating system design. Our testbed targets a multilevel quantum system (qudit) to expand the accessible Hilbert space of a simple-to-manufacture quantum device and focuses on quantum simulation, typically implemented through custom gates designed with quantum optimal control methods, rather than on a universal computing framework with a fixed gate set. We leverage the Lab’s high-performance computing (HPC) program and related expertise to simulate quantum systems, develop hybrid algorithms, and generate gates optimized for given simulations. Additionally, we have adopted a co-design philosophy from the HPC community in designing new hardware, so that the systems we develop are optimized for the specific physics simulations we plan to use them for.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A Two-Fold Structural Classification Method for Determining the Accurate Ensemble of Protein Structures

Atomic-level structural characterization of flexible proteins, such as intrinsically disordered proteins and multi-domain proteins connected by flexible linkers, is challenging as they possess distinct conformations in physiological conditions. Significant efforts have been made to develop integrated approaches by combining small angle neutron/X-ray scattering experiments with molecular simulations to reveal the distinct atomic structures and the corresponding populations for these flexible proteins. One widely used method, the basis-set supported ensemble method, classifies the simulation-generated protein conformations into a set of structural basis and then derives the corresponding populations by fitting to the experimental data. This method makes an implicit assumption that protein conformations of similar structures have similar small angle scattering profiles.The present work demonstrates that, for various protein systems ranging from compact globular proteins and flexible multidomain proteins through to intrinsically disordered proteins, this method provides inaccurate assessment of the structural ensemble of the protein molecules due to the breakdown of the assumption made. To alleviate this problem, a two-fold-clustering method is developed to cluster the simulation-generated protein structures using information on both 3D structure and scattering profiles. As benchmarked by both simulation and experimental results, this new method yields much more accurate populations of structural basis of protein molecules.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Hypersonic blunt body wake computations using DSMC and Navier-Stokes solvers

Numerical results obtained with direct simulation Monte Carlo (DSMC) and Navier-Stokes methods are presented for Mach 20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how does the wake structure change as a function of rarefaction, what are the afterbody levels of heating, and to what limits are continuum models realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Moss, James N.↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

WAKE STRUCTURE↗
NASA NTRS1998

Bounded Error Schemes for the Wave Equation on Complex Domains

This paper considers the application of the method of boundary penalty terms ("SAT") to the numerical solution of the wave equation on complex shapes with Dirichlet boundary conditions. A theory is developed, in a semi-discrete setting, that allows the use of a Cartesian grid on complex geometries, yet maintains the order of accuracy with only a linear temporal error-bound. A numerical example, involving the solution of Maxwell's equations inside a 2-D circular wave-guide demonstrates the efficacy of this method in comparison to others (e.g. the staggered Yee scheme) - we achieve a decrease of two orders of magnitude in the level of the L2-error.

Abarbanel, Saul↗
NASA NTRS2004

Inner Radiation Belt Representation of the Energetic Electron Environment: Model and Data Synthesis Using the Salammbo Radiation Belt Transport Code and Los Alamos Geosynchronous and GPS Energetic Particle Data

The highly energetic electron environment in the inner magnetosphere (GEO inward) has received a lot of research attention in resent years, as the dynamics of relativistic electron acceleration and transport are not yet fully understood. These electrons can cause deep dielectric charging in any space hardware in the MEO to GEO region. We use a new and novel approach to obtain a global representation of the inner magnetospheric energetic electron environment, which can reproduce the absolute environment (flux) for any spacecraft orbit in that region to within a factor of 2 for the energy range of 100 KeV to 5 MeV electrons, for any levels of magnetospheric activity. We combine the extensive set of inner magnetospheric energetic electron observations available at Los Alamos with the physics based Salammbo transport code, using the data assimilation technique of "nudging". This in effect input in-situ data into the code and allows the diffusion mechanisms in the code to interpolate the data into regions and times of no data availability. We present here details of the methods used, both in the data assimilation process and in the necessary inter-calibration of the input data used. We will present sample runs of the model/data code and compare the results to test spacecraft data not used in the data assimilation process.

Friedel, R. H. W.↗