Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Exact solution of two simple non-equilibrium electron-phonon and electron-electron coupled systems: The atomic limit of the Holstein-Hubbard model and the generalized Hatsugai-Komoto model

One of the challenges in many-body physics is determining the effects of phonons on strongly correlated electrons. The difficulty arises from strong correlations at differing energy scales—for band metals, Migdal-Eliashberg theory accurately determines electron-phonon coupling effects due to the absence of vertex corrections, but strongly correlated electrons require a more complex description, and the standard Migdal-Eliashberg approach does not necessarily apply. In this work, we solve for the atomic limit Green's function of the Holstein-Hubbard model with both time-dependent electron-electron and electron-phonon couplings. Then, we examine the photoemission spectra of this model in and out of equilibrium. Next, we use similar methods to exactly solve an extended version of the Hatsugai-Komoto model and examine its behavior in and out of equilibrium. These calculations lead us to propose using the first moment of the photoemission spectra to signal non-equilibrium changes in electron-electron and electron-phonon couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dedicated beam position monitor pair for model-independent lattice characterization at NSLS-II

This paper reports recent lattice characterization results obtained at the National Synchrotron Light Source II (NSLS-II) storage ring, conducted without reliance on a lattice model. A pair of beam position monitors (BPMs) with bunch-by-bunch (B×B) resolution, were recently installed in a section of the storage ring free of magnetic fields. The new BPM pair measured the beam, or bunch’s transverse Poincar´e map precisely after the beam was excited. Linear one-turn-matrices (OTM) were then derived, and from these, the 4-dimensional coupled Twiss parameters were extracted at the locations of the BPM pair. By normalizing beam oscillation amplitudes with the Twiss parameters, the global action-variables were obtained. These action-variables facilitated the measurement of the local Twiss parameters observed by other BPMs independent on lattice model. This method is general, and particularly useful in certain scenarios such as a round beam mode in a diffraction-limited light source ring. We applied it to assess both weakly and strongly coupled lattices at the NSLS-II ring. Through analysis of the strongly coupled lattice, the quadrupole tilt errors were estimated to be less than 400 µrad. Utilizing the BPMs’ B×B resolution, for the first time we observed the variations of the linear lattice along a long bunch-train.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Dedicated beam position monitor pair for model-independent lattice characterization at NSLS-II

This paper reports recent lattice characterization results obtained at the National Synchrotron Light Source II (NSLS-II) storage ring, conducted without reliance on a lattice model. A pair of beam position monitors (BPMs) with bunch-by-bunch (B$\times$B) resolution, were recently installed in a section of the storage ring free of magnetic fields. The new BPM pair measured the beam, or bunch’s transverse Poincaré map precisely after the beam was excited. Linear one-turn-matrices (OTM) were then derived, and from these, the 4-dimensional coupled Twiss parameters were extracted at the locations of the BPM pair. By normalizing beam oscillation amplitudes with the Twiss parameters, the global action-variables were obtained. Additionally, these action-variables facilitated the measurement of the local Twiss parameters observed by other BPMs independent on lattice model. This method is general, and particularly useful in certain scenarios such as a round beam mode in a diffraction-limited light source ring. We applied it to assess both weakly and strongly coupled lattices at the NSLS-II ring. Through analysis of the strongly coupled lattice, the quadrupole tilt errors were estimated to be less than 400 μrad. Utilizing the BPMs’ B$\times$B resolution, for the first time we observed the variations of the linear lattice along a long bunch-train.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Microscopic Origins of the Nonlinear Behavior of Particle-Filled Rubber Probed with Dynamic Strain XPCS

