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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum simulation of charge and exciton transfer in multi-mode models using engineered reservoirs

Quantum simulation enables studies of open-system dynamics in non-perturbative regimes by programming electronic, vibrational, and environmental interactions on comparable energy scales. Trapped ions offer this capability, combining spins, phonons, and tunable dissipation on one platform. We demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Using tailored spin-phonon interactions with reservoir engineering, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor sites and track its non-equilibrium dynamics. We continuously tune the system from the charge transfer regime to the vibrationally assisted exciton transfer regime and find that degenerate modes enhance transfer rates at large energy gaps, while non-degenerate modes activate pathways that reduce the energy-gap dependence. Thus, the presence of one additional vibration introduces interfering pathways and reshapes non-perturbative excitation transfer. Our results establish a scalable, hardware-efficient route to simulate vibronic processes with engineered environments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Further rotation reversal studies in C-Mod L-mode plasmas

Studies of core toroidal rotation reversal phenomenology in C-Mod deuterium L-mode plasmas have been expanded to include details of the dependences on plasma current and toroidal magnetic field. Rotation reversal occurs at a critical density, and universal scaling indicates that the product of ncritq95R ∼ BT/2, with ncrit in 1020/m3, R in m, and BT in T. Measurements in H and He plasmas exhibit similar behavior, including a connection with the linear Ohmic confinement/saturated Ohmic confinement transition and the cutoff for non-diffusive heat transport. Electron density and ion cyclotron range of frequencies power modulation experiments suggest that the collisionality ν* is a unifying parameter. Strong impurity puffing causes the critical density to increase, indicating that the situation is more complicated than only collisionality, perhaps involving the details of the effects of dilution on ion temperature gradient mode stability.

Physics↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗
DOE OSTI2023

$μ^+$ Knight shift in UTe 2 : Evidence for relocalization in a Kondo lattice

The local magnetic susceptibility of the spin-triplet superconductor UTe 2 has been investigated by positive muon (μ + ) Knight shift measurements in the normal state. Three distinct μ + Knight shift components are observed for a magnetic field applied parallel to the c axis. Two of these exhibit a breakdown in the linear relationship with the bulk magnetic susceptibility (χ) below a temperature T*~30K, which points to a gradual emergence of a correlated Kondo liquid. Below T r ~12K linearity is gradually restored, indicating relocalization of the Kondo liquid quasiparticles. The third Knight shift component is two orders of magnitude larger, and despite the c-axis alignment of the external field, scales with the a axis χ above T r ~12K. We conjecture that this component is associated with magnetic clusters and the change in the temperature dependence of all three Knight shift components below T r is associated with a change in magnetic correlations. Finally, our findings indicate that prior to the onset of superconductivity the development of the itinerant heavy-electron fluid is halted by a gradual development of local U 5⁢ƒ-moment fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

RM3 Wave Tank Validation Model

An approximately 1/75th scale point absorber wave energy absorber was built to validate the testing systems of a 16k gallon single paddle wave tank. The model was build based on the RM3 design and incorporated a linear position sensor, a force transducer, and wetness detection sensors. The data set also includes motion tracking data of the device's two bodies acquired from 4x Qualisys cameras. The tank wave spectrum is measured by 4 ultrasonic water height sensors.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

A Constraint on Primordial B-modes from the First Flight of the Spider Balloon-borne Telescope

We present the first linear polarization measurements from the 2015 long-duration balloon flight of Spider, which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in the Spider bands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis of Spider and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and r parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, finding r < 0.11 and r < 0.19, respectively. Roughly half the uncertainty in r derives from noise associated with the template subtraction. New data at 280 GHz from Spider's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

High-energy micrometre-scale pixel direct conversion X-ray detector

The objective of this work was to fabricate and characterize a new X-ray imaging detector with micrometre-scale pixel dimensions (7.8 µm) and high detection efficiency for hard X-ray energies above 20 keV. A key technology component consists of a monolithic hybrid detector built by direct deposition of an amorphous selenium film on a custom designed CMOS readout integrated circuit. Characterization was carried out at the synchrotron beamline 1-BM-B at the Advanced Photon Source of Argonne National Laboratory. The direct conversion detector demonstrated micrometre-scale spatial resolution with a 63 keV modulation transfer function of 10% at Nyquist frequency. In addition, spatial resolving power down to 8 µm was determined by imaging a transmission bar target at 21 keV. X-ray signal linearity, responsivity and lag were also characterized in the same energy range. Finally, phase contrast edge enhancement was observed in a phase object placed in the beam path. This amorphous selenium/CMOS detector technology can address gaps in commercially available X-ray detectors which limit their usefulness for existing synchrotron applications at energies greater than 50 keV; for example, phase contrast tomography and high-resolution imaging of nanoscale lattice distortions in bulk crystalline materials using Bragg coherent diffraction imaging. Here, the technology will also facilitate the creation of novel synchrotron imaging applications for X-ray energies at or above 20 keV.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multi-service battery energy storage system optimization and control

