Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Schwarz 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 55 records · Page 3

Artifact identification in X-ray diffraction data using machine learning methods

In situ synchrotron high-energy X-ray powder diffraction (XRD) is highly utilized by researchers to analyze the crystallographic structures of materials in functional devices ( e.g. battery materials) or in complex sample environments ( e.g. diamond anvil cells or syntheses reactors). An atomic structure of a material can be identified by its diffraction pattern along with a detailed analysis of the Rietveld refinement which yields rich information on the structure and the material, such as crystallite size, microstrain and defects. For in situ experiments, a series of XRD images is usually collected on the same sample under different conditions ( e.g. adiabatic conditions) yielding different states of matter, or is simply collected continuously as a function of time to track the change of a sample during a chemical or physical process. In situ experiments are usually performed with area detectors and collect images composed of diffraction patterns. For an ideal powder, the diffraction pattern should be a series of concentric Debye–Scherrer rings with evenly distributed intensities in each ring. For a realistic sample, one may observe different characteristics other than the typical ring pattern, such as textures or preferred orientations and single-crystal diffraction spots. Textures or preferred orientations usually have several parts of a ring that are more intense than the rest, whereas single-crystal diffraction spots are localized intense spots owing to diffraction of large crystals, typically >10 µm. In this work, an investigation of machine learning methods is presented for fast and reliable identification and separation of the single-crystal diffraction spots in XRD images. The exclusion of artifacts during an XRD image integration process allows a precise analysis of the powder diffraction rings of interest. When it is trained with small subsets of highly diverse datasets, the gradient boosting method can consistently produce high-accuracy results. The method dramatically decreases the amount of time spent identifying and separating single-crystal diffraction spots in comparison with the conventional method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cambium 2022 Scenario Descriptions and Documentation

The National Renewable Energy Laboratory’s (NREL’s) Cambium data sets are annually released sets of simulated hourly emission, cost, and operational data for a range of modeled futures of the U.S. electric sector with metrics designed to be useful for long-term decision making. The 2022 Cambium data set is the third annual release. The data sets are a companion product to NREL’s Standard Scenarios, which are likewise released annually and are a set of projections of how the U.S. electric sector could evolve across a suite of different potential futures, but covering more scenarios with less temporal granularity. In this documentation, we describe Cambium 2022’s scenarios (Section 3), define the metrics (Section 5), and document the Cambium-specific methods for calculating those metrics (Section 6).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

Cambium 2023 Scenario Descriptions and Documentation

The National Renewable Energy Laboratory's (NREL's) Cambium data sets are annually released sets of simulated hourly emission, cost, and operational data for a range of modeled futures of the U.S. electric sector with metrics designed to be useful for long-term decision- making. The 2023 Cambium data set is the fourth annual release. The data sets are a companion product to NREL's Standard Scenarios, which are likewise released annually and are a set of projections of how the U.S. electric sector could evolve across a suite of different potential futures, but covering more scenarios with less temporal granularity (Gagnon et al. 2024). Information about Cambium and related publications can be found at https://www.nrel.gov/analysis/cambium.html, and the Cambium data sets can be viewed and downloaded at https://scenarioviewer.nrel.gov/. In this documentation, we describe Cambium 2023's scenarios, define the metrics, and document the Cambium-specific methods for calculating those metrics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thiol-based defect healing of WSe 2 and WS 2

Recent research on two-dimensional (2D) transition metal dichalcogenides (TMDCs) has led to remarkable discoveries of fundamental phenomena and to device applications with technological potential. Large-scale TMDCs grown by chemical vapor deposition (CVD) are now available at continuously improving quality, but native defects and natural degradation in these materials still present significant challenges. Spectral hysteresis in gate-biased photoluminescence (PL) measurements of WSe 2 further revealed long-term trapping issues of charge carriers in intrinsic defect states. To address these issues, we apply here a two-step treatment with organic molecules, demonstrating the “healing” of native defects in CVD-grown WSe 2 and WS 2 by substituting atomic sulfur into chalcogen vacancies. We uncover that the adsorption of thiols provides only partial defect passivation, even for high adsorption quality, and that thiol adsorption is fundamentally limited in eliminating charge traps. However, as soon as the molecular backbone is trimmed and atomic sulfur is released to the crystal, both bonds of the sulfur are recruited to passivate the divalent defect and the semiconductor quality improves drastically. Time-dependent X-ray photoelectron spectroscopy (XPS) is applied here together with other methods for the characterization of defects, their healing, leading energies and occupation. First-principles calculations support a unified picture of the electronic passivation of sulfur-healed WSe 2 and WS 2 . This work provides a simple and efficient method for improving the quality of 2D semiconductors and has the potential to impact device performance even after natural degradation.

