SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2016

TU-H-CAMPUS-JeP1-04: Deformable Image Registration Performances in Pelvis Patients: Impact of CBCT Image Quality

Purpose: To investigate the accuracy and robustness, against image noise and artifacts (typical of CBCT images), of a commercial algorithm for deformable image registration (DIR), to propagate regions of interest (ROIs) in computational phantoms based on real prostate patient images. Methods: The Anaconda DIR algorithm, implemented in RayStation was tested. Two specific Deformation Vector Fields (DVFs) were applied to the reference data set (CTref) using the ImSimQA software, obtaining two deformed CTs. For each dataset twenty-four different level of noise and/or capping artifacts were applied to simulate CBCT images. DIR was performed between CTref and each deformed CTs and CBCTs. In order to investigate the relationship between image quality parameters and the DIR results (expressed by a logit transform of the Dice Index) a bilinear regression was defined. Results: More than 550 DIR-mapped ROIs were analyzed. The Statistical analysis states that deformation strenght and artifacts were significant prognostic factors of DIR performances, while noise appeared to have a minor role in DIR process as implemented in RayStation as expected by the image similarity metric built in the registration algorithm. Capping artifacts reveals a determinant role for the accuracy of DIR results. Two optimal values for capping artifacts were found to obtain acceptable DIR results (DICE> 075/ 0.85). Various clinical CBCT acquisition protocol were reported to evaluate the significance of the study. Conclusion: This work illustrates the impact of image quality on DIR performance. Clinical issues like Adaptive Radiation Therapy (ART) and Dose Accumulation need accurate and robust DIR software. The RayStation DIR algorithm resulted robust against noise, but sensitive to image artifacts. This result highlights the need of robustness quality assurance against image noise and artifacts in the commissioning of a DIR commercial system and underlines the importance to adopt optimized protocols for CBCT image acquisitions in ART clinical implementation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Rapid QSTS Simulations for High-Resolution Comprehensive Assessment of Distributed PV

The rapid increase in penetration of distributed energy resources on the electric power distribution system has created a need for more comprehensive interconnection modeling and impact analysis. Unlike conventional scenario-based studies, quasi-static time-series (QSTS) simulations can realistically model time-dependent voltage controllers and the diversity of potential impacts that can occur at different times of year. However, to accurately model a distribution system with all its controllable devices, a yearlong simulation at 1-second resolution is often required, which could take conventional computers a computational time of 10 to 120 hours when an actual unbalanced distribution feeder is modeled. This computational burden is a clear limitation to the adoption of QSTS simulations in interconnection studies and for determining optimal control solutions for utility operations. The solutions we developed include accurate and computationally efficient QSTS methods that could be implemented in existing open-source and commercial software used by utilities and the development of methods to create high-resolution proxy data sets. This project demonstrated multiple pathways for speeding up the QSTS computation using new and innovative methods for advanced time-series analysis, faster power flow solvers, parallel processing of power flow solutions and circuit reduction. The target performance level for this project was achieved with year-long high-resolution time series solutions run in less than 5 minutes within an acceptable error.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

How Bayesian methods can improve R -matrix analyses of data: The example of the d t reaction

The 3 H(d, n) 4 He reaction is of significant interest in nuclear astrophysics and nuclear applications. It is an important, early step in big-bang nucleosynthesis and a key process in nuclear fusion reactors. We use one- and two-level R-matrix approximations to analyze data on the cross section for this reaction at center-of-mass energies below 215 keV. We critically examine the data sets using a Bayesian statistical model that allows for both common-mode and additional point-to-point un- certainties. We use Markov Chain Monte Carlo sampling to evaluate this R-matrix-plus-statistical model and find two-level R-matrix results that are stable with respect to variations in the channel radii. The S factor at 40 keV evaluates to 25.36(19) MeV b (68% credibility interval). We discuss our Bayesian analysis in detail and provide guidance for future applications of Bayesian methods to R-matrix analyses. We also discuss possible paths to further reduction of the S-factor uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

A sharp interface model for deterministic simulation of dendrite growth

A high order level set model is developed for deterministic simulation of dendritic growth in unstable solidifying systems. The model captures motion of the front implicitly on a structured finite difference grid, enables calculation of its geometric properties and also applies boundary conditions on the immersed interface. Interfacial capillary effect is incorporated in the model through the Gibbs-Thomson condition. Canonical problems for evaluating grid convergence of the numerical method and validation tests for stability of a growing nucleus in the presence of isotropic surface tension are presented. The growth morphology of solidifying nuclei in undercooled metallic melts is quantitatively analyzed. Effects of crystal anisotropy and melt undercooling on the front geometry, propagation speed and formation of branched dendritic structures are examined. In conclusion, the complex morphological changes, such as remelting of secondary perturbations under specific conditions of undercooling, are also captured during the quantitative analysis.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Sensitivity Analysis of Occupant Preferences on Energy Usage in Residential Buildings: Preprint