Battery energy storage systems (BESS) have become fundamental part of modern power systems due to their capability to provide multiple grid services. As the renewable penetration increases, BESS procurement is also expected to increase where it is envisioned to play a systematic and strategical role in power systems planning and operation. Hence, in this paper we present a multiple grid service procurement and operation for BESS - ranging from energy arbitrage, reserve/regulation services, power factor correction, and demand management. The proposed framework considers an optimal multi-temporal dimension, designed to be operable for both planning and real-time operation. Moreover, non-linearity inherent to BESS services and uncertainty associated to market forecasts variables are addressed using techniques such as polyhedral norms and robust optimization approaches. Here, the developed model is tested using a utility-scaled BESS and the obtained results show the effectiveness of the systematic BESS multi-service planning and operation approach.

25 ENERGY STORAGE↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A MeerKAT-meets-LOFAR study of MS 1455.0 + 2232: a 590 kiloparsec ‘mini’-halo in a sloshing cool-core cluster

ABSTRACT Radio mini-haloes are poorly understood, moderately extended diffuse radio sources that trace the presence of magnetic fields and relativistic electrons on scales of hundreds of kiloparsecs, predominantly in relaxed clusters. With relatively few confirmed detections to-date, many questions remain unanswered. This paper presents new radio observations of the galaxy cluster MS 1455.0+2232 performed with MeerKAT (covering the frequency range 872−1712 MHz) and LOFAR (covering 120−168 MHz), the first results from a homogeneously selected mini-halo census. We find that this mini-halo extends for ∼590 kpc at 1283 MHz, significantly larger than previously believed, and has a flatter spectral index (α = −0.97 ± 0.05) than typically expected. Our X-ray analysis clearly reveals a large-scale (254 kpc) sloshing spiral in the intracluster medium. We perform a point-to-point analysis, finding a tight single correlation between radio and X-ray surface brightness with a super-linear slope of $b_{\rm 1283\, MHz} = 1.16^{+0.06}_{-0.07}$ and $b_{\rm 145\, MHz} = 1.15^{+0.09}_{-0.08}$; this indicates a strong link between the thermal and non-thermal components of the intracluster medium. Conversely, in the spectral index/X-ray surface brightness plane, we find that regions inside and outside the sloshing spiral follow different correlations. We find compelling evidence for multiple sub-components in this mini-halo for the first time. While both the turbulent (re-)acceleration and hadronic scenarios are able to explain some observed properties of the mini-halo in MS 1455.0+2232, neither scenario is able to account for all the evidence presented by our analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optical magnons with dominant bond-directional exchange interactions in the honeycomb lattice iridate α — Li 2 IrO 3

We have used resonant inelastic x-ray scattering to reveal optical magnons in a honeycomb lattice iridate α - Li 2 IrO 3 . The spectrum in the energy region 20–25 meV exhibits momentum dependence, of which energy is highest at the location of the magnetic Bragg peak, ( h , k ) = ( ± 0.32 , 0 ) , and lowered toward ( 0 , 0 ) and ( ± 1 , 0 ) . We compare our data with a linear spin-wave theory based on a generic nearest-neighbor spin model. We find that a dominant bond-directional Kitaev interaction of order 20 meV is required to explain the energy scale observed in our study. The observed excitations are understood as stemming from optical magnon modes whose intensity is modulated by a structure factor, resulting in the apparent momentum dependence. We also observed diffuse magnetic scattering arising from the short-range magnetic correlation well above T N . In contrast to Na 2 IrO 3 , this diffuse scattering lacks the C 3 rotational symmetry of the honeycomb lattice, suggesting that the bond anisotropy is far from negligible in α - Li 2 IrO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Semi-Dirac Fermions in a Topological Metal

Topological semimetals with massless Dirac and Weyl fermions represent the forefront of quantum materials research. In two dimensions, a peculiar class of fermions that are massless in one direction and massive in the perpendicular direction was predicted 16 years ago. These highly exotic quasiparticles—the semi-Dirac fermions—ignited intense theoretical and experimental interest but remain undetected. Using magneto-optical spectroscopy, we demonstrate the defining feature of semi-Dirac fermions— B 2 / 3 scaling of Landau levels—in a prototypical nodal-line metal ZrSiS. In topological metals, including ZrSiS, nodal lines extend the band degeneracies from isolated points to lines, loops, or even chains in the momentum space. With calculations and theoretical modeling, we pinpoint the observed semi-Dirac spectrum to the crossing points of nodal lines in ZrSiS. Crossing nodal lines exhibit a continuum absorption spectrum but with singularities that scale as B 2 / 3 at the crossing. Our work sheds light on the hidden quasiparticles emerging from the intricate topology of crossing nodal lines and highlights the potential to explore quantum geometry with linear optical responses. Published by the American Physical Society 2024

Shao, Yinming (ORCID:0000000228910028)↗
DOE OSTI2020

s-blas

This package includes the implementation for four sparse linear algebra kernels: Sparse-Matrix-Vector-Multiplication (SpMV), Sparse-Triangular-Solve (SpTRSV), Sparse-Matrix-Transposition (SpTrans) and Sparse-Matrix-Matrix-Multiplication (SpMM) for Single-node Multi-GPU (scale-up) platforms such as NVIDIA DGX-1 and DGX-2.