36 MATERIALS SCIENCE↗

A Fast Temporal Decomposition Procedure for Long-Horizon Nonlinear Dynamic Programming

We propose a fast temporal decomposition procedure for solving long-horizon nonlinear dynamic programs. The core of the procedure is sequential quadratic programming (SQP) that utilizes a differentiable exact augmented Lagrangian as the merit function. Within each SQP iteration, we approximately solve the Newton system using an overlapping temporal decomposition strategy. We show that the approximate search direction is still a descent direction of the augmented Lagrangian provided the overlap size and penalty parameters are suitably chosen, which allows us to establish the global convergence. Moreover, we show that a unit step size is accepted locally for the approximate search direction and further establish a uniform, local linear convergence over stages. This local convergence rate matches the rate of the recent Schwarz scheme (Na et al. 2022). However, the Schwarz scheme has to solve nonlinear subproblems to optimality in each iteration, whereas we only perform a single Newton step instead. Numerical experiments validate our theories and demonstrate the superiority of our method.

97 MATHEMATICS AND COMPUTING↗

Rapid detection of rare events from in situ X-ray diffraction data using machine learning

High-energy X-ray diffraction methods can non-destructively map the 3D microstructure and associated attributes of metallic polycrystalline engineering materials in their bulk form. These methods are often combined with external stimuli such as thermo-mechanical loading to take snapshots of the evolving microstructure and attributes over time. However, the extreme data volumes and the high costs of traditional data acquisition and reduction approaches pose a barrier to quickly extracting actionable insights and improving the temporal resolution of these snapshots. This article presents a fully automated technique capable of rapidly detecting the onset of plasticity in high-energy X-ray microscopy data. The technique is computationally faster by at least 50 times than the traditional approaches and works for data sets that are up to nine times sparser than a full data set. This new technique leverages self-supervised image representation learning and clustering to transform massive data sets into compact, semantic-rich representations of visually salient characteristics ( e.g. peak shapes). These characteristics can rapidly indicate anomalous events, such as changes in diffraction peak shapes. It is anticipated that this technique will provide just-in-time actionable information to drive smarter experiments that effectively deploy multi-modal X-ray diffraction methods spanning many decades of length scales.

Zheng, Weijian↗

Helmholtz vibrations in bowed strings

We report for almost 160 years, it has been known that Helmholtz oscillations, unique to vibrating strings in bowed instruments (violin, cello, etc.), have two distinct regimes: “slip” and “stick.” During the slip regime, the force at the bow-string interaction is attributed to friction between the sliding bow hair and the vibrating string, with a friction coefficient that decreases with increasing relative velocity. Yet the hair-string interaction during the stick regime is less understood. We propose that the interaction force during the stick regime is proportional to the product of the longitudinal acoustic impedance of the bow hair to the relative bow-string velocity. We validate this hypothesis by solving the string's differential equation of motion, including an enhanced formulation to avoid parasitic high-frequency oscillations. This physical model enables us to analyze, in real time, the characteristics of the Helmholtz oscillations, including the string shape, excitation of harmonics, Schelleng ripples, and string energy, showing that the bowed string gains energy during the stick regime and loses energy during the slip regime.

36 MATERIALS SCIENCE↗

Direct geometric probe of singularities in band structure

