Search NASA⌕ Search

SEARCH · Search NASA

Results for “level set method”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Flowing fluid electrical conductivity logging of a deep borehole during and following drilling: estimation of transmissivity, water salinity and hydraulic head of conductive zones

Flowing fluid electrical conductivity (FFEC) logging is a hydrogeologic testing method that is usually conducted in an existing borehole. However, for the 2,500-m deep COSC-1 borehole, drilled at Åre, central Sweden, it was done within the drilling period during a scheduled 1-day break, thus having a negligible impact on the drilling schedule, yet providing important information on depths of hydraulically conductive zones and their transmissivities and salinities. In this paper, we present a reanalysis of this set of data together with a new FFEC logging data set obtained soon after drilling was completed, also over a period of 1 day, but with a different pumping rate and water-level drawdown. Their joint analysis not only results in better estimates of transmissivity and salinity in the conducting fractures intercepted by the borehole, but also yields the hydraulic head values of these fractures, an important piece of information for the understanding of hydraulic structure of the subsurface. Two additional FFEC logging tests were done about 1 year later, and are used to confirm and refine this analysis. Results show that from 250 to 2,000 m depths, there are seven distinct hydraulically conductive zones with different hydraulic heads and low transmissivity values. For the final test, conducted with a much smaller water-level drawdown, inflow ceased from some of the conductive zones, confirming that their hydraulic heads are below the hydraulic head measured in the wellbore under non-pumped conditions. Lastly, the challenges accompanying 1-day FFEC logging are summarized, along with lessons learned in addressing them.

58 GEOSCIENCES↗
DOE OSTI2024

Contributions From Cloud Morphological Changes to the Interannual Shortwave Cloud Feedback Based on MODIS and ISCCP Satellite Observations

The surface temperature-mediated change in cloud properties, referred to as the cloud feedback, continues to dominate the uncertainty in climate projections. A larger number of contemporary global climate models (GCMs) project a higher degree of warming than the previous generation of GCMs. This greater projected warming has been attributed to a less negative cloud feedback in the Southern Ocean. Here, we apply a novel “double decomposition method” that employs the “cloud radiative kernel” and “cloud regime” concepts, to two data sets of satellite observations to decompose the interannual cloud feedback into contributions arising from changes within and shifts between cloud morphologies. Our results show that contributions from the latter to the cloud feedback are large for certain regimes. We then focus on interpreting how both changes within and between cloud morphologies impact the shortwave cloud optical depth feedback over the Southern Ocean in light of additional observations. Results from the former cloud morphological changes reveal the importance of the wind response to warming increases low- and mid-level cloud optical thickness in the same region. Results from the latter cloud morphological changes reveal that a general shift from thick storm-track clouds to thinner oceanic low-level clouds contributes to a positive feedback over the Southern Ocean that is offset by shifts from thinner broken clouds to thicker mid- and low-level clouds. Our novel analysis can be applied to evaluate GCMs and potentially diagnose shortcomings pertaining to their physical parameterizations of particular cloud morphologies.

58 GEOSCIENCES↗
DOE OSTI2018

Dissection of the multichannel reaction of acetylene with atomic oxygen: from the global potential energy surface to rate coefficients and branching dynamics

The O( 3 P) + C 2 H 2 reaction is the first step in acetylene oxidation. The accurate kinetic data and the understanding of the reaction dynamics is of great importance. To this end, in this study a full-dimensional global potential energy surface (PES) for the ground triplet state of the O( 3 P) + C 2 H 2 reaction is constructed based on approximately 85,000 ab initio points calculated at the level of explicitly correlated unrestricted coupled cluster single, double, and perturbative triple excitations with the explicitly correlated polarized valence triple zeta basis set (UCCSD(T)-F12b/VTZ-F12). The PES is fit using the permutation invariant polynomial-neural network (PIP-NN) approach with a total root mean square error of 0.21 kcal mol –1 . The key topographic features of the PES, including multiple potential wells and saddle points along different reaction pathways, are well represented by this fit PES. The kinetics and dynamics of the O( 3 P) + C 2 H 2 reaction are investigated using the quasi-classical trajectory (QCT) method. The calculated rate coefficients are in good agreement with experimental data over a wide temperature range, especially when the temperature is lower than 1500 K. The product branch ratio has also been determined, which indicates the H + HCCO channel as the dominant reaction pathway at 298–3000 K, accounting for 80–90% of the overall rate coefficient, in agreement with experimental observations. The dynamics of the reaction is analyzed in detail.

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

