Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear scaling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2020

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗
DOE OSTI2023

Interaction of extended dislocations with nanovoid clusters

Voids of nanoscale dimensions in irradiated metals can act as obstacles to dislocation motion and cause strengthening. In this work, nanovoid strengthening and the influences of void size, void spacing and material properties, such as stacking fault energies, on dislocation bypass mechanisms are investigated using Phase Field Dislocation Dynamics, a three-dimensional mesoscale model that predicts the minimum energy pathway taken by discrete dislocations. A broad range of face centered cubic metals (copper, nickel, silver, rhodium, and platinum) and nanovoid sizes and spacings are treated, altogether spanning void size–to–dislocation stacking fault width ratios from less than unity to ten. Material γ-surfaces, calculated from ab initio methods, are input directly into the formulation. The analysis reveals that the critical bypass stress scales linearly with the linear void fraction, effective isotropic shear modulus, and ratio of the intrinsic to unstable stacking fault energies. With only a few exceptions, the critical stress is controlled by the stress required for the leading partial to impinge the voids (to move within range of the attractive image stress field of the void). When the void diameter is nearly an order of magnitude greater than the stacking fault width, the mechanism determining critical strength shifts to the stress for the dislocation to breakaway after partially cutting the void. Furthermore, this situation corresponds to that treated by line tension models and is realized here for Pt, with a sub-nanometer stacking fault width.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electrochemically induced fracture in LLZO: How the interplay between flaw density and electrostatic potential affects operability

Fracture and short circuit in the Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolyte are two key issues that prevent its adoption in battery cells. In this paper, we utilize phase-field simulations that couple electrochemistry and fracture to evaluate the maximum electric potential that LLZO electrolytes can support as a function of crack density. In the case of a single crack, we find that the applied potential at the onset of crack propagation exhibits inverse square root scaling with respect to crack length, analogous to classical fracture mechanics. Here, we further find that the short-circuit potential scales linearly with crack length. In the realistic case where the solid electrolyte contains multiple cracks, we reveal that failure fits the Weibull model. The failure distributions shift to favor failure at lower overpotentials as areal crack density increases. Furthermore, when flawless interfacial buffers are placed between the applied potential and the bulk of the electrolyte, failure is mitigated. When constant currents are applied, current focuses in near-surface flaws, leading to crack propagation and short circuit. We find that buffered samples sustain larger currents without reaching unstable overpotentials and without failing. Our findings suggest several mitigation strategies for improving the ability of LLZO to support larger currents and improve operability.

25 ENERGY STORAGE↗
DOE OSTI2024

Multiplicity dependence of ϒ production at forward rapidity in pp collisions at $\sqrt{s}$ = 13 TeV

The measurement of ϒ(1S), ϒ(2S), and ϒ(3S) yields as a function of the charged-particle multiplicity density, d N ch /d η , using the ALICE experiment at the LHC, is reported in pp collisions at √s = 13 TeV. The Υ meson yields are measured at forward rapidity (2.5 < y < 4) in the dimuon decay channel, whereas the charged-particle multiplicity is defined at central rapidity (| η | <1). Both quantities are divided by their average value in minimum bias events to compute the self-normalized quantities. The increase of the self-normalized ϒ(1S), ϒ(2S), and ϒ(3S) yields is found to be compatible with a linear scaling with the self-normalized d N ch /d η , within the uncertainties. The self-normalized yield ratios of excited-to-ground ϒ states are compatible with unity within uncertainties. Similarly, the measured double ratio of the self-normalized ϒ(1S) to the self-normalized J/ ψ yields, both measured at forward rapidity, is compatible with unity for self-normalized charged-particle multiplicities beyond one. The measurements are compared with theoretical predictions incorporating initial or final state effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Femtoscopic study of the proton-proton and proton-deuteron systems in heavy-ion collisions at the LHC

This work reports femtoscopic correlations of p – – p ($\overline{p}$ – – $\overline{p}$) and p – – d ($\overline{p}$ – – $\overline{d}$) pairs measured in Pb–Pb collisions at center-of-mass energy per nucleon $\sqrt{s_{NN}}$ = 5.02 TeV in the ALICE Collaboration. A fit to the measured proton-proton correlation functions allows one to extract the dependence of the nucleon femtoscopic radius of the particle-emitting source on the pair transverse mass (m T ) and on the average charge particle multiplicity $\langle$dN ch /dη$\rangle$ 1/3 for three centrality intervals (0–10%, 10 – – 30 %, 30 – – 50 %). In both cases, the expected power-law and linear scalings are observed, respectively. The measured p–d correlations can be described by both two- and three-body calculations, indicating that the femtoscopy observable is not sensitive to the short-distance features of the dynamics of the p-(p-n) system, due to the large inter-particle distances in Pb–Pb collisions at the LHC. Indeed, in this study, the minimum measured femtoscopic source sizes for protons and deuterons have a minimum value at 2.73$^{+0.05}_{–0.05}$ and 3.10$^{+1.04}_{–0.86}$ fm, respectively, for the 30–50% centrality collisions. Moreover, the m T -scaling obtained for the p–p and p–d systems is compatible within 1σ of the uncertainties. These findings provide new input for fundamental studies on the production of light (anti)nuclei under extreme conditions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Spectral emission characteristics near 646 nm and plasma properties of laser-induced plasma of gaseous uranium hexafluoride

