SEARCH · Engineering Papers

Results for “Parallel processing”

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 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Towards a Circular Economy for PET Bottles in the U.S. - 4P Framework

The United States generates the most plastic waste of any country. Along with that GHG emissions from the global plastic economy are expected to increase to 15% of the global carbon budget by 2050. It is imperative that plastic recycling is made a reality to reduce both plastic pollution in the environment and GHG emissions. A portfolio of end-of-life strategies must be implemented to minimize environmental impacts and retain valuable plastic material, but it is challenging to compare options that generate products with different utility and lifetime. Plastic use reduction, reuse and recycling are thus increasingly important, but making informed policy and research decisions within this space can be challenging given the diverse range of available solutions. The novel analysis framework, Plastic Parallel Pathways Platform (4P) has been equipped with consequential life cycle assessment, techno-economic analysis, and a plastic circularity indicator to estimate the greenhouse gas (GHG) emissions, circularity, and cost of polyethylene terephthalate (PET) down-cycling to lower-quality resin, closed-loop recycling to food-grade PET bottles, up-cycling to fiber-reinforced plastic (FRP), and conversion to non-plastic products (electricity, oil) on a United States economy-wide basis. Integrating system dynamics into this robust plastics model that already incorporates techno-economics, circularity, and environmental impacts will enable identification of key bottlenecks between manufacturers, waste sorters, and reclaimers that currently prevent rapid decarbonization of the plastics economy. System dynamics (SD) explore the evolution of activities and technologies based on changed macro parameters such as plastic demand and supply, market shifts, and cross-sectoral interactions. This project particularly aims to explore the interplay between, waste collection, plastic waste sorting, recycling, and manufacturing, as well as the effect of plastic bale quality and plastic reuse initiatives on the surrounding process stages. This functionality will facilitate combinatory analysis in which a portfolio of end-of-life pathways are assessed simultaneously, with the exact makeup of that portfolio affected by parameters such as technology scales, resource constraints, and waste mitigation efforts. Integrating SD with the 4P framework enables analyzing the effect of increased revenue and reinvestment into improving process efficiencies, sorting and collection quantities. Through that, market effects of increased recycled resin availability can be studied for the plastics systems model for the US. The results will help identify technical or economic bottlenecks that currently limit efforts to decarbonize the U.S. plastics economy.

circular economy↗

DOE OSTI2024

Systems and methods of adaptive two-wavelength single-camera imaging thermography (ATSIT) for accurate and smart in-situ process temperature measurement during metal additive manufacturing

A two-wavelength, single-camera imaging thermography system for in-situ temperature measurement of a target, comprising: a target light path inlet conduit for receiving a target light beam reflected from the target; a beam splitter installed in a splitter housing at a distal end of the target light path conduit, wherein the beam splitter divides the target light beam into a first light beam and a second light beam; a first light path conduit emanating from the splitter housing comprising a first aperture iris installed within the first light path conduit for aligning the first light beam; a first band pass filter installed within the first light path conduit for regulating the first light beam to a first wavelength λ1 and an optional half waveplate installed within the first light path conduit to modulate a polarization ratio of the first light beam of λ1 wavelength; a second light path conduit emanating from the splitter housing comprising a second aperture iris installed within the second light path conduit for aligning the second light beam; a second band pass filter installed within the second light path conduit for regulating the second light beam to a second wavelength λ2; a junction housing, wherein distal ends of each of the first and second light path conduits are connected to the junction housing; a polarizing beam splitter installed in the junction housing, wherein the polarizing beam splitter reflects the first light beam of λ1 wavelength along the same path or a parallel path of the second light beam of λ2 wavelength that passes directly through the polarizing beam splitter unreflected to create a merged light beam comprising light of λ1 and λ2 wavelengths; and a light path outlet conduit connected to the junction for directing the merged beam to a high-speed camera for imaging.

Zhao, Xiayun↗

DOE OSTI2023

Separation of rare earth element radioisotopes by reverse-phase high-speed counter-current chromatography

