Engineering Papers⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Noise robust detection of quantum phase transitions

Quantum computing allows for the manipulation of highly correlated states whose properties quickly go beyond the capacity of any classical method to calculate. Thus one natural problem which could lend itself to quantum advantage is the study of ground-states of condensed matter models, and the transitions between them. However, current levels of hardware noise can require extensive application of error-mitigation techniques to achieve reliable computations. In this work, we use several IBM devices to explore a finite-size spin model with multiple “phaselike” regions characterized by distinct ground-state configurations. Using preoptimized Variational Quantum Eigensolver (VQE) solutions, we demonstrate that in contrast to calculating the energy, where zero-noise extrapolation is required in order to obtain qualitatively accurate yet still unreliable results, calculations of the energy derivative, two-site spin correlation functions, and the fidelity susceptibility yield accurate behavior across multiple regions, even with minimal or no application of error-mitigation approaches. Taken together, these sets of observables could be used to identify level crossings in a simple, noise-robust manner which is agnostic to the method of ground state preparation. This work shows promising potential for near-term application to identifying quantum phase transitions, including avoided crossings and nonadiabatic conical intersections in electronic structure calculations. Published by the American Physical Society 2024

Lively, Kevin (ORCID:0000000320981494)↗
DOE OSTI2023

Searches for lepton-flavour-violating decays of the Higgs boson into $eτ$ and $μτ$ in $ \sqrt{s}$ = 13 TeV pp collisions with the ATLAS detector

This paper presents direct searches for lepton flavour violation in Higgs boson decays, H → eτ and H → μτ, performed using data collected with the ATLAS detector at the LHC. The searches are based on a data sample of proton-proton collisions at a centre-of-mass energy $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb –1 . Leptonic (τ → ℓν ℓ ν τ ) and hadronic (τ → hadrons ν τ ) decays of the τ-lepton are considered. Two background estimation techniques are employed: the MC-template method, based on data-corrected simulation samples, and the Symmetry method, based on exploiting the symmetry between electrons and muons in the Standard Model backgrounds. No significant excess of events is observed and the results are interpreted as upper limits on lepton-flavour-violating branching ratios of the Higgs boson. The observed (expected) upper limits set on the branching ratios at 95% confidence level, $\mathcal{B}$(H → eτ) < 0.20% (0.12%) and $\mathcal{B}$(H → μτ ) < 0.18% (0.09%), are obtained with the MC-template method from a simultaneous measurement of potential H → eτ and H → μτ signals. The best-fit branching ratio difference, $\mathcal{B}$(H → μτ) → $\mathcal{B}$(H → eτ), measured with the Symmetry method in the channel where the τ-lepton decays to leptons, is (0.25 ± 0.10)%, compatible with a value of zero within 2.5σ.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

On the Efficient Evaluation of the Exchange Correlation Potential on Graphics Processing Unit Clusters

The predominance of Kohn–Sham density functional theory (KS-DFT) for the theoretical treatment of large experimentally relevant systems in molecular chemistry and materials science relies primarily on the existence of efficient software implementations which are capable of leveraging the latest advances in modern high-performance computing (HPC). With recent trends in HPC leading toward increasing reliance on heterogeneous accelerator-based architectures such as graphics processing units (GPU), existing code bases must embrace these architectural advances to maintain the high levels of performance that have come to be expected for these methods. In this work, we purpose a three-level parallelism scheme for the distributed numerical integration of the exchange-correlation (XC) potential in the Gaussian basis set discretization of the Kohn–Sham equations on large computing clusters consisting of multiple GPUs per compute node. In addition, we purpose and demonstrate the efficacy of the use of batched kernels, including batched level-3 BLAS operations, in achieving high levels of performance on the GPU. We demonstrate the performance and scalability of the implementation of the purposed method in the NWChemEx software package by comparing to the existing scalable CPU XC integration in NWChem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Bench Scale Study of Coating Formulations for Abatement of Mercury Vapor and Recovery of Mercury Spills