A quantum system’s energy landscape may have points where multiple energy surfaces are degenerate and that exhibit singular geometry of the wave function manifold, with major consequences for the system’s properties. Ultracold atoms in optical lattices have been used to indirectly characterize such points in the band structure. Here we measured the non-Abelian transformation produced by transport directly through the singularities. We accelerated atoms along a quasi-momentum trajectory that enters, turns, and then exits the singularities at linear and quadratic band-touching points of a honeycomb lattice. Measurements after transport identified the topological winding numbers of these singularities to be 1 and 2, respectively. Our work introduces a distinct method for probing singularities that enables the study of non-Dirac singularities in ultracold-atom quantum simulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lasergate: A windowless gas target for enhanced laser preheat in magnetized liner inertial fusion

At the Z Facility at Sandia National Laboratories, the magnetized liner inertial fusion (MagLIF) program aims to study the inertial confinement fusion in deuterium-filled gas cells by implementing a three-step process on the fuel: premagnetization, laser preheat, and Z-pinch compression. In the laser preheat stage, the Z-Beamlet laser focuses through a thin polyimide window to enter the gas cell and heat the fusion fuel. However, it is known that the presence of the few μm thick window reduces the amount of laser energy that enters the gas and causes window material to mix into the fuel. These effects are detrimental to achieving fusion; therefore, a windowless target is desired. The Lasergate concept is designed to accomplish this by “cutting” the window and allowing the interior gas pressure to push the window material out of the beam path just before the heating laser arrives. In this work, we present the proof-of-principle experiments to evaluate a laser-cutting approach to Lasergate and explore the subsequent window and gas dynamics. Further, an experimental comparison of gas preheat with and without Lasergate gives clear indications of an energy deposition advantage using the Lasergate concept, as well as other observed and hypothesized benefits. While Lasergate was conceived with MagLIF in mind, the method is applicable to any laser or diagnostic application requiring direct line of sight to the interior of gas cell targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle-in-cell techniques for the study of space charge effects in the Advanced Cryogenic Gas Stopper

The Facility for Rare Isotope Beams, FRIB, will significantly increase the rates and varieties of rare isotopes compared to those presently available at existing rare-isotope facilities. The production of stopped beams via the current generation of linear gas stopping devices available is not well matched to the beam intensities that will be available at FRIB. Many experiments that make use of reaccelerated beams, from the ReA facility, greatly benefit from having the highest-intensity, rare-isotope beams possible. Some experiments would also benefit greatly from stopped beams with improved purity. In this paper, the intensity bottlenecks of FRIB’s newest linear gas stopper are studied via experiment and simulation in support of a future detailed conceptual design of a next-generation linear gas stopper that can provide rare-isotope beams with rates of 10^8 pps, a factor of 50 higher than the current state-of-the-art system, with improved beam purity.

43 PARTICLE ACCELERATORS↗

Identification of hadronic tau lepton decays using a deep neural network

A new algorithm is presented to discriminate reconstructed hadronic decays of tau leptons ($\tau_\mathrm{h}$) that originate from genuine tau leptons in the CMS detector against $\tau_\mathrm{h}$ candidates that originate from quark or gluon jets, electrons, or muons. The algorithm inputs information from all reconstructed particles in the vicinity of a $\tau_\mathrm{h}$ candidate and employs a deep neural network with convolutional layers to efficiently process the inputs. This algorithm leads to a significantly improved performance compared with the previously used one. For example, the efficiency for a genuine $\tau_\mathrm{h}$ to pass the discriminator against jets increases by 10-30% for a given efficiency for quark and gluon jets. Furthermore, a more efficient $\tau_\mathrm{h}$ reconstruction is introduced that incorporates additional hadronic decay modes. The superior performance of the new algorithm to discriminate against jets, electrons, and muons and the improved $\tau_\mathrm{h}$ reconstruction method are validated with LHC proton-proton collision data at $\sqrt{s} =$ 13 TeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Removal of spurious data in Bragg coherent diffraction imaging: an algorithm for automated data preprocessing