We present a detailed investigation of the temporal evolution of line and continuum emissions in laser-induced breakdown spectroscopy (LIBS) of gaseous uranium hexafluoride (UF 6 ), focusing on the spectral region near 646 nm. Spectral emission features, signal-to-background ratios (SBRs) of selected uranium lines, and spectral linewidths were examined under varying UF 6 pressures (15–60 Torr) and laser pulse energies (10–60 mJ). Higher pressures and pulse energies enhanced continuum emission and reduced SBRs but did not cause significant spectral congestion. Additional studies with the use of different laser systems, including nanosecond-pulsed Nd:YAG lasers (at fundamental and various harmonics) and a femtosecond-pulsed Ti:sapphire laser, revealed long-lived plasma continua in all cases. This persistent continuum is attributed to a pseudo-continuum from overlapping molecular emissions, as its intensity scales linearly with electron number density, deviating from the expected quadratic dependence of classical ion–electron interactions for free-free and free-bound continuum emission. Based on plasma persistence time, SBR, U II/U I intensity ratios, and electronic excitation temperature, no fundamental advantage was found for femtosecond-pulsed lasers over conventional nanosecond-pulsed ones for UF 6 enrichment assay with direct LIBS measurement.

Laser-induced breakdown spectroscopy↗
DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗
DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗
DOE OSTI2022

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

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

Ground State Energy Functional with Hartree–Fock Efficiency and Chemical Accuracy

We introduce the deep post Hartree–Fock (DeePHF) method, a machine learning-based scheme for constructing accurate and transferable models for the ground-state energy of electronic structure problems. DeePHF predicts the energy difference between results of highly accurate models such as the coupled cluster method and low accuracy models such as the Hartree–Fock (HF) method, using the ground-state electronic orbitals as the input. It preserves all the symmetries of the original high accuracy model. The added computational cost is less than that of the reference HF or DFT and scales linearly with respect to system size. We examine the performance of DeePHF on organic molecular systems using publicly available data sets and obtain the state-of-art performance, particularly on large data sets.

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

Auger Dynamics in InP/ZnSe/ZnS Quantum Dots Having Pure and Doped Shells

In this study, we report the negative trion and biexciton Auger dynamics in very high-quality InP/ZnSe/ZnS quantum dots (QDs) having varying amounts of indium-based traps in the ZnSe shell. Negative trion times are determined by time-correlated photon-counting measurements on QDs that have been photoreduced with lithium triethylborohydride. We find that the Auger times vary from 280 to 425 ps and scale linearly with the total particle volume. Excess indium in the ZnSe shell gives rise to hole traps that are transiently populated following photoexcitation. Comparing stoichiometric and nonstoichiometric particles, we find that the negative trion lifetimes are independent of the presence of indium in the ZnSe shell. Biexciton dynamics are determined from transient absorption (TA) bleach recovery measurements. Absorption of two photons results in the formation of two types of biexcitons: those having both holes in the InP core (XX state) and those having one hole trapped in the ZnSe shell (XT state). XX state Auger times are measured in stoichiometric QDs, and for these particles, the Auger times are ~80 ps. Nonstoichiometric QDs show an additional long TA decay component, assigned to Auger recombination of the XT state. Here, a previous study (J. Phys. Chem. C 2021, 125, 4110–4118) found that hole trapping at indium-based shell traps results in a slow photoluminescence (PL) rise component. We find that the fraction of the XT state formed by the absorption of two photons correlates with the product of the fraction of slow PL rise and the rise time constant following one-photon excitation.

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

Effects of Tethered Polymers on Dynamics of Nanoparticles in Unentangled Polymer Melts

Polymer-tethered nanoparticles (NPs) are commonly added to a polymer matrix to improve the material properties. Critical to the fabrication and processing of such composites is the mobility of the tethered NPs. Here, we study the motion of tethered NPs in unentangled polymer melts using molecular dynamics simulations, which offer a precise control of the grafted chain length N g and the number z of grafted chains per particle. As Ng increases, there is a crossover from particle-dominated to tethered-chain-dominated terminal diffusion of NPs with the same z. The mean squared displacement of loosely tethered NPs in the case of tethered-chain-dominated terminal diffusion exhibits two subdiffusive regimes at intermediate time scales for small z. The first one at shorter time scales arises from the dynamical coupling of the particle and matrix chains, while the one at longer time scales is due to the participation of the particle in the dynamics of the tethered chains. The friction of loosely grafted chains in unentangled melts scales linearly with the total number of monomers in the chains, as the friction of individual monomers is additive in the absence of hydrodynamic coupling. As more chains are grafted to a particle, hydrodynamic interactions between grafted chains emerge. As a result, there is a nondraining layer of hydrodynamically coupled chain segments surrounding the bare particle. Outside the nondraining layer is a free-draining layer of grafted chain segments with no hydrodynamic coupling. The boundary of the two layers is the stick surface where the shear stress due to the relative melt flow is balanced by the friction between the grafted and melt chains in the interpenetration layer. The stick surface is located further away from the bare surface of the particle with higher grafting density.