Residential buildings, accounting for 37% of the total electricity consumption in the United States, are suitable for demand-side management (DSM) programs to support effective and economical operation of the power system. A home energy management system (HEMS) enables residential buildings to participate in such programs. It is important to account for occupant preferences in HEMS to ensure occupant satisfaction while participating in DSM programs. For example, people who prefer a higher thermal comfort level are likely to consume more energy. In this study, we used foresee™, a HEMS developed by the National Renewable Energy Lab (NREL), to perform a sensitivity analysis of occupant preferences with the following objectives: minimize utility cost, minimize carbon footprint, and maximize thermal comfort. To incorporate the preferences into the HEMS, the SMARTER method was used to derive a set of weighting factors for each objective. We performed week-long building energy simulations using a model of a home in Fort Collins, Colorado, where there is mandatory time-of-use electricity rate structure. The foreseeTM HEMS was used to control the home with six different sets of occupant preferences. The study shows that occupant preferences can have a significant impact and is important to consider when modeling residential buildings. Results show that the HEMS could achieve energy reduction ranging from 3% to 21%, cost savings ranging from 5% to 24%, and carbon emission reduction ranging from 3% to 21%, while maintaining a low thermal discomfort level ranging from 0.78 K-hour to 6.47 K-hour in a one-week period during winter. These outcomes quantify the impact of varying occupant preferences and will be useful for controlling the electrical grid and developing HEMS solutions.

carbon footprint↗

DOE OSTI2018

Quantum Chemical Approach for Calculating Stability Constants of Mercury Complexes

Stability constants are central to the multi-scale modeling of the thermodynamic speciation, cycling, and transport of mercury (Hg) and other contaminants in aquatic environments. However, for Hg, experimental values for many relevant complexes are not available, and for others can span ranges in excess of 10 log units. The missing data and large uncertainties mean there are significant knowledge gaps in predictions of thermodynamic speciation. As an alternative to experimental measurements, thermodynamic quantities can be calculated with quantum chemical methods. Among these, density functional theory (DFT) with a polarizable continuum solvent combines accuracy with practicability. In this paper, we present an accurate and quick approach in which we use DFT with continuum solvation to calculate stability constants of Hg complexes with inorganic and low molecular-weight organic ligands in aqueous solution. Specifically, we use the M06/[SDD]6-31+G(d,p) level of theory in combination with a modified version of the SMD solvent model in which the solute radii are re-optimized with a scaled solvent-accessible surface (sSAS) approach. For the set of 37 Hg complexes used for optimization, which contain environmentally relevant functional groups and have reliable experimental stability constants, we obtain a mean unsigned error of 1.4 log units. Testing the method on an independent set of 12 Hg complexes reproduces the experimental stability constants to a mean unsigned error of 1.6 log units. Finally, this approach is a substantial step toward generally applicable rapid stability constant derivation for a wide range of Hg complexes, including those present in dissolved organic matter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

A Halotyrosine Antibody that Detects Increased Protein Modifications in Asthma Patients

Background-Airway inflammation plays an important pathophysiological role in asthma. Eosinophils produce hypobromite and bromotyrosine while neutrophils produce hypochlorite and chlorotyrosine. Objective-To evaluate halotyrosine modifications of individual airway proteins as a marker of inflammation in asthma using an antibody-based assay. Methods-We developed a novel monoclonal antibody (BTK-94C) that binds halogenated tyrosine residues, and used this antibody in a custom enzyme-linked immunosorbent assay (ELISA) microarray platform to examine halotyrosine levels in 23 proteins in three independent sets of sputum samples (52 samples total). Results-In 15 subjects with either no asthma, or with asthma characterized by high or low sputum eosinophil counts, we found associations between increased halotyrosine levels of at least three proteins and severity of airway hyperresponsiveness (AHR). Treatment with mepolizumab in 17 patients with sputum eosinophilia markedly reduced the sputum eosinophilia and significantly reduced halotyrosine levels in one sputum protein. Further analysis of 10 subjects with neutrophilic asthma and 10 health controls demonstrated a broad increase in halotyrosine in the patients with airway neutrophilia. Conclusions-Significantly higher levels of halotyrosine are associated with asthma in the asthma phenotypes we examined. The halotyrosine levels correlated with indirect AHR in the form of exercise-induced bronchoconstriction. Clinical Implication-An antibody-based assay for tyrosine halogenation in specific proteins may prove useful for assessing airway inflammation in asthma. Capsule Summary-An antibody to measure protein monobrominated tyrosine and other halotyrosine modifications was developed and used to evaluate halogenation in specific proteins in the airways for the first time. Associations were found between levels of halotyrosine and exercise-induced bronchoconstriction, and eosinophil and neutrophil inflammation in sputum from asthma patients. The ELISA developed here can be easily translated into a clinical test.