A sealant formulation has been developed that effectively mitigates human exposure to elemental mercury vapor evolving from spilled mercury. The sealant may be dispensed using a fogging methodology (for industrial settings) or a common spray bottle (for home or school settings). In parallel, the potential for Hg recovery using strippable coatings was evaluated. Four liquid coatings were evaluated: one fogging fixative (FX2) developed at INL for the containment of surface and airborne contamination and three commercially available strippable decontamination gels (Carboset 441, Encor 449, and Stripcoat TLC Free). Bench scale mercury vapor abatement tests were conducted in an enclosure using FX2 with the addition of a mercury absorbent and a mercury reactant. Testing with FX2 showed the reactant contributed significantly to reducing the rate of mercury evaporation. The reactant had a minimal impact on the mercury evaporation rate for the strippable coatings. The potential of strippable coatings to remove/recover mercury was inferred from the adhesive and cohesive performance of sprayable dilutions of the coatings applied to porous and non-porous surfaces. Testing provided the following results: • Mercury vapor concentration was reduced by a factor of 40 using the reactant enhanced FX2 solution, as well as with dilutions of Carboset and Encor. Dilute Stripcoat was less effective at reducing Hg evaporation. • A 100× reduction in Hg evaporation rate was achieved using a coating of FX + the Hg reactant (a.k.a. FX Hg), compared to uncoated Hg. • A 2× reduction in the Hg vapor accumulation rate was achieved using FX Hg, compared to the strippable coatings. • Strippability and adhesion tests show that dilutions of Carboset and Stripcoat can be stripped from a variety of indoor surfaces. Dilute Encor was not strippable. Results suggest FX Hg and dilute Carboset will likely reduce Hg vapor concentrations, in an indoor setting, to a level at or below the Occupational Safety and Health Administration 8-hour human exposure limit (0.1 mg/m3). FX Hg dispensed via fogging would increase worker safety during remediation and demolition of industrial facilities with high levels of elemental mercury contamination. When applied via spray bottle, FX Hg may provide a safe and convenient method for the mitigation of Hg vapor in residential or commercial settings (e.g., offices, homes, schools). This approach would be most beneficial where recovery of the spilled Hg is impractical. A sprayed then stripped application of dilute Carboset may provide a safe and convenient method for recovery of small Hg spills in residential and commercial settings. Further testing is required to determine if liquid Hg sufficiently adheres to the strippable coating to be recovered in this way. If it can be dispensed as a fog, dilute Carboset could be employed to recover Hg contamination in industrial settings. Testing of the foggability of dilute Carboset is required to assess this potential.

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

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study

The metastable ring structure of the ozone 1{sup 1}A{sub 1} ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two {sup 1}A{sub 1} states. In the present work, valence correlated energies of the 1{sup 1}A{sub 1} state and the 2{sup 1}A{sub 1} state were calculated at the 1{sup 1}A{sub 1} open minimum, the 1{sup 1}A{sub 1} ring minimum, the transition state between these two minima, the minimum of the 2{sup 1}A{sub 1} state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1{sup 1}A{sub 1} state, the present calculations yield the estimates of (ring minimum—open minimum) ∼45–50 mh and (transition state—open minimum) ∼85–90 mh. For the (2{sup 1}A{sub 1}–{sup 1}A{sub 1}) excitation energy, the estimate of ∼130–170 mh is found at the open minimum and 270–310 mh at the ring minimum. At the transition state, the difference (2{sup 1}A{sub 1}–{sup 1}A{sub 1}) is found to be between 1 and 10 mh. The geometry of the transition state on the 1{sup 1}A{sub 1} surface and that of the minimum on the 2{sup 1}A{sub 1} surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

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

Calculating the grain boundary inclination of voxelated grain structures using a smoothing algorithm