Results of the first implementation of RF phase signature matching at LANSCE

The LINAC at the Los Alamos Neutron Science Center (LANSCE) has been utilizing the Delta-t method to match the RF cavities to the design acceleration parameters since its commissioning in 1972. The differences in time-of-flight between two subsequent Beam Position and Phase Monitors (BPPMs) are measured with both accelerated and drifting beams, depending on the whether the module is set to on or off. The algorithm optimizes the module amplitude and phase via iterative measurements if the initial phase is in the vicinity of the design value. With an upgrade to a faster readout system, a scan over the whole RF cavity phase range requires relatively less time than the classical optimization procedure. The Phase Scan Signature Matching (PSSM) method provides a time-efficient method that ensures the phase selection lands on the bunching side and empowers future analyses to build module-specific models. The PSSM also utilizes a direct model to determine the correct amplitude to sub-percent level instead of using linearized matrices. Furthermore, lacking a reliable energy measurement method in the LINAC, we measure the beam phases at two downstream locations to increase the precision of energy measurements. In this letter, we also discuss the sensitivities of PSSM, error propagation, and the implementation results for the 2019 and 2020 beam cycles.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

SU-F-P-45: Clinical Experience with Radiation Dose Reduction of CT Examinations Using Iterative Reconstruction Algorithms

Purpose: Iterative reconstruction (IR) algorithms have been adopted by medical centers in the past several years. IR has a potential to substantially reduce patient dose while maintaining or improving image quality. This study characterizes dose reductions in clinical settings for CT examinations using IR. Methods: We retrospectively analyzed dose information from patients who underwent abdomen/pelvis CT examinations with and without contrast media in multiple locations of our Healthcare system. A total of 743 patients scanned with ASIR on 64 slice GE lightspeed VCTs at three sites, and 30 patients scanned with SAFIRE on a Siemens 128 slice Definition Flash in one site was retrieved. For comparison, patient data (n=291) from a GE scanner and patient data (n=61) from two Siemens scanners where filtered back-projection (FBP) was used was collected retrospectively. 30% and 10% ASIR, and SAFIRE Level 2 was used. CTDIvol, Dose-length-product (DLP), weight and height from all patients was recorded. Body mass index (BMI) was calculated accordingly. To convert CTDIvol to SSDE, AP and lateral dimensions at the mid-liver level was measured for each patient. Results: Compared with FBP, 30% ASIR reduces dose by 44.1% (SSDE: 12.19mGy vs. 21.83mGy), while 10% ASIR reduced dose by 20.6% (SSDE 17.32mGy vs. 21.83). Use of SAFIRE reduced dose by 61.4% (SSDE: 8.77mGy vs. 22.7mGy). The geometric mean for patients scanned with ASIR was larger than for patients scanned with FBP (geometric mean is 297.48 mmm vs. 284.76 mm). The same trend was observed for the Siemens scanner where SAFIRE was used (geometric mean: 316 mm with SAFIRE vs. 239 mm with FBP). Patient size differences suggest that further dose reduction is possible. Conclusion: Our data confirmed that in clinical practice IR can significantly reduce dose to patients who undergo CT examinations, while meeting diagnostic requirements for image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-527: A Novel Dynamic Multileaf Collimator Leaf-Sequencing Algorithm in Radiation Therapy

Purpose: A novel leaf-sequencing algorithm is developed for generating arbitrary beam intensity profiles in discrete levels using dynamic multileaf collimator (MLC). The efficiency of this dynamic MLC leaf-sequencing method was evaluated using external beam treatment plans delivered by intensity modulated radiation therapy technique. Methods: To qualify and validate this algorithm, integral test for the beam segment of MLC generated by the CORVUS treatment planning system was performed with clinical intensity map experiments. The treatment plans were optimized and the fluence maps for all photon beams were determined. This algorithm started with the algebraic expression for the area under the beam profile. The coefficients in the expression can be transformed into the specifications for the leaf-setting sequence. The leaf optimization procedure was then applied and analyzed for clinical relevant intensity profiles in cancer treatment. Results: The macrophysical effect of this method can be described by volumetric plan evaluation tools such as dose-volume histograms (DVHs). The DVH results are in good agreement compared to those from the CORVUS treatment planning system. Conclusion: We developed a dynamic MLC method to examine the stability of leaf speed including effects of acceleration and deceleration of leaf motion in order to make sure the stability of leaf speed did not affect the intensity profile generated. It was found that the mechanical requirements were better satisfied using this method. The Project is sponsored by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Special Analysis: Impact of Updated GSA Flow Model on E-Area Low-Level Waste Facility Groundwater Performance