Asthma, sputum, halotyrosine, airway hyperresponsi↗

DOE OSTI2015

Depth-resolved registration of transesophageal echo to x-ray fluoroscopy using an inverse geometry fluoroscopy system

Purpose: Image registration between standard x-ray fluoroscopy and transesophageal echocardiography (TEE) has recently been proposed. Scanning-beam digital x-ray (SBDX) is an inverse geometry fluoroscopy system designed for cardiac procedures. This study presents a method for 3D registration of SBDX and TEE images based on the tomosynthesis and 3D tracking capabilities of SBDX. Methods: The registration algorithm utilizes the stack of tomosynthetic planes produced by the SBDX system to estimate the physical 3D coordinates of salient key-points on the TEE probe. The key-points are used to arrive at an initial estimate of the probe pose, which is then refined using a 2D/3D registration method adapted for inverse geometry fluoroscopy. A phantom study was conducted to evaluate probe pose estimation accuracy relative to the ground truth, as defined by a set of coregistered fiducial markers. This experiment was conducted with varying probe poses and levels of signal difference-to-noise ratio (SDNR). Additional phantom and in vivo studies were performed to evaluate the correspondence of catheter tip positions in TEE and x-ray images following registration of the two modalities. Results: Target registration error (TRE) was used to characterize both pose estimation and registration accuracy. In the study of pose estimation accuracy, successful pose estimates (3D TRE < 5.0 mm) were obtained in 97% of cases when the SDNR was 5.9 or higher in seven out of eight poses. Under these conditions, 3D TRE was 2.32 ± 1.88 mm, and 2D (projection) TRE was 1.61 ± 1.36 mm. Probe localization error along the source-detector axis was 0.87 ± 1.31 mm. For the in vivo experiments, mean 3D TRE ranged from 2.6 to 4.6 mm and mean 2D TRE ranged from 1.1 to 1.6 mm. Anatomy extracted from the echo images appeared well aligned when projected onto the SBDX images. Conclusions: Full 6 DOF image registration between SBDX and TEE is feasible and accurate to within 5 mm. Future studies will focus on real-time implementation and application-specific analysis.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Digital Application Specific Integrated Circuit Design

Single radiation effects can be transient and static. The charge collected from radiation has a prompt (drift) and a slow (diffusion) component. This current spike causes disparate results: SEU, SET, SEB, MBU, and more. A system level SEE mitigation technique, used to overcome, is the triple module redundancy (TMR). TMR for system Verilog (TMR4SV), written in Python, is the first step towards design IP triplication. It allows various triplication methods for flip flops, clocks, and logic.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Protein Conformational States—A First Principles Bayesian Method

Automated identification of protein conformational states from simulation of an ensemble of structures is a hard problem because it requires teaching a computer to recognize shapes. We adapt the naïve Bayes classifier from the machine learning community for use on atom-to-atom pairwise contacts. The result is an unsupervised learning algorithm that samples a ‘distribution’ over potential classification schemes. We apply the classifier to a series of test structures and one real protein, showing that it identifies the conformational transition with >95% accuracy in most cases. A nontrivial feature of our adaptation is a new connection to information entropy that allows us to vary the level of structural detail without spoiling the categorization. This is confirmed by comparing results as the number of atoms and time-samples are varied over 1.5 orders of magnitude. Further, the method’s derivation from Bayesian analysis on the set of inter-atomic contacts makes it easy to understand and extend to more complex cases.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Developing Performance Cost Index Targets for ASHRAE Standard 90.1 Appendix G – Performance Rating Method