42 ENGINEERING↗
DOE OSTI2024

Selected Area Manipulation of MoS 2 via Focused Electron Beam-Induced Etching for Nanoscale Device Editing

Here we demonstrate direct-write patterning of single and multilayer MoS 2 via a focused electron beam-induced etching (FEBIE) process mediated with the XeF 2 precursor. MoS 2 etching is performed at various currents, areal doses, on different substrates, and characterized using scanning electron and atomic force microscopies as well as Raman and photoluminescence spectroscopies. Scanning transmission electron microscopy reveals a sub-40 nm etching resolution and the progression of point defects and lateral etching of the consequent unsaturated bonds. The results confirm that the electron beam-induced etching process is minimally invasive to the underlying material in comparison to ion beam techniques, which damage the subsurface material. Single-layer MoS 2 field-effect transistors are fabricated, and device characteristics are compared for channels that are edited via the selected area etching process. The source–drain current at constant gate and source–drain voltage scale linearly with the edited channel width. Moreover, the mobility of the narrowest channel width decreases, suggesting that backscattered and secondary electrons collaterally affect the periphery of the removed area. Focused electron beam doses on single-layer transistors below the etching threshold were also explored as a means to modify/thin the channel layer. The FEBIE exposures showed demonstrative effects via the transistor transfer characteristics, photoluminescence spectroscopy, and Raman spectroscopy. While strategies to minimize backscattered and secondary electron interactions outside of the scanned regions require further investigation, here, we show that FEBIE is a viable approach for selective nanoscale editing of MoS 2 devices.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting Adhesion Energies of Metal Nanoparticles to Support Surfaces, Which Determines Metal Chemical Potential versus Particle Size and Thus Catalyst Performance

Improved catalysts and electrocatalysts composed of transition metal nanoparticles dispersed on high-area supports are essential for energy and environmental technologies. The chemical potential of the metal atoms in these supported nanoparticles is an important descriptor that correlates with both their catalytic activity and deactivation rate. This descriptor (μ M ) is predictably determined by the particle size and the adhesion energy per unit area at the metal/support interface (E adh ). We show here that the adhesion energies for different metals on a given support scale linearly with a simple property of the metal: for oxides, it is proportional to the metal oxophilicity, and for the carbon support, it increases linearly with metal carbophilicity (both divided by the area per metal atom). Furthermore, these relationships allow predicting E adh for other metal/support combinations, thus allowing estimation of μ M versus particle size and thereby better structure-based predictions of catalysts’ performance, which can aid in designing improved catalysts.

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

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

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

Nanoscale View of Engineered Massive Dirac Quasiparticles in Lithographic Superstructures

Massive Dirac fermions are low-energy electronic excitations characterized by a hyperbolic band dispersion. They play a central role in several emerging physical phenomena such as topological phase transitions, anomalous Hall effects, and superconductivity. This work demonstrates that massive Dirac fermions can be controllably induced by lithographically patterning superstructures of nanoscale holes in a graphene device. Their band dispersion is systematically visualized using angle-resolved photoemission spectroscopy with nanoscale spatial resolution. A linear scaling of effective mass with feature sizes is reported, underlining the Dirac nature of the superstructures. In situ electrostatic doping dramatically enhances the effective hole mass and leads to the direct observation of an electronic band gap that results in a peak-to-peak band separation of 0.64 ± 0.03 eV, which is shown via first-principles calculations to be strongly renormalized by carrier-induced screening. The methodology demonstrates band structure engineering guided by directly viewing structurally and electrically tunable massive Dirac quasiparticles in lithographic superstructures at the nanoscale.

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

Scalable freeform optimization of wide-aperture 3D metalenses by zoned discrete axisymmetry

We introduce a novel framework for design and optimization of 3D freeform metalenses that attains nearly linear scaling of computational cost with diameter, by breaking the lens into a sequence of radial “zones” with 𝑛-fold discrete axisymmetry, where 𝑛 increases with radius. This allows vastly more design freedom than imposing continuous axisymmetry, while avoiding the compromises of the locally periodic approximation (LPA) or scalar diffraction theory. Using a GPU-accelerated finite-difference time-domain (FDTD) solver in cylindrical coordinates, we perform full-wave simulation and topology optimization within each supra-wavelength zone. We validate our approach by designing millimeter and centimeter-scale, poly-achromatic, 3D freeform metalenses which outperform the state of the art. By demonstrating the scalability and resulting optical performance enabled by our “zoned discrete axisymmetry” (ZDA) and supra-wavelength domain decomposition, we highlight the potential of our framework to advance large-scale meta-optics and next-generation photonic technologies.

Sun, Mengdi [Wesleyan University]↗