In March 2018, the Savannah River National Laboratory (SRNL) notified Solid Waste Management (SWM) of new information that could potentially impact groundwater (GW) disposal limits and possibly require temporary protective measures (Crowley 2018). GW flow directions in the E-Area Low-Level Waste Facility (LLWF) Slit Trenches (STs), Engineered Trenches (ETs), and the Low Activity Waste Vault (LAWV) have notably changed in the new 2018 General Separations Area (GSA) flow model (Flach 2018a). The 2018 GSA flow model had been updated using 20 years of new hydrologic field data and model calibration methods that employed mathematical optimization software. It was later determined that the primary reason for the change in GW flow directions was due to the influence of low permeability caps that had been placed over the Old Burial Grounds and the Low Level Radioactive Waste Disposal Facility (LLRWDF). Updated flow directions in the model produce a higher degree of plume overlap for disposal units (DU’s) in the southeastern portion of E-Area than had been predicted in the 2008 E-Area LLWF Performance Assessment (PA). SWM and SRNL outlined the scope for this Special Analysis (SA) and proposed a set of interim measures (IM’s) to protect trench operations from exceeding Performance Objectives (PO’s) during the preparation of this SA, both of which were approved by the SWM Performance Assessment Review Committee (PARC) (Mooneyhan 2018).

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

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗
DOE OSTI2014

The HCO{sup +}–H{sub 2} van der Waals interaction: Potential energy and scattering

We compute the rigid-body, four-dimensional interaction potential between HCO{sup +} and H{sub 2}. The ab initio energies are obtained at the coupled-cluster single double triple level of theory, corrected for Basis Set Superposition Errors. The ab initio points are fit onto the spherical basis relevant for quantum scattering. We present elastic and rotationally inelastic coupled channels scattering between low lying rotational levels of HCO{sup +} and para-/ortho-H{sub 2}. Results are compared with similar earlier computations with He or isotropic para-H{sub 2} as the projectile. Computations agree with earlier pressure broadening measurements.

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

Rapid prototyping of arbitrary 2D and 3D wireframe DNA origami

Wireframe DNA origami assemblies can now be programmed automatically from the top-down using simple wireframe target geometries, or meshes, in 2D and 3D, using either rigid, six-helix bundle (6HB) or more compliant, two-helix bundle (DX) edges. While these assemblies have numerous applications in nanoscale materials fabrication due to their nanoscale spatial addressability and high degree of customization, no easy-to-use graphical user interface software yet exists to deploy these algorithmic approaches within a single, standalone interface. Further, top-down sequence design of 3D DX-based objects previously enabled by DAEDALUS was limited to discrete edge lengths and uniform vertex angles, limiting the scope of objects that can be designed. Here, we introduce the open-source software package ATHENA with a graphical user interface that automatically renders single-stranded DNA scaffold routing and staple strand sequences for any target wireframe DNA origami using DX or 6HB edges, including irregular, asymmetric DX-based polyhedra with variable edge lengths and vertices demonstrated experimentally, which significantly expands the set of possible 3D DNA-based assemblies that can be designed. ATHENA also enables external editing of sequences using caDNAno, demonstrated using asymmetric nanoscale positioning of gold nanoparticles, as well as providing atomic-level models for molecular dynamics, coarse-grained dynamics with oxDNA, and other computational chemistry simulation approaches.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Modeling material interfaces with hybrid adhesion method