Analytical scale purification of rare earth element (REE) radioisotopes is typically accomplished using cation-exchange resins (e.g. AG 50W-X8) and high-performance liquid chromatography (HPLC). Despite the variety of improvements made since the development of this separation process in the 1950s, nearest neighbor separations remain a challenge, as does the issue of irreversible sample adsorption. Herein, we report a study that evaluates the potential of high-speed counter-current chromatography (HSCCC) as an alternative method for purifying REE elements, with specific reference to separations of fission product REE of interest to nuclear forensics. Complementary HSCCC REE separation experiments, one spiked with radiotracer and REE fission product activity, allowed for in depth analysis of resulting fractions from both an elemental (inductively coupled plasma atomic emission spectroscopy, ICP-AES) and radiological (gamma-ray spectrometry, beta counting) purity perspective. The highly reproducible nature of separation profiles generated from HSCCC instruments was leveraged to simplify work-up of samples containing radioisotopes. Subsequent radioanalytical evaluation revealed minimal carryover of Eu into neighboring Sm and Tb fractions (as indicated by presence of 150Eu), and trace contamination of the Tb fraction with Y (as indicated by presence of 91Y). Subtle differences in stationary phase retention across the two columns were reflected in significant variations in decontamination factors of duplicate parallel separations. Furthermore, these differences paired with obtained distribution of radioisotopes provided valuable insights into future improvements. Collectively, this study represents a significant step forward in development of HSCCC technology for task specific REE radioisotope purification.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2026

The high level trigger and express data production at STAR

To meet the demands of the Beam Energy Scan phase-II (BES-II) program, the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) developed a dual real-time framework consisting of a High Level Trigger (HLT) and an Express Data Production system (xProduction). The HLT operates online within the Data Acquisition (DAQ) chain on a dedicated multi-core CPU cluster with the option to offload compute-intensive kernels to Xeon Phi coprocessors. It uses parallelized algorithms, such as the Cellular Automaton (CA) Track Finder, to perform rapid tracking, vertexing, and event filtering. This allows it to select events of interest in real time and provide immediate feedback on detector and beam conditions. In contrast, the xProduction workflow runs concurrently and independently of the DAQ loop. It applies near offline-quality calibration and reconstruction within hours of data collection. The xProduction input is the express data stream, whose content can be enriched by HLT trigger/priority selections under DAQ/HLT resource constraints, and it uses the STAR calibration/conditions framework, incorporating online calibration/QA information when available. This enables early preliminary physics analysis, including the reconstruction of rare signals, such as hyperons and hypernuclei. It also provides collaboration-wide access to analysis-ready datasets. Together, the HLT and xProduction systems form a complementary architecture: the HLT performs online event selection while the xProduction chain delivers high-quality results within a short amount of time. This integrated framework has enabled the prompt reconstruction of the $^5_Λ$ He hypernucleus with high statistical significance and the efficient processing of hundreds of millions of heavy-ion collision events. In conclusion, its demonstrated scalability and robustness establish a model for future high-luminosity experiments requiring both online event filtering and rapid access to analysis-quality data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Crystallographic Texture, Structure, and Stress Transmission in Nugget Sandstone Examined With X‐Ray Tomography and Diffraction Microscopy

Subsurface processes in sandstones are controlled by porosity, permeability, and deformation mechanisms, all of which are controlled by a complex interplay of crystallographic rock texture, structure, and micromechanics. Texture, structure, and micromechanics have historically been studied using optical and electron microscopy of thin-sections. Here, we employed a new combination of in situ X-ray tomography and ray diffraction microscopy to study crystallographic texture, structure, and grain stresses in 3D. We examined these features in a sample of Nugget sandstone, a sandstone constituting hydrocarbon reservoirs across the American West. Our aims are threefold. First, we demonstrate the utility of X-ray diffraction microscopy probes for revealing texture, structure, and stress transmission in 3D. Second, we apply these techniques to Nugget sandstone and discuss findings in the context of prior work. Third, we study grain stress tensor evolution during mechanical compression to examine whether their heterogeneity and orientation evolution reflect that of inter-particle forces in granular materials. Our results show: (a) larger grains featured higher intra-granular misorientations, possibly from an increased prevalence of cements; (b) pores closed parallel to the loading direction and opened normal to loading; (c) grain stresses featured heterogeneity and orientations similar to inter-particle forces in non-cohesive granular materials; (d) grains featured compressive stresses in the loading direction and tensile stresses orthogonal to the loading direction, the latter resisting sample dilation and grain separation. Our work demonstrates the first known application of multi-modal X-ray tomography and diffraction microscopy to sandstone, providing new 3D insight into the nature of quartz cement and stress evolution.

high energy diffraction microscopy↗

DOE OSTI2025

Fast and Accurate Greenberger-Horne-Zeilinger Encoding Using All-to-All Interactions