The underlying microscopic response of filler networks in reinforced rubber to dynamic strain is not well understood due to the experimental difficulty of directly measuring filler network behavior in samples undergoing dynamic strain. This difficulty can be overcome with in situ X-ray photon correlation spectroscopy (XPCS) measurements. The contrast between the silica filler and the rubber matrix for X-ray scattering allows us to isolate the filler network behavior from the overall response of the rubber. This in situ XPCS technique probes the microscopic breakdown and reforming of the filler network structure, which are responsible for the nonlinear dependence of modulus on strain, known in the rubber science community as the Payne effect. These microscopic changes in the filler network structure have consequences for the macroscopic material performance, especially for the fuel efficiency of tire tread compounds. In this work, we elucidate the behavior with in situ dynamic strain XPCS experiments on industrially relevant, vulcanized rubbers filled (13 vol %) with novel air-milled silica of ultrahigh-surface area (UHSA) (250 m 2 /g). The addition of a silane coupling agent to rubber containing this silica causes an unexpected and counterintuitive increase in the Payne effect and decrease in energy dissipation. For this rubber, we observe a nearly two-fold enhancement of the storage modulus and virtually equivalent loss tangent compared to a rubber containing a coupling agent and conventional silica. Interpretation of our in situ XPCS results simultaneously with interpretation of traditional dynamic mechanical analysis (DMA) strain sweep experiments reveals that the debonding or yielding of bridged bound rubber layers is key to understanding the behavior of rubber formulations containing the silane coupling agent and high-surface area silica. These results demonstrate that the combination of XPCS and DMA is a powerful method for unraveling the microscale filler response to strain which dictates the dynamic mechanical properties of reinforced soft matter composites. With this combination of techniques, we have elucidated the great promise of UHSA silica when used in concert with a silane coupling agent in filled rubber. Such composites simultaneously exhibit large moduli and low hysteresis under dynamic strain.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cumulant methods for electron-phonon problems. I. Perturbative expansions

Here in this work, we investigate the ability of the cumulant expansion (CE) to capture one-particle spectral information in electron-phonon coupled systems at both zero and finite temperatures. In particular, we present a comprehensive study of the second- and fourth-order CEs for the one-dimensional Holstein model as compared with numerically exact methods. We investigate both finite sized systems as well as the approach to the thermodynamic limit, drawing distinctions, and connections between the behavior of systems in and away from the thermodynamic limit that enable a greater understanding of the ability of the CE to capture real-frequency information across the full range of wave vectors. We find that for zero electronic momentum, the spectral function is well described by the second-order CE at low and high temperatures. However, for nonzero electronic momenta, the CE is only accurate at high temperature. We analyze the fourth-order cumulant and find that while it improves the description of the short-time dynamics encoded in the one-particle Green's function, it can introduce divergences in the time domain as well as unphysical negative spectral weight in the spectral function. When well-behaved, the fourth-order CE does provide notable accurate corrections to the second-order CE. Finally, we use our results to comment on the use of the CE as a tool for calculating transport behavior in the realistic ab initio modeling of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Perspective on Many-Body Methods for Molecular Polaritonic Systems

Recent advances in strong light–matter interactions have revealed a wealth of new physical phenomena in molecules embedded in optical cavities, including modified chemical reactivity, altered excitation spectra, and novel quantum correlations. To describe these effects from first-principles, the field of ab initio quantum electrodynamics (QED) has emerged as a compelling extension of quantum chemistry that treats electronic and photonic degrees of freedom on equal footing. In this Perspective, we review the growing landscape of many-body QED methods, including Hartree–Fock, density functional theory (QEDFT), time-dependent DFT (QED-TDDFT), configuration interaction (QED-CI), complete active space (QED-CASSCF), coupled cluster (QED-CC), quantum Monte Carlo (QED-QMC), and density matrix renormalization group (QED-DMRG), highlighting recent developments and implementations. We further explore real-time methods, gradient and Hessian formalisms, and the integration of nonadiabatic nuclear dynamics. Applications range from benchmark simulations of polaritonic chemistry to quantum simulations on emerging quantum hardware. We conclude by outlining future directions for theory development and interdisciplinary efforts at the interface of quantum chemistry, condensed matter, and quantum optics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Epoxy-based infrared filter assembly and associated fabrication devices and method

Epoxy-based inline infrared (IR) filter assembly, and manufacture and use of the same. Co-axial infrared filter assemblies comprise a substantially cylindrical filter body forming a central cavity characterized by opposing holes at each end. The filter body forms an outer conductor, and SMA connectors coupled to the opposing holes at each end of the body are electrically coupled to form an inner conductor positioned along a long axis of the filter body. An infrared absorbing material (such as castable epoxy resin) fills the central cavity of the filter body. Methods for producing the co-axial infrared filter include pressing SMA connectors into the respective ends of the filter body, electrically coupling the SMA connectors, and filling the filter body with epoxy. Electronic systems for operating a dark matter detector include a feedline comprising a coaxial filter configured to advantageously block infrared noise.

Spahn, Gabriel↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Selective One-Dimensional 13 C– 13 C Spin-Diffusion Solid-State Nuclear Magnetic Resonance Methods to Probe Spatial Arrangements in Biopolymers Including Plant Cell Walls, Peptides, and Spider Silk