Li, Ang↗
DOE OSTI2021

Topographical shifts in fine root lifespan in a mixed, mesic temperate forest

Root lifespan, often is estimated in landscape- and ecosystem-level carbon models using linear approximations. In water manipulation experiments, fine root lifespan can vary with soil water content. Soil water content is generally structured by complex topography, which is largely unaccounted for in landscape- and ecosystem-scale carbon models. Topography governs the range of soil water content experienced by roots which may impact their lifespan. We hypothesized that root lifespan varied nonlinearly across a temperate, mesic, forested catchment due to differences in soil water content associated with topographic position. We expected regions of the landscape that were too wet or too dry would have soils that were not optimal for roots and thus result in shorter root lifespans. Specifically, we hypothesized that root lifespan would be longest in areas that consistently had soil water content in the middle of the soil water content spectrum, while in soils at either very low or very high soil water content, root lifespan would be relatively short. We tested this hypothesis by collecting and analyzing two years of minirhizotron and soil moisture data in plots widely distributed in the Shale Hills catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Pennsylvania. We found that fine root lifespans were longer in traditionally wetter topographic regions, but detected no short term (biweekly) effect of soil moisture on root lifespan. Additionally, depth in soil, soil series, slope face orientation, and season of birth strongly affected root lifespans across the catchment. In contrast, lifespan was unaffected by root diameter or mycorrhizal association. Failure to account for these variables could result in erroneous estimates of fine root lifespan and, consequentially, carbon flux in temperate forested regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Methods for Computing Physically Realistic Estimates of Electric Water Heater Demand Response Resource Suitable for Bulk Power System Planning Models

Demand response is commonly called on to reduce load during system peak times or to respond to contingency events. In future power systems with higher shares of wind and solar generation (which we describe together as variable generation [VG]), demand response could have more opportunities to provide energy shifting or operating reserve services. This report evaluates the ability of residential electric water heaters, both electric resistance water heaters (ERWHs) and heat pump water heaters (HPWHs), to provide such services starting from detailed whole-building energy models that realistically represent New England single family home stock. We use a parsimonious surrogate model to represent operational flexibility in a form suitable for linear and mixed integer programming. This enables relatively fast determination of aggregate contingency reserve resource, price-taking energy shifting outcomes, and in some cases the determination of aggregate models at the megawatt (MW) scale that can be directly included in large-scale grid models. After selecting modeling methods and parameters through various computational experiments, we find interquartile ranges of contingency reserve resource in ISO-NE for about 603,400 ERWHs of 45 MW - 69 MW for Claim10 (50 minute responses provided with 10 minutes of advanced notification) and 65 MW - 102 MW for Claim30 (30 minute responses provided with 30 minutes of advanced notification), and for about 619,000 HPWHs of 48 MW - 88 MW for Claim10 and 52 MW - 90 MW for Claim30. The overall reserve resource is up to 32% of total load for ERWHs providing Claim10 service, 47% for ERWHs providing Claim30 service, 93% for HPWHs providing Claim10 service, and 97% for HPWHs providing Claim30 service. More work is required to determine if HPWHs are inherently more suitable than ERWHs for providing contingency reserve or if these results reflect idiosyncrasies of the single family home stock model used in this study. The value of this contingency resource in a Near-term VG model of ISO-NE is $\$ 0.40$ to $\$1.20$ per water heater-year, and significantly larger, $\$ 3.80$ to $\$ 5.30$ per water heater-year in a Mid-term VG model of ISONE. Aggregating surrogate models to the MW-scale for energy shifting service is more challenging than for contingency service and we only present such results for ERWHs, because we were unable to determine satisfactory ways to deal with HPWHs' time-varying and path dependent operational characteristics. Individual surrogate models suitable for evaluating the energy shifting resource from both ERWHs and HPWHs are created, however, and dispatched against day-ahead prices from the Near-Term VG and Mid-Term VG models of ISO-NE. The individual surrogate models are able to access and potentially shift all 640 GWh of HPWH load and 1,547 GWh of ERWH load we modeled in two different single family home stock models. In contrast, the most effective model of aggregate ERWH shifting resource we created only captured 34.7% of the total ERWH load. Energy shifting affected by price-taking dispatch against modeled day-ahead energy prices produces per water heater year profits of $\$19.44$ - $\$22.93$ for individual HPWHs, $\$39.11$ - $\$40.54$ for individual ERWHs, and up to $\$4.00$ - $\$4.24$ for aggregated ERWHs, with the variations mainly due to grid conditions (more or less VG). When the supply-side response to these changes is accounted for, the per water heater year production cost savings for ISO-NE are $\$7.50$ to $\$17.70$ for the most effective set of endogenously dispatched aggregate ERWHs, $\$15.60$ to $\$15.70$ for individual ERWHs dispatched against the DA prices, and $\$10.70$ to $\$11.20$ for individual HPWHs dispatched against DA prices. Those ranges primarily represent the difference between Near-Term VG and Mid-Term VG grid conditions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