The 𝑁-qubit Greenberger-Horne-Zeilinger (GHZ) state is an important resource for quantum technologies. Here, we consider the task of GHZ encoding using all-to-all interactions, which prepares the GHZ state in a special case, and is furthermore useful for quantum error correction, interaction-rate enhancement, and transmitting information using power-law interactions. The naive protocol based on parallelizing CNOT gates takes O(1)-time of Hamiltonian evolution. In this work, we propose a fast protocol that achieves GHZ encoding with high accuracy. The evolution time O⁡(log 2 ⁡𝑁/𝑁) almost saturates the theoretical limit Ω⁡(log⁡𝑁/𝑁). Moreover, the final state is close to the ideal encoded one with high fidelity >1–10 −3 , up to large system sizes 𝑁 ≲ 2000. The protocol only requires a few stages of time-independent Hamiltonian evolution; the key idea is to use the data qubit as control, and to use fast spin-squeezing dynamics generated by e.g., two-axis twisting.

quantum computation↗

DOE OSTI2024

Evaluating Performance Portability with the CMS Heterogeneous Pixel Reconstruction code

In the past years the landscape of tools for expressing parallel algorithms in a portable way across various compute accelerators has continued to evolve significantly. There are many technologies on the market that provide portability between CPU, GPUs from several vendors, and in some cases even FPGAs. These technologies include C++ libraries such as Alpaka and Kokkos, compiler directives such as OpenMP, the SYCL open specification that can be implemented as a library or in a compiler, and standard C++ where the compiler is solely responsible for the offloading. Given this developing landscape, users have to choose the technology that best fits their applications and constraints. For example, in the CMS experiment the experience so far in heterogeneous reconstruction algorithms suggests that the full application contains a large number of relatively short computational kernels and memory transfer operations. In this work we use a stand-alone version of the CMS heterogeneous pixel reconstruction code as a realistic use case of HEP reconstruction software that is capable of leveraging GPUs effectively. We summarize the experience of porting this code base from CUDA to Alpaka, Kokkos, SYCL, std::par, and OpenMP offloading. We compare the event processing throughput achieved by each version on NVIDIA and AMD GPUs as well as on a CPU, and compare those to what a native version of the code achieves on each platform.

Andriotis, Nikolaos↗

DOE OSTI2024

CCSI Toolset 3.21 Release

CCSI Toolset 3.21 Release Highlights Parallelization support was added for Sequential Design of Experiments (SDOE) computations using Dask (preliminary). Input type dependent ordering capability was added to the SDOE module. With this implementation the user can specify the level of difficulty to change an input (Easy or Hard) and FOQUS will generate the appropriate ordered design depending on the input difficulty combination. Python version support was extended. FOQUS is now compatible with Python 3.8 through 3.12. Platforms used for automated testing were expanded to include macOS ARM (Apple Silicon). Updates to the FOQUS documentation to include information on how to set paths for SimSinter and TurbineLite. Turbine configuration section was added to Debugging Documentation.

AS↗

DOE OSTI2026

Experimental and simulation study of target biasing effects on plasma transport in linear plasma device MPS-LD

Linear plasma devices (LPDs) are important experimental platforms for investigating plasma–material interactions (PMI). In PMI experiments, it has been found that applying a target bias not only effectively modifies the incident ion energy, but also induces significant changes in the electron density and electron temperature, whereby the evolution of these plasma parameters is primarily governed by plasma transport processes. However, at present, the physical process and mechanism underlying such bias-induced variations remain unclear. In this work, biasing experiments under argon plasma discharge conditions were first carried out on the MPS-LD device. For the corresponding experiments, an electric potential model was newly developed based on the BOUT++ LPD module, enabling self-consistent simulations of plasma transport under biased conditions. Numerical simulations were then performed to reproduce the experimental results and to validate the accuracy of the proposed model. Finally, by combining experimental measurements with numerical simulations, a bias-voltage scan was performed to investigate how the electron density and electron temperature vary with the bias voltage (U bias ). The results show that applying negative bias decreases the target electron density (n e,T ) while increasing the target electron temperature (T e,T ). In contrast, positive bias increases both n e,T and T e,T ; however, at high positive bias, n e,T first reaches a maximum and subsequently decreases with further increases in U bias . The underlying physical mechanisms are analyzed using particle flux, momentum, and energy conservation. It indicates that the applied bias regulates the parallel electric field, thereby changing ion and electron velocities, and consequently affecting the electron density. At high positive bias, the ion velocity is further influenced by ion viscosity, leading to the reversal in n e,T . Meanwhile, the enhanced parallel electric field drives stronger currents, significantly increasing ion–electron frictional work and converting the input bias power into electron energy, which raises the electron temperature. In conclusion, these results contribute to a deeper understanding of the effects and mechanisms of biasing on plasma transport in the MPS-LD device.