Appendix G, the Performance Rating Method in ASHRAE Standard 90.1 has been updated to make two significant changes for the 2016 edition, to be published in October of 2016. First, it allows Appendix G to be used as a third path for compliance with the standard in addition to rating beyond code building performance. This prevents modelers from having to develop separate building models for code compliance and beyond code programs. Using this new version of Appendix G to show compliance with the 2016 edition of the standard, the proposed building design needs to have a performance cost index (PCI) less than targets shown in a new table based on building type and climate zone. The second change is that the baseline design is now fixed at a stable level of performance set approximately equal to the 2004 code. Rather than changing the stringency of the baseline with each subsequent edition of the standard, compliance with new editions will simply require a reduced PCI (a PCI of zero is a net-zero building). Using this approach, buildings of any era can be rated using the same method. The intent is that any building energy code or beyond code program can use this methodology and merely set the appropriate PCI target for their needs. This report discusses the process used to set performance criteria for compliance with ASHRAE Standard 90.1-2016 and suggests a method for demonstrating compliance with other codes and beyond code programs.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2016

Developing Performance Cost Index Targets for ASHRAE Standard 90.1 Appendix G – Performance Rating Method - Rev.1

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2018

Machine learning properties of binary wurtzite superlattices

The burgeoning paradigm of high-throughput computations and materials informatics brings new opportunities in terms of targeted materials design and discovery. The discovery process can be significantly accelerated and streamlined if one can learn effectively from available knowledge and past data to predict materials properties efficiently. Indeed, a very active area in materials science research is to develop machine learning based methods that can deliver automated and cross-validated predictive models using either already available materials data or new data generated in a targeted manner. In the present paper, we show that fast and accurate predictions of a wide range of properties of binary wurtzite superlattices, formed by a diverse set of chemistries, can be made by employing state-of-the-art statistical learning methods trained on quantum mechanical computations in combination with a judiciously chosen numerical representation to encode materials’ similarity. These surrogate learning models then allow for efficient screening of vast chemical spaces by providing instant predictions of the targeted properties. Moreover, the models can be systematically improved in an adaptive manner, incorporate properties computed at different levels of fidelities and are naturally amenable to inverse materials design strategies. Finally, while the learning approach to make predictions for a wide range of properties (including structural, elastic and electronic properties) is demonstrated here for a specific example set containing more than 1200 binary wurtzite superlattices, the adopted framework is equally applicable to other classes of materials as well.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Data-driven approach to parameterize SCAN + U for an accurate description of 3 d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the DFT + U method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. Here we develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-U correction, and (3) DFT/DFT + U compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for 3d transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve SCAN/SCAN + U compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize DFT + U performance for other target physical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

On the applicability of one- and many-electron quantum chemistry models for hydrated electron clusters

We evaluate the applicability of a hierarchy of quantum models in characterizing the binding energy of excess electrons to water clusters. In particular, we calculate the vertical detachment energy of an excess electron from water cluster anions with methods that include one-electron pseudopotential calculations, density functional theory (DFT) based calculations, and ab initio quantum chemistry using MP2 and eom-EA-CCSD levels of theory. The examined clusters range from the smallest cluster size (n = 2) up to nearly nanosize clusters with n = 1000 molecules. The examined cluster configurations are extracted from mixed quantum-classical molecular dynamics trajectories of cluster anions with n = 1000 water molecules using two different one-electron pseudopotenial models. We find that while MP2 calculations with large diffuse basis set provide a reasonable description for the hydrated electron system, DFT methods should be used with precaution and only after careful benchmarking. Strictly tested one-electron psudopotentials can still be considered as reasonable alternatives to DFT methods, especially in large systems. The results of quantum chemistry calculations performed on configurations, that represent possible excess electron binding motifs in the clusters, appear to be consistent with the results using a cavity structure preferring one-electron pseudopotential for the hydrated electron, while they are in sharp disagreement with the structural predictions of a non-cavity model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Hydrogen Bond Benchmark: Focal‐Point Analysis and Assessment of DFT Functionals

We performed a hierarchical, convergent ab initio benchmark study and systematically analyzed the performance of density functional approximations for describing hydrogen bonds in small neutral, cationic, and anionic complexes, as well as in larger systems involving amide, urea, deltamide, and squaramide moieties. Focal point analyses (FPA), extrapolating to the ab initio limit, were carried out using correlated wave function methods up to CCSDT(Q) for the small complexes and CCSD(T) for the larger systems, together with correlation-consistent Gaussian basis sets up to the complete basis set limit. Optimized geometries and vibrational frequencies were obtained at the CCSD(T) level. The resulting FPA hydrogen-bond energies converge within a few tenths of a kcal mol −1 . These reference data were used to evaluate 60 density functionals (including 12 dispersion-corrected), spanning the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. Overall, the meta-hybrid M06-2X provides the best performance for both hydrogen bond energies and geometries, while the dispersion-corrected GGAs BLYP-D3(BJ) and BLYP-D4 also yield accurate hydrogen-bond data and can serve as cost-effective options for studying large and complex systems.

coupled cluster theory↗