We have developed a flexible method for calculating the grain boundary (GB) inclinations of voxelated grain structure data using smoothing algorithms. We compared the performance of four algorithms: the linear interpolation, Allen–Cahn, level-set, and vertex algorithms. We assessed their accuracy using 2D and 3D cases with known inclinations. The vertex algorithm provided the best balance between accuracy and efficiency for 2D structures while the linear interpolation algorithm provided the best balance for 3D structures. We compared the GB inclinations calculated using our smoothing method on a 3D high energy X-ray diffraction microscopy (HEDM) dataset to those determined by meshing the GBs. The two approaches determined similar GB plane distributions, though they varied significantly at triple junctions. In conclusion, the smoothing method was demonstrated for two sources of 3D voxelated grain structures: HEDM data and results from Monte Carlo Potts grain growth simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Active immunization using hand-push emulsification method increases the operator's risk of transcutaneous immunization

Highlights: • Researchers have the right to be told the relevant risks of immunological experiments. • Susceptible individuals to allergic diseases are suggested to be more careful when they perform immunological experiments. Setting up an animal model by using active immunization methods is a common means of studying immune-related diseases or producing antibodies with high titer and high activities. However, the security during the process of pathogen emulsification remains unclear. In a physical examination, we unexpectedly noticed high levels of angiotensin II type 1 receptor autoantibody (AT1-AA) specific to the immunizing antigen in the sera of some researchers who had participated in setting up active immunization animal models, and we were puzzled about the cause of AT1-AA production. In this study, we intended to investigate whether the emulsified antigen was the source of infection in these researchers, and if so, how to prevent it from occurring. AT1-AA was detected by advanced ELISA method. The participants presented higher levels of AT1-AA compared with non-participants of the same laboratory. This phenomenon remained that some factors during the process of rat model establishment may contribute to AT1-AA production. Animal and glove penetration studies indicated the emulsified antigen infection was attributed to neither aerosol or fur touch nor penetrating through gloves. However, AT1-AA level was largely decreased in the participants after they used an automatic emulsification device. Because of the strong permeability of the adjuvant, we speculated that emulsified antigen might get access to the unprotected skin of the participants accidentally during the immunization process. These results demonstrated that accidental contacts of emulsified antigens may infect researchers during the process of traditional hand-push emulsification, resulting in high specific autoantibody levels, which can be prevented by using appropriate tools.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Character of intermolecular interaction in pyridine-argon complex: Ab initio potential energy surface, internal dynamics, and interrelations between SAPT energy components

The pyridine-Ar (PAr) van der Waals (vdW) complex is studied using a high level ab initio method. Its structure, binding energy, and intermolecular vibrational states are determined from the analytical potential energy surface constructed from interaction energy (IE) values computed at the coupled cluster level of theory with single, double, and perturbatively included triple excitations with the augmented correlation consistent polarized valence double-ζ (aug-cc-pVDZ) basis set complemented by midbond functions. The structure of the complex at its global minimum with Ar at a distance of 3.509 Å from the pyridine plane and shifted by 0.218 Å from the center of mass towards nitrogen agrees well with the corresponding equilibrium structure derived previously from the rotational spectrum of PAr. The PAr binding energy D{sub e} of 392 cm{sup −1} is close to that of 387 cm{sup −1} calculated earlier at the same ab initio level for the prototypical benzene-Ar (BAr) complex. However, under an extension of the basis set, D{sub e} for PAr becomes slightly lower than D{sub e} for BAr. The ab initio vdW vibrational energy levels allow us to estimate the reliability of the methods for the determination of the vdW fundamentals from the rotational spectra. To disclose the character of the intermolecular interaction in PAr, the symmetry-adapted perturbation theory (SAPT) is employed for the analysis of different physical contributions to IE. It is found that SAPT components of IE can be approximately expressed in the binding region by only two of them: the exchange repulsion and dispersion energy. The total induction effect is negligible. The interrelations between various SAPT components found for PAr are fulfilled for a few other complexes involving aromatic molecules and Ar or Ne, which indicates that they are valid for all rare gas (Rg) atoms and aromatics.

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

Beyond pinball loss: Quantile methods for calibrated uncertainty quantification