Two-dimensional (2D) and 3D through-space 13 C– 13 C homonuclear spin-diffusion techniques are powerful solid-state nuclear magnetic resonance (NMR) tools for extracting structural information from 13 C-enriched biomolecules, but necessarily long acquisition times restrict their applications. In this work, we explore the broad utility and underutilized power of a chemical shift-selective one-dimensional (1D) version of a 2D 13 C– 13 C spin-diffusion solid-state NMR technique. The method, which is called 1D dipolar-assisted rotational resonance (DARR) difference, is applied to a variety of biomaterials including lignocellulosic plant cell walls, microcrystalline peptide fMLF, and black widow dragline spider silk. 1D 13 C– 13 C spin-diffusion methods described here apply in select cases in which the 1D 13 C solid-state NMR spectrum displays chemical shift-resolved moieties. Furthermore, this is analogous to the selective 1D nuclear Overhauser effect spectroscopy (NOESY) experiment utilized in liquid-state NMR as a faster (1D instead of 2D) and often less ambiguous (direct sampling of the time domain data, coupled with increased signal averaging) alternative to 2D NOESY. Selective 1D 13 C– 13 C spin-diffusion methods are more time-efficient than their 2D counterparts such as proton-driven spin diffusion (PDSD) and dipolar-assisted rotational resonance. The additional time gained enables measurements of 13 C– 13 C spin-diffusion buildup curves and extraction of spin-diffusion time constants TSD, yielding detailed structural information. Specifically, selective 1D DARR difference buildup curves applied to 13 C-enriched hybrid poplar woody stems confirm strong spatial interaction between lignin and acetylated xylan polymers within poplar plant secondary cell walls, and an interpolymer distance of ~0.45–0.5 nm was estimated. Additionally, Tyr/Gly long-range correlations were observed on isotopically enriched black widow spider dragline silks.

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

Serpent - Bison - THM Preliminary Multiphysics Modeling of a Nuclear Thermal Propulsion System Fuel Assembly

This work demonstrates the Monte Carlo neutronic and Thermo-Hydraulic coupling scheme using the Serpent code and MOOSE application Bison and Thermo Hydraulic Module. The coupling scheme is then applied to the reference BWX Technologies Nuclear Ther- mal Propulsion system at he fuel assembly level where it’s used to perform an analysis of the isothermal material coefficients and potential material reactivity worth. A method is developed to isolate which feedback effects should be considered for proceeding with reduced order deterministic neutronic modeling where branch off analysis must be con- ducted. The convergence behavior of the coupling scheme is demonstrated where it fol- lows the standard Picard iteration approach. Verification studies for the method of deduc- ing relevant feedback effects is also demonstrated.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Localization dynamics in a centrally coupled system

In systems in which interactions couple a central degree of freedom and a bath, one would expect signatures of the bath's phase to be reflected in the dynamics of the central degree of freedom. This has been recently explored in connection with many-body localized baths coupled with a central qubit or a single-cavity mode - systems with growing experimental relevance in various platforms. Such models also have an interesting connection with Floquet many-body localization via quantizing the external drive, although this has been relatively unexplored. Here we adapt the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) method, a well-known tree tensor network algorithm, to numerically simulate the dynamics of a central degree of freedom, represented by a d-level system (qudit), coupled to a disordered interacting one-dimensional spin bath. ML-MCTDH allows us to reach ≈10 2 lattice sites, a far larger system size than what is feasible with exact diagonalization or kernel polynomial methods. From the intermediate time dynamics, we find a well-defined thermodynamic limit for the qudit dynamics upon appropriate rescaling of the system-bath coupling. The spin system shows similar scaling collapse in the Edward-Anderson spin-glass order parameter or entanglement entropy at relatively short times. At longer timescales, we see slow growth of the entanglement, which may arise from dephasing mechanisms in the localized system or long-range interactions mediated by the central degree of freedom. Similar signs of localization are shown to appear as well with unscaled system-bath coupling.

1-dimensional spin chains↗
DOE OSTI2021

Probing Clustering Dynamics between Silicon and PAA or LiPAA Slurries under Processing Conditions

This work explores the complex interplay between slurry aggregation, agglomeration, and conformation (i.e., shape) of poly(acrylic acid) (PAA)- and lithiated PAA-based silicon slurries as a function of the shear rate and the resulting slurry homogeneity. These values were measured by small-angle neutron scattering (SANS) and rheology-coupled ultra-SANS at conditions relevant to battery electrode casting. Different binder solution preparation methods, either a ball milling (BM) process or a planetary centrifugal mixing (PCM) process, dramatically modify the resulting polymer dynamics and organization around a silicon material. This is due to the different energy profiles of mixing where the more violent and higher energy PCM causes extensive breakdown and reformation of the binder, which is now likely in a branched conformation, while the lower energy BM results in simply lower-molecular weight linear polymers. In conclusion, the breakdown and reorganization of the polymer structure affect silicon slurry homogeneity, which affects subsequent electrode architecture.

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

Neutrinoless double β decay with light sterile neutrinos: the contact terms

We study neutrinoless double-beta decay in extensions of the Standard Model that include n right-handed neutrino singlets, with masses m s below the GeV scale. Generalizing recently developed matching methods, we determine the m s dependence of the short-range nn → pp couplings that appear to leading order in the chiral effective field theory description of neutrinoless double beta decay. We focus on two scenarios, corresponding to the minimal νSM and left-right symmetric models. We illustrate the impact of our new results in the case of the νSM, showing a significant impact on the neutrinoless double-beta decay half-life when m s is in the 200–800 MeV range.

effective field theories↗
DOE OSTI2024

Long duration battery sizing, siting, and operation under wildfire risk using progressive hedging

Battery sizing and siting problems are computationally challenging due to the need to make long-term planning decisions that are cognizant of short-term operational decisions. This paper considers sizing, siting, and operating batteries in a power grid to maximize their benefits, including price arbitrage and load shed mitigation, during both normal operations and periods with high wildfire ignition risk. Here we formulate a multi-scenario optimization problem for long duration battery storage while considering the possibility of load shedding during Public Safety Power Shutoff (PSPS) events that de-energize lines to mitigate severe wildfire ignition risk. To enable a computationally scalable solution of this problem with many scenarios of wildfire risk and power injection variability, we develop a customized temporal decomposition method based on a progressive hedging framework. Extending traditional progressive hedging techniques, we consider coupling in both placement variables across all scenarios and state-of-charge variables at temporal boundaries. This enforces consistency across scenarios while enabling parallel computations despite both spatial and temporal coupling. The proposed decomposition facilitates efficient and scalable modeling of a full year of hourly operational decisions to inform the sizing and siting of batteries. With this decomposition, we model a year of hourly operational decisions to inform optimal battery placement for a 240-bus WECC model in under 70 min of wall-clock time.

25 ENERGY STORAGE↗
DOE OSTI2024

A FFT-based mesoscale continuum dislocation mechanics with defect energy: Applications to composites and polycrystals

A crystal plasticity elastoviscoplastic FFT (fast Fourier transform) formulation with a mesoscale continuum field dislocation mechanics model is presented, which incorporates a defect energy density that depends on GND densities and an associated material length scale. This allows to thermodynamically derive internal length scale dependent intra-crystalline backstress and Peach–Koehler force acting on GND densities. The model considers GND density evolution through a filtered numerical spectral approach, which is coupled with stress equilibrium through the elastoviscoplastic FFT algorithm. The discrete Fourier transform (DFT) method together with finite difference (FD) schemes is applied to solve both the backstress tensor and the Fourier–Green operator. Numerical results are first reported for two-phase laminate composites with plastic single crystal channels and elastic precipitates for shear loadings. Channel size effects are simulated and analyzed on the overall and local hardening behaviors during monotonous loadings. In addition, the evolutions of GND densities and the role of their associated backstress on size effects are examined during reversible shear loading. In a second part, the role of the defect energy internal length scale on polycrystal’s hardening during tension–compression is discussed. The results are compared to those obtained using FFT-based continuum field dislocation mechanics without defect energy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Local modeling for FRF estimation with noisy input measurements

The frequency response function (FRF) is an essential means by which dynamic systems are qualified. In recent years, local modeling approaches have been extensively researched and shown to significantly outperform traditional FRF estimators. However, the standard local modeling approach assumes a perfectly-known system input, which results in biased FRF estimates in the presence of input noise. This paper derives a simple adjustment that can be used to improve FRF estimation for systems subjected to random excitation with noisy input data. This improvement can be implemented with little modification to standard local modeling algorithms and with little additional computational burden. The adjustment is coupled with a model selection procedure to avoid underfitting and overfitting. In conclusion, the methods presented in this paper are validated on a simulation, and they are shown to reduce bias due to input noise.

47 OTHER INSTRUMENTATION↗