BOUT++ simulation↗

DOE OSTI2020

Intrinsic Ion Transport Properties of Block Copolymer Electrolytes

Knowledge of intrinsic properties is of central importance for materials design and assessing suitability for specific applications. Self-assembling block copolymer electrolytes (BCEs) are of great interest for applications in solid-state energy storage devices. A fundamental understanding of ion transport properties, however, is hindered by the difficulty in deconvoluting extrinsic factors, such as defects, from intrinsic factors, such as the presence of interfaces between the domains. Here, we quantify the intrinsic ion transport properties of a model BCE system consisting of poly(styrene-block-ethylene oxide) (SEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt using a generalizable strategy of depositing thin films on interdigitated electrodes and self-assembling fully connected parallel lamellar structures throughout the films. Comparison between conductivity in homopolymer poly(ethylene oxide) (PEO)-LiTFSI electrolytes and the analogous conducting material in SEO over a range of salt concentrations (r, molar ratio of lithium ion to ethylene oxide repeat units) and temperatures reveals that between 20% and 50% of the PEO in SEO is inactive. Using mean-field theory calculations of the domain structure and monomer concentration profiles at domain interfaces-both of which vary substantially with salt concentration-the fraction of inactive PEO in the SEO, as derived from conductivity measurements, can be quantitatively reconciled with the fraction of PEO that is mixed with greater than a few volume percent of polystyrene. Despite the detrimental interfacial effects for ion transport in BCEs, the intrinsic conductivity of the SEO studied here (ca. 10 -3 S/cm at 90 degrees C, r = 0.085) is an order of magnitude higher than reported values from bulk samples of similar molecular weight SEO (ca. 10 -4 S/cm at 90 degrees C, r = 0.085). Finally, this work provides motivation and methods for pursuing improved BCE chemical design, interfacial engineering, and processing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

A domain wall-magnetic tunnel junction artificial synapse with notched geometry for accurate and efficient training of deep neural networks

Inspired by the parallelism and efficiency of the brain, several candidates for artificial synapse devices have been developed for neuromorphic computing, yet a nonlinear and asymmetric synaptic response curve precludes their use for backpropagation, the foundation of modern supervised learning. Spintronic devices—which benefit from high endurance, low power consumption, low latency, and CMOS compatibility—are a promising technology for memory, and domain-wall magnetic tunnel junction (DW-MTJ) devices have been shown to implement synaptic functions such as long-term potentiation and spike-timing dependent plasticity. In this work, we propose a notched DW-MTJ synapse as a candidate for supervised learning. Using micromagnetic simulations at room temperature, we show that notched synapses ensure the non-volatility of the synaptic weight and allow for highly linear, symmetric, and reproducible weight updates using either spin transfer torque (STT) or spin–orbit torque (SOT) mechanisms of DW propagation. We use lookup tables constructed from micromagnetics simulations to model the training of neural networks built with DW-MTJ synapses on both the MNIST and Fashion-MNIST image classification tasks. Accounting for thermal noise and realistic process variations, the DW-MTJ devices achieve classification accuracy close to ideal floating-point updates using both STT and SOT devices at room temperature and at 400 K. Our work establishes the basis for a magnetic artificial synapse that can eventually lead to hardware neural networks with fully spintronic matrix operations implementing machine learning.

42 ENGINEERING↗

DOE OSTI2025

Entanglement engineering of optomechanical systems by reinforcement learning

Entanglement is fundamental to quantum information science and technology, yet controlling and manipulating entanglement—so-called entanglement engineering—for arbitrary quantum systems remains a formidable challenge. There are two difficulties: the fragility of quantum entanglement and its experimental characterization. We develop a model-free deep reinforcement-learning (RL) approach to entanglement engineering, in which feedback control together with weak continuous measurement and partial state observation is exploited to generate and maintain desired entanglement. We employ quantum optomechanical systems with linear or nonlinear photon–phonon interactions to demonstrate the workings of our machine-learning-based entanglement engineering protocol. In particular, the RL agent sequentially interacts with one or multiple parallel quantum optomechanical environments, collects trajectories, and updates the policy to maximize the accumulated reward to create and stabilize quantum entanglement over an arbitrary amount of time. The machine-learning-based model-free control principle is applicable to the entanglement engineering of experimental quantum systems in general.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗

DOE OSTI2026

Subsurface mechanical damage of fused silica glass during grinding by various sub-aperture tools with and without ultrasonics

The subsurface mechanical damage (SSD) depth after grinding fused silica glass with a comprehensive set of sub-aperture fixed abrasive grinding tools [cup, wheel, belt, pad, and rotary face mill (with and without ultrasonics)] and process parameters has been statistically measured using the taper wedge technique and evaluated. Consistent with a previously reported grinding model [J. Non-Cryst. Solids 352, 5601–5617 (2006) Crossref , Materials Science and Technology of Optical Fabrication (Wiley & Sons, 2018)], based on the sliding indentation fracture by sliding particles or asperities where the normal load per particle determines the depth of the fracture (and ultimately the overall SSD depth), the dominant factor controlling SSD depth was found to be the abrasive size regardless of the tool type and process conditions. Compared to full aperture grinding methods, the overall SSD depth was higher using the sub-aperture tools, likely due to the higher effective pressure and higher load per particle distribution. Here, in addition to abrasive size, a significant reduction in SSD depth was achieved by: (1) reducing the load distribution on the abrasive particles via increase in contact area and/or decrease in mechanical loading; (2) using a more compliant host tool medium; and (3) in what we believe is a more novel way, using ultrasonics. Combining low abrasive size, larger contact area, and a compliant host, the 6 µm diamond in a resin matrix (Trizact) on a foam pad led to very low SSD depth (~ 4.6 µm), relatively fast grinding rate (186 mm 3 /h), and little or no figure degradation. This grinding tool/process is an attractive choice for final grind, resulting in significantly reduced polish out (i.e., “grey out”) times. With the rotary face mill tool, the use of ultrasonics consistently led to a SSD depth reduction (ranging from 17%–34%). A new fracture mechanics-based model, to the best of our knowledge, where the relevant normal load is parallel to the feed direction, has been developed to explain how ultrasonics leads to lower SSD depth. The key factors, supported by finite element stress analysis and load measurements, are (1) the initiation of fractures at higher z heights during the tool’s ultrasonic vertical oscillations, thus propagating less deep into workpiece; (2) reduction in load (and therefore reduction in fracture propagation distance) due to smaller tool-workpiece feed direction contact area (again caused by higher heights relative to depth of cut); (3) upward movement of the tool during oscillation leads to fracturing toward the surface instead into the depth; and finally (4) at tool’s lowest point of oscillation cycle, there may not be enough time for the fracture to propagate to its full length.

Optics and optical instruments↗

DOE OSTI2026

Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods

Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.

Axisymmetric system↗

DOE OSTI2023

Multiscale Ecosystem for solving Maxwell-Schrodinger equations of open quantum systems (OpenMS)

Light-matter interactions play an important role in many branches of physics, chemistry, energy, and materials science. In the strong coupling regime, light-matter interactions are able to tune the materials properties via the formation of new quasiparticles, such as plasmons and polaritons. However, theoretical and numerical modeling of the light-matter interaction-mediated processes remain a big challenge because light-matter interactions are fundamentally multiscale and multiphysics problems involving multiple interactions between electrons, nuclei, and photons at different time/length scales. In the current community, light-matter interactions were treated at different levels of theoretical complexity in quantum chemistry and quantum optics. In quantum optics or quantum photonics, the matter is usually simplified as a few-level system, and the light is treated quantum-mechanically. On the other hand, quantum chemistry explores first-principles methods, including both single-particle and many-body-based techniques, to describe the electronic properties of matter in detail. However, light is usually prescribed as a classical electromagnetic field, and the light-matter interaction is taken into account as an external potential via classical approximations. This software is designed to fill current modeling shortcomings by delivering a first-ever scalable multiscale platform for simulating light-matter interactions in realistic electromagnetic environments. The software solves Maxwell and Schrodinger equations self-consistent on the heterogeneous platforms. It implements HF/DFT, TDDFT, and coupled-cluster counterparts for light-matter interactions and adopts modular programming to offload massively parallel algorithms on a large number of CPU/GPUs.

Zhang, Yu↗