Amongthemanywaysofquantifying uncertainty in a regression setting, specifying the full quantile function is attractive, as quantiles are amenable to interpretation and evaluation. A model that predicts the true conditional quantiles for each input, at all quantile levels, presents a correct and efficient representation of the underlying uncertainty. To achieve this, many current quantile-based methods focus on optimizing the pinball loss. However, this loss restricts the scope of applicable regression models, limits the ability to target many desirable properties (e.g. calibration, sharpness, centered intervals), and may produce poor conditional quantiles. In this work, we develop new quantile methods that address these shortcomings. In particular, we propose methods that can apply to any class of regression model, select an explicit balance between calibration and sharpness, optimize for calibration of centered intervals, and produce more accurate conditional quantiles. We provide a thorough experimental evaluation of our methods, which includes a high dimensional uncertainty quantification task in nuclear fusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Artificial Intelligence Benchmarking

AI benchmarking is a method for evaluating the effectiveness of an AI model using a set of standardized metrics, for example, high school-level math exams. These benchmarks and their results will enable ranking various AI models based on their effectiveness in performing a specific task.

Krishnan, Anjay [Fermilab]↗
DOE OSTI2017

Single–Column Model Simulations of Subtropical Marine Boundary–Layer Cloud Transitions Under Weakening Inversions

Results are presented of the GASS/EUCLIPSE single-column model inter-comparison study on the subtropical marine low-level cloud transition. A central goal is to establish the performance of state-of-the-art boundary-layer schemes for weather and climate models for this cloud regime, using large-eddy simulations of the same scenes as a reference. A novelty is that the comparison covers four different cases instead of one, in order to broaden the covered parameter space. Three cases are situated in the North-Eastern Pacific, while one reflects conditions in the North-Eastern Atlantic. A set of variables is considered that reflects key aspects of the transition process, making use of simple metrics to establish the model performance. Using this method some longstanding problems in low level cloud representation are identified. Considerable spread exists among models concerning the cloud amount, its vertical structure and the associated impact on radiative transfer. Here, the sign and amplitude of these biases differ somewhat per case, depending on how far the transition has progressed. After cloud breakup the ensemble median exhibits the well-known “too few too bright” problem. The boundary layer deepening rate and its state of decoupling are both underestimated, while the representation of the thin capping cloud layer appears complicated by a lack of vertical resolution. Encouragingly, some models are successful in representing the full set of variables, in particular the vertical structure and diurnal cycle of the cloud layer in transition. An intriguing result is that the median of the model ensemble performs best, inspiring a new approach in subgrid parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Non-Gaussian estimates of tensions in cosmological parameters

We discuss how to efficiently and reliably estimate the level of agreement and disagreement on parameter determinations from different experiments, fully taking into account non-Gaussianities in the parameter posteriors. We develop two families of scalable algorithms that allow us to perform this type of calculations in increasing number of dimensions and for different levels of tensions. One family of algorithms rely on kernel density estimates of posterior distributions while the other relies on machine learning modeling of the posterior distribution with normalizing flows. We showcase their effectiveness and accuracy with a set of benchmark examples and find both methods agree with each other and the true tension within 0.5σ in difficult cases and generally to 0.2σ or better. This allows us to study the level of internal agreement between different measurements of the clustering of cosmological structures from the Dark Energy Survey and their agreement with measurements of the Cosmic Microwave Background from the Planck satellite.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bias-Adjustment Methods for Future Subdaily Precipitation Extremes Consistent Across Durations

Model output from climate projections often requires bias-adjustment to compensate for systematic model errors. A bias-adjustment method for extreme precipitation intensity is proposed that preserves the scaling equation for different accumulation levels from hourly to daily, using intensity-duration-frequency (IDF) modeling. A validation is performed within a pseudo-reality setting, based on hourly precipitation from 28 regional climate model projections of the EURO-CORDEX ensemble over Belgium. The scaling-based adjustment methods improve upon previous methods, an optimal method is identified, and, analytical quantile mapping methods must be avoided due to three identified problems. The ensemble mean of the adjusted extreme precipitation intensity obeys the above-mentioned scale-invariance property, which is consistent with observed extreme intensities. We thus show that IDF modeling provides added value in the context of bias-adjustment, and, that the particular IDF model proposed balances well between accuracy and the preservation of desired properties such as scale invariance and consistency among rainfall durations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Few-shot leaf segmentation

In this work, we use few-shot learning to segment the body and vein architecture of Populus trichocarpa leaves from high-resolution scans obtained in the UC Davis common garden. Leaf and vein segmentation are formulated as separate tasks, in which convolutional neural networks (CNNs) are used to iteratively expand partial segmentations until reaching stopping criteria. Our leaf and vein segmentation approaches use just 50 and 8 manually traced images for training, respectively, and are applied to a set of 2,634 top and bottom leaf scans. We show that both methods achieve high segmentation accuracy, in some cases exceeding even human-level segmentation. The leaf and vein segmentations are subsequently used to extract 68 morphological traits using traditional open-source image processing tools, which are validated using real-world physical measurements. For a biological perspective, we perform a genome-wide association study using the "vein density" trait to discover novel genetic architectures associated with multiple physiological processes relating to leaf development and function. In addition to sharing all of the few-shot learning code, we are releasing all images, manual segmentations, model predictions, 68 extracted leaf phenotypes, and a new set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes.

Lagergren, John [Oak Ridge National Lab. (ORNL), O↗
DOE OSTI2023

Hybrid Cyber-attack Detection in Photovoltaic Farms

Here, to address the cyber-physical security in PV farms, a hybrid cyber-attack detection is proposed in this manuscript. To secure PV farms, the proposed method integrates model-based and data-driven methods by fusing the detection score at the device and system levels. First, a model-based cyber-attack detection method is developed for each PV inverter. A residual between the estimation of the Kalman filter and measurement is calculated. By leveraging the calculated residual from all inverters, a squared Mahalanobis distance is developed for device detection score generation. At the system level, a convolutional neural network (CNN) is proposed to detect cyber-attack using the waveform data at the point of common coupling (PCC) in PV farms. To improve the CNN detection accuracy, a set of well-designed features are extracted from the raw waveform data. Finally, a weighted detection score fusion method is proposed to combine device and system detection scores by using their complementary strength. The feasibility and robustness of the proposed method are validated by testing cases and a comparative experiment.

14 SOLAR ENERGY↗
DOE OSTI2015

SU-F-BRB-05: Collision Avoidance Mapping Using Consumer 3D Camera

Purpose: To develop a fast and economical method of scanning a patient’s full body contour for use in collision avoidance mapping without the use of ionizing radiation. Methods: Two consumer level 3D cameras used in electronic gaming were placed in a CT simulator room to scan a phantom patient set up in a high collision probability position. A registration pattern and computer vision algorithms were used to transform the scan into the appropriate coordinate systems. The cameras were then used to scan the surface of a gantry in the treatment vault. Each scan was converted into a polygon mesh for collision testing in a general purpose polygon interference algorithm. All clinically relevant transforms were applied to the gantry and patient support to create a map of all possible collisions. The map was then tested for accuracy by physically testing the collisions with the phantom in the vault. Results: The scanning fidelity of both the gantry and patient was sufficient to produce a collision prediction accuracy of 97.1% with 64620 geometry states tested in 11.5 s. The total scanning time including computation, transformation, and generation was 22.3 seconds. Conclusion: Our results demonstrate an economical system to generate collision avoidance maps. Future work includes testing the speed of the framework in real-time collision avoidance scenarios. Research partially supported by a grant from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Time delay lens modelling challenge

ABSTRACT In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant H0. However, published state-of-the-art analyses require of order 1 yr of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge, we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated data sets. The results in Rungs 1 and 2 show that methods that use only the point source positions tend to have lower precision ($10\!-\!20{{\ \rm per\ cent}}$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic data sets can recover H0 within the target accuracy (|A| < 2 per cent) and precision (<6 per cent per system), even in the presence of a poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

Ding, X.↗