A molecular dynamics simulation approach is presented to approximate layered material structures using discrete interatomic potentials through classical mechanics and the underlying principles of quantum mechanics. This method isolates the energetic contributions of the system into two pure material layers and an interfacial region used to simulate the adhesive properties of the diffused interface. The strength relationship of the adhesion contribution is calculated through small-scale separation calculations and applied to the molecular surfaces through an inter-layer bond criterion. By segregating the contributions into three regions and accounting for the interfacial excess energies through the adhesive surface bonds, it is possible to model each material with an independent potential while maintaining an acceptable level of accuracy in the calculation of mechanical properties. This method is intended for the atomistic study of the delamination mechanics, typically observed in thin-film applications. Therefore, the work presented in this paper focuses on mechanical tensile behaviors, with observations in the elastic modulus and the delamination failure mode. To introduce the hybrid adhesion method, we apply the approach to an ideal bulk copper sample, where an interface is created by disassociating the force potential in the middle of the structure. Various mechanical behaviors are compared to a standard EAM control model to demonstrate the adequacy of this approach in a simple setting. In addition, we demonstrate the robustness of this approach by applying it on (1) a Cu-Cu 2 O interface with interactions between two atom types, and (2) an Al-Cu interface with two dissimilar FCC lattices. These additional examples are verified against EAM and COMB control models to demonstrate the accurate simulation of failure through delamination, and the formation and propagation of dislocations under loads. Finally, the results conclude that by modeling the energy contributions of an interface using hybrid adhesion bonds, we can provide an accurate approximation method for studies of large-scale mechanical properties, as well as the representation of various delamination phenomena at the atomic scale.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Comparison of methods for geologic storage of carbon dioxide in saline formations

Preliminary estimates of CO{sub 2} storage potential in geologic formations provide critical information related to Carbon Capture, Utilization, and Storage (CCUS) technologies to mitigate CO{sub 2} emissions. Currently multiple methods to estimate CO{sub 2} storage and multiple storage estimates for saline formations have been published, leading to potential uncertainty when comparing estimates from different studies. In this work, carbon dioxide storage estimates are compared by applying several commonly used methods to general saline formation data sets to assess the impact that the choice of method has on the results. Specifically, six CO{sub 2} storage methods were applied to thirteen saline formation data sets which were based on formations across the United States with adaptations to provide the geologic inputs required by each method. Methods applied include those by (1) international efforts – the Carbon Sequestration Leadership Forum (Bachu et al., 2007); (2) United States government agencies – U.S. Department of Energy – National Energy Technology Laboratory (US-DOE-NETL, 2012) and United States Geological Survey (Brennan et al., 2010); and (3) the peer-reviewed scientific community – Szulczewski et al. (2012) and Zhou et al. (2008). A statistical analysis of the estimates generated by multiple methods revealed that assessments of CO{sub 2} storage potential made at the prospective level were often statistically indistinguishable from each other, implying that the differences in methodologies are small with respect to the uncertainties in the geologic properties of storage rock in the absence of detailed site-specific characterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Dark Energy Survey Year 3 results: Curved-sky weak lensing mass map reconstruction

ABSTRACT We present reconstructed convergence maps, mass maps, from the Dark Energy Survey (DES) third year (Y3) weak gravitational lensing data set. The mass maps are weighted projections of the density field (primarily dark matter) in the foreground of the observed galaxies. We use four reconstruction methods, each is a maximum a posteriori estimate with a different model for the prior probability of the map: Kaiser–Squires, null B-mode prior, Gaussian prior, and a sparsity prior. All methods are implemented on the celestial sphere to accommodate the large sky coverage of the DES Y3 data. We compare the methods using realistic ΛCDM simulations with mock data that are closely matched to the DES Y3 data. We quantify the performance of the methods at the map level and then apply the reconstruction methods to the DES Y3 data, performing tests for systematic error effects. The maps are compared with optical foreground cosmic-web structures and are used to evaluate the lensing signal from cosmic-void profiles. The recovered dark matter map covers the largest sky fraction of any galaxy weak lensing map to date.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Multichannel spectral estimation in acoustics: A state-space approach

Spectral estimation is a necessary methodology to analyze the frequency content of noisy data sets especially in acoustic applications. Many spectral techniques have evolved starting with the classical Fourier transform methods based on the well-known Wiener-Khintchine relationship relating the covariance-to-spectral density as a transform pair culminating with more elegant model-based parametric techniques that apply prior knowledge of the data to produce a high-resolution spectral estimate. Multichannel spectral representations are a class of both nonparametric, as well as parametric, estimators that provide improved spectral estimates. In any case, classical nonparametric multichannel techniques can provide reasonable estimates when coupled with peak-peaking methods as long as the signal levels are reasonably high. Furthermore, parametric multichannel methods can perform quite well in low signal level environments even when applying simple peak-picking techniques. In this paper, the performance of both nonparametric (periodogram) and parametric (state-space) multichannel spectral estimation methods are investigated when applied to both synthesized noisy structural vibration data as well as data obtained from a sounding rocket flight. It is demonstrated that for the multichannel problem, state-space techniques provide improved performance, offering a parametric alternative compared to classical methods.