Bragg coherent diffraction imaging (BCDI) provides a powerful tool for obtaining high-resolution structural information from nanocrystalline materials. Here a BCDI sample consisting of a large number of randomly oriented nanoscale crystals is considered. Ideally, only one crystal is oriented to produce a Bragg peak on the detector. However, diffraction from other crystals often produces additional signals on the detector. Before the measured diffraction patterns can be processed into structural images, scientists routinely need to manually identify and remove the `alien' intensities from sources other than the intended crystal. With the development of modern high-coherence storage rings, such as the upgraded Advanced Photon Source (APS), the already slow process of manual preprocessing will be untenable for the large volumes of data that will be produced. An automated method of identifying and deleting alien intensities is proposed. This method exploits the fact that BCDI of a perfect crystal produces diffraction data with inversion symmetry around the Bragg peak. This approach uses the machine learning clustering method DBSCAN to distinguish between diffraction from multiple sources, and then calculates cluster size and inversion symmetry to assess whether clusters of intensity belong to desired data or alien signals. This approach can dramatically reduce the amount of time spent manually processing data, allowing BCDI data processing capabilities to keep pace with the technological advances of fourth-generation synchrotron light sources.

36 MATERIALS SCIENCE↗

Data on erosion and hydrogen fuel retention in Beryllium plasma-facing materials

ITER will use beryllium as a plasma-facing material in the main chamber, covering a total surface area of about 620 m 2 . Given the importance of beryllium erosion and co-deposition for tritium retention in ITER, significant efforts have been made to understand the behaviour of beryllium under fusion-relevant conditions with high particle and heat loads. This paper provides a comprehensive report on the state of knowledge of beryllium behaviour under fusion-relevant conditions: the erosion mechanisms and their consequences, beryllium migration in JET, fuel retention and dust generation. The paper reviews basic laboratory studies, advanced computer simulations and experience from laboratory plasma experiments in linear simulators of plasma–wall interactions and in controlled fusion devices using beryllium plasma-facing components. A critical assessment of analytical methods and simulation codes used in beryllium studies is given. The overall objective is to review the existing set of data with a broad literature survey and to identify gaps and research needs to broaden the database for ITER.

Erosion–deposition↗

Optimal Methods for Estimating Cactus Pear Biomass Using Cladode Dimensions of Morphologically Diverse Accessions

Current allometric methods for photosynthetic-stem (cladode) plants, such as cactus pear (Opuntia spp.), require refinement to be used in field settings in which diverse accessions are grown. We analysed cladode dimensional data using 14 accessions representing four species and two hybrids to quantify statistically significant morphological differences among accessions and derived cross-accession models to approximate cladode fresh weight. A Box model using cladode dimensions (e.g., length, width, thickness and diameter) and factorial combinations of these measures (e.g., length*width*thickness*diameter vs. fresh weight) resulted in the highest coefficient of determination (R 2 = 0.95 general fit) across all accessions for estimating fresh weight along with parsimony estimates using the Schwarz–Bayes Criterion (SBC), which assesses the most consistent performance on individual accessions. A Fitting-box modelling approach used the measured cladode area captured using ImageJ (R 2 = 0.93 general fit). Lastly, an Elliptical model used an elliptical approximation for the measured area and performed well over all accessions (R 2 = 0.94 general fit) while avoiding extensive manual measurements. These models meet or exceed the performance of previously published approaches when applied across morphologically diverse accessions, providing efficient tools for nondestructive estimation of cactus pear biomass under the conditions tested.

Opuntia↗

Luminosity determination using Z boson production at the CMS experiment

The measurement of Z boson production is presented as a method to determine the integrated luminosity of CMS data sets. The analysis uses proton–proton collision data, recorded by the CMS experiment at the CERN LHC in 2017 at a center-of-mass energy of 13 TeV. Events with Z bosons decaying into a pair of muons are selected. The total number of Z bosons produced in a fiducial volume is determined, together with the identification efficiencies and correlations from the same data set, in small intervals of 20 pb –1 of integrated luminosity, thus facilitating the efficiency and rate measurement as a function of time and instantaneous luminosity. Using the ratio of the efficiency-corrected numbers of Z bosons, the precisely measured integrated luminosity of one data set is used to determine the luminosity of another. For the first time, a full quantitative uncertainty analysis of the use of Z bosons for the integrated luminosity measurement is performed. The uncertainty in the extrapolation between two data sets, recorded in 2017 at low and high instantaneous luminosity, is less than 0.5%. We show that the Z boson rate measurement constitutes a precise method, complementary to traditional methods, with the potential to improve the measurement of the integrated luminosity.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