SEARCH · Engineering Papers
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.
We report the results of the first search for the decay B$^{0}_{s}$ → π 0 π 0 using 121.4 fb -1 of data collected at the Y(5S) resonance with the Belle detector at the KEKB asymmetric-energy e + e - collider. We observe no signal and set a 90% confidence level upper limit of 7.7 ×10 -6 on the B$^{0}_{s}$ → π 0 π 0 decay branching fraction.
Particle-based simulations are ubiquitous throughout many fields of computational science and engineering, spanning the atomistic level with molecular dynamics (MD), to mesoscale particle-in-cell (PIC) simulations for solid mechanics, device-scale modeling with PIC methods for plasma physics, and massive N-body cosmology simulations of galaxy structures, with many other methods in between (Hockney & Eastwood, 1989). While these methods use particles to represent significantly different entities with completely different physical models, many low-level details are shared including performant algorithms for short- and/or long-range particle interactions, multi-node particle communication patterns, and other data management tasks such as particle sorting and neighbor list construction. Cabana is a performance portable library for particle-based simulations, developed as part of the Co-Design Center for Particle Applications (CoPA) within the Exascale Computing Project (ECP) (Alexander et al., 2020). The CoPA project and its full development scope, including ECP partner applications, algorithm development, and similar software libraries for quantum MD, is described in (Mniszewski et al., 2021). Cabana uses the Kokkos library for on-node parallelism (Edwards et al., 2014; Trott et al., 2022), enabling simulation on multi-core CPU and GPU architectures, and MPI for GPU-aware, multi-node communication. Cabana provides particle simulation capabilities on almost all current Kokkos backends, including serial execution, OpenMP (including OpenMP-Target for GPUs), CUDA (NVIDIA GPUs), HIP (AMD GPUs), and SYCL (Intel GPUs), providing a clear path for the coming generation of accelerator-based exascale hardware. Cabana builds on Kokkos by providing new particle data structures and particle algorithms resulting in a similar execution policy-based, node-level programming model that is intended to be used in addition to the core Kokkos library within an application. Cabana is designed as an application and physics agnostic, but particle-specific toolkit which can either be used to generate a new application, or to be used as needed in existing applications at various levels of invasiveness including through interfaces that wrap user memory in existing data structures.
Euler-Lagrange point-particle simulation has emerged as a premier methodology for studying dispersed particle-laden flows. This method's popularity stems from its ability to resolve fine-scale fluid structures while also tracking individual particles at reduced cost using an appropriate particle acceleration model. However, the point-particle model has known convergence issues in that refinement of the fluid grid can lead to changes in the predicted statistics. The reasons for non-convergence are twofold: the point-particle two-way coupling force in the Navier-Stokes equations requires a numerical regularization and, without careful implementation, yields a singular force on the fluid with grid refinement. The second factor that yields grid-dependent statistics is that the point-particle force model in general depends on the undisturbed fluid velocity. When the undisturbed fluid velocity is not robustly modeled in a grid-insensitive way, the calculated force for both particles and fluid will be grid-dependent, contaminating their respective statistics. While the first issue regarding regularizing the point-particle source term has received attention in the literature, the consequences of robustly modeling the undisturbed velocity in the context of grid refinement of turbulence has received little attention. Here, we consider decaying homogeneous isotropic turbulence laden with particles at different Stokes numbers. For a given Stokes number, we systematically refine the grid and demonstrate that explicitly modeling the undisturbed fluid velocity yields relative grid insensitivity for the energy of the particle and fluid phases, as well as acceleration of the particles. We also demonstrate that an appropriately defined dissipation rate is also grid-insensitive when an undisturbed fluid velocity correction is used. In contrast, when the undisturbed fluid velocity is modeled using the conventional approach of interpolating the local fluid velocity to the particle location, we show this procedure yields divergent statistics with grid refinement. In particular, we show that higher-order interpolation of the fluid velocity in two-way coupled problems is worse than lower-order interpolation, in the absence of a correction procedure to estimate the undisturbed fluid velocity. We also examine velocity derivative statistics of the fluid phase and demonstrate that these statistics are not in general convergent even when the undisturbed fluid velocity is explicitly modeled. Collectively, the observations in this paper are used to present a philosophy on the types of questions which are answerable with point-particle methods.
Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.
Fabrication of functional polymer-based particles by crosslinking UV-curable polymer drops in mid-air and collecting crosslinked particles in a solid container, a liquid suspension, or an air flow. The particles can contain different phases in the form or layered structures that contain one to multiple cores, or structures that are blended with dissolved or emulsified smaller domains. A curing system produces ultraviolet rays that are directed onto the particles in the jet stream from one side. A reflector positioned on other side of the jet stream reflects the ultraviolet rays back onto the particles in the jet stream.
Fabrication of functional polymer-based particles by crosslinking UV-curable polymer drops in mid-air and collecting crosslinked particles in a solid container, a liquid suspension, or an air flow. The particles can contain different phases in the form or layered structures that contain one to multiple cores, or structures that are blended with dissolved or emulsified smaller domains. A curing system produces ultraviolet rays that are directed onto the particles in the jet stream from one side. A reflector positioned on other side of the jet stream reflects the ultraviolet rays back onto the particles in the jet stream.
We explored a novel framework for simulating the small-scale processes that drive the evolution of aerosol, cloud, and precipitation particles, which are a critical gap in the predictive understanding of weather and climate. Particle-based methods have emerged as an effective tool for modeling aerosol-cloud-precipitation interactions, but existing particle-based models are computationally too expensive to simulate the large domains relevant for the atmosphere or to represent the full suite of relevant processes. The lack of a comprehensive and efficient reference model is a critical bottleneck in our understanding of cloud and precipitation processes and our ability to parameterize these processes for regional- and global-scale simulations. To address this need, we explored an approach to accelerate and expand particle-based models using a new machine learning approach, graph network-based simulators (GNS). Rather than modeling the evolution of the system by numerically integrating continuity equations, the GNS represents dynamics through learned message passing. Our aim was to develop fast and accurate surrogate models for particle-based simulations. We explored applying GNS to simulate cloud droplet transport, growth, and evaporation under turbulent conditions, but we found the GNS over-smoothed the simulations. We then applied the GNS to simulate aerosol dynamics through gas condensation and found the GNS was able to reproduce the benchmark, physics-based simulation with high accuracy.
A novel technique based on machine learning is introduced to reconstruct the decays of highly Lorentz-boosted particles. Using an end-to-end deep learning strategy, the technique bypasses existing rule-based particle reconstruction methods typically used in high energy physics analyses. It uses minimally processed detector data as input and directly outputs particle properties of interest. The new technique is demonstrated for the reconstruction of the invariant mass of particles decaying in the CMS detector. The decay of a hypothetical scalar particle A into two photons, A → γ γ , is chosen as a benchmark decay. Lorentz boosts γ L = 60 – 600 are considered, ranging from regimes where both photons are resolved to those where the photons are closely merged as one object. A training method using domain continuation is introduced, enabling the invariant mass reconstruction of unresolved photon pairs in a novel way. The new technique is validated using π 0 → γ γ decays in LHC collision data.
Synthetic apolipoproteins based on native/naturally occurring homolog proteins can be prepared using solid-phase peptide synthesis approaches combined with native chemical ligation methods to create analogs of full length apolipoproteins. The chemical synthesis is expected to allow introduction of non-natural amino acids, e.g., α,α′-dialkyl amino acids, with a periodicity that encourages both helix formation and amphipathicity. Such apolipoprotein analogs are expected to encourage, in some embodiments, facile and more complete NLP formation, enabling consideration of full spectrum of nanoparticle-based biotechnology applications ranging from therapeutic sequestration and delivery to energy/biofuel production to biopolymer production.
Most methods of particle-size analysis involve particle fractionation which, is often done by filtration or the use of screens. This paper provides proof of concept testing for new in-pipe particle fractionation technologies that can be used at laboratory to industrial scales. These technologies are extensions of recent advancements in bump arrays as well as testing of new inline angled-to-the-flow slats, and advanced boycott-based Avalanche separation technologies. Experiments were conducted for three targeted particle fractionation technologies to evaluate the in-pipe removal of oversized sands in Newtonian slurries prior to the introduction of bump arrays. The mesofluidic system uses an array of staggered posts as configured in a bump array for Newtonian slurry conditions. The methods and results of these experiments are described in this paper. This information could be applied in most industrial and laboratory systems to rapidly remove oversized particles from a flowing slurry.
During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.
The paper considers methods of particle cooling applicable to beam cooling in high energy hadron colliders at the collision energy. Presently, there are two major methods of the cooling the electron cooling and stochastic cooling. The later, in application to colliders, requires exceptionally large frequency band of cooling system. Presently two methods are considered. They are the optical stochastic cooling (OSC) and the coherent electron cooling (CEC). OSC and CEC are essentially extensions of microwave stochastic cooling, operating in 1-10 GHz frequency range, to the optical frequencies enabling bands up to 30-300 THz. The OSC uses undulators as a pickup and a kicker, and an optical amplifier for signal amplification, while the CEC uses an electron beam for all these functions. We discuss major limitations, advantages and disadvantages of electron and stochastic cooling systems.
Time projection chambers are useful for high energy particle physics, nuclear physics, and astronomy. To enhance the particle detection efficiency and performance of the projection chambers functional bilayer thin film coatings based on the atomic layer deposition method are utilized. Coating material selection is based on Auger neutralization process ion induced electron emission from metallic surfaces (e.g., Mo or W) combined with a high secondary electron emission coefficient. Application of high secondary electron emission materials (e.g., MgO and CaF2) enhances the multiplication of these emitted electrons from ion induction processes. Therefore, using suitable bilayer coatings the overall TPC signal detection efficiency can be increased.
Methods and lessons learned in developing GUI-enabled radiative modeling suite for particle-based simulations including accelerated MCRT, distribution factor surrogate-models, and volumetric participating media
Not Available
We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.
Single particle cryo-electron microscopy is a vital tool for 3D characterization of protein structures. A typical workflow involves acquiring projection images of a collection of randomly oriented particles, picking and classifying individual particle projections by orientation, and finally using the individual particle projections to reconstruct a 3D map of the electron density profile. The reconstruction is challenging because of the low signal-to-noise ratio of the data, the unknown orientation of the particles, and the sparsity of data especially when dealing with flexible proteins where there may not be sufficient data corresponding to each class to obtain an accurate reconstruction using standard algorithms. In this paper we present a model-based image reconstruction technique that uses a regularized cost function to reconstruct the 3D density map by assuming known orientations for the particles. Our method casts the reconstruction as minimizing a cost function involving a novel forward model term that accounts for the contrast transfer function of the microscope, the orientation of the particles and the center of rotation offsets. We combine the forward model term with a regularizer that enforces desirable properties in the volume to be reconstructed. Using simulated data, we demonstrate how our method can significantly improve upon the typically used approach.
Double- and single-differential cross sections for inclusive charged-current 𝜈 𝜇 -nucleus scattering are reported for the kinematic domain 0 to 2 GeV/𝑐 in three-momentum transfer and 0 to 2 GeV in available energy, at a mean 𝜈 𝜇 energy of 1.86 GeV. The measurements are based on an estimated 995,760 𝜈 𝜇 charged-current (CC) interactions in the scintillator medium of the NOvA Near Detector. The subdomain populated by 2-particle-2-hole (2p2h) reactions is identified by the cross section excess relative to predictions for 𝜈 𝜇 -nucleus scattering that are constrained by a data control sample. Models for 2-particle-2-hole processes are rated by 𝜒 2 comparisons of the predicted-versus-measured 𝜈 𝜇 CC inclusive cross section over the full phase space and in the restricted subdomain. Shortfalls are observed in neutrino generator predictions obtained using the theory-based València and SuSAv2 2p2h models.