42 ENGINEERING↗
DOE OSTI2019

Logical memory address regions

Systems, apparatuses, and methods for implementing logical memory address regions in a computing system. The physical memory address space of a computing system may be partitioned into a plurality of logical memory address regions. Each logical memory address region may be dynamically configured at run-time to meet changing application needs of the system. Each logical memory address region may also be configured separately from the other logical memory address regions. Each logical memory address region may have associated parameters that identify region start address, region size, cell-level mode, physical-to-device mapping scheme, address masks, access permissions, wear-leveling data, encryption settings, and compression settings. These parameters may be stored in a table which may be used when processing memory access requests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A geostatistics-informed hierarchical sensitivity analysis method for complex groundwater flow and transport modeling

Sensitivity analysis is an important tool for development and improvement of mathematical models, especially for complex systems with a high dimension of spatially correlated parameters. Variance-based global sensitivity analysis has gained popularity because it can quantify the relative contribution of uncertainty from different sources. However, its computational cost increases dramatically with the complexity of the considered model and the dimension of model parameters. As such, in this study, we developed a new sensitivity analysis method that integrates the concept of variance-based method with a hierarchical uncertainty quantification framework. Different uncertain inputs are grouped and organized into a multilayer framework based on their characteristics and dependency relationships to reduce the dimensionality of the sensitivity analysis. A set of new sensitivity indices are defined for the grouped inputs using the variance decomposition method. Using this methodology, we identified the most important uncertainty source for a dynamic groundwater flow and solute transport model at the Department of Energy (DOE) Hanford site. The results indicate that boundary conditions and permeability field contribute the most uncertainty to the simulated head field and tracer plume, respectively. The relative contribution from each source varied spatially and temporally. By using a geostatistical approach to reduce the number of realizations needed for the sensitivity analysis, the computational cost of implementing the developed method was reduced to a practically manageable level. The developed sensitivity analysis method is generally applicable to a wide range of hydrologic and environmental problems that deal with high-dimensional spatially distributed input variables.

58 GEOSCIENCES↗
DOE OSTI2020

In Silico Guidance for In Vitro Androgen and Glucocorticoid Receptor ToxCast Assays

Molecular initiating events (MIEs) are key events in adverse outcome pathways (AOPs) that link molecular chemistry to target biology. As they are based in chemistry, these interactions are excellent targets for computational chemistry approaches to in silico modelling. In this work, we aim to link ligand chemical structure to MIEs for androgen receptor (AR) and glucocorticoid receptor (GR) binding using ToxCast data. This has been done using an automated computational algorithm to perform maximal common substructure searches on chemical binders for each target from the ToxCast dataset. The models developed show a high level of accuracy, correctly assigning 87.20% of AR binders and 96.81% of GR binders in a 25% test set using holdout cross-validation. The 2D structural alerts developed can be used as in silico models to predict these MIEs, and as guidance for in vitro ToxCast assays to confirm hits. These models can target such experimental work, reducing the number of assays to be performed to gain required toxicological insight. Development of these models has also allowed some structural alerts to be identified as predictors for agonist or antagonist behavior at the receptor target. This work represents a first step in using computational methods to guide and target experimental approaches.

Allen, Timothy H.↗
DOE OSTI2017

A Workflow for Visual Diagnostics of Binary Classifiers using Instance-Level Explanations

In this paper, we propose a visual analytics workflow to help data scientists explore, diagnose, and understand the decisions made by a binary classifier. The approach leverages “instance-level explanations”, measures of local feature relevance that explain single instances and uses them to build a set of visual representations that guide the users in their investigation. The workflow is based on three main visual representations and steps: one based on aggregate statistics to see how data distributes across correct or incorrect decisions; one based on explanations to understand which features are used to make these decisions; and one based on raw data, to derive insights on potential root causes for the observed patterns. The work is validated through a long-term collaboration with a group of healthcare professionals who used our method to make sense of machine learning models they developed. The case studies from this collaboration demonstrate that the proposed method helps experts derive useful knowledge about the model and the phenomena it describes and also generates useful hypotheses on how a model can be improved.

Krause, Josua↗