Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Perturbative methods”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-local slicing approaches for NNLO QCD in MCFM

We present the implementation of several processes at Next-to-Next-to Leading Order (NNLO) accuracy in QCD in the parton-level Monte Carlo program MCFM. The processes treated are pp → H, W ± , Z, W ± H, ZH, W ± γ, Zγ and γγ and, for the first time in the code, W + W - , W ± Z and ZZ. Decays of the unstable bosons are fully included, resulting in a flexible fully differential Monte Carlo code. The NNLO corrections have been calculated using two non-local slicing approaches, isolating the doubly unresolved region by cutting on the zero-jettiness, $\mathscr{T}_0$, or on q T , the transverse momentum of the colour singlet final-state particles. We find that for most, but not all processes the q T slicing method leads to smaller power corrections for equal computational burden.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

3D Electron Diffraction Structure Determination of Terrylene, a Promising Candidate for Intermolecular Singlet Fission

Herein we demonstrate the prowess of the 3D electron diffraction approach by unveiling the structure of terrylene, the third member in the series of peri-condensed naphthalene analogues, which has eluded structure determination for 65 years. The structure was determined by direct methods using electron diffraction data and corroborated by dispersion-inclusive density functional theory optimizations. Terrylene crystalizes in the monoclinic space group P2 1 /a, arranging in a sandwich-herringbone packing motif, similar to analogous compounds. Having solved the crystal structure, we use many-body perturbation theory to evaluate the excited-state properties of terrylene in the solid-state. We find that terrylene is a promising candidate for intermolecular singlet fission, comparable to tetracene and rubrene.

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

Superconducting magnesium diboride films for levitation of laser targets

Inertial confinement fusion targets would benefit from being levitated inside hohlraums to avoid capsule support-related implosion perturbations. Levitation inside a magnetic trap requires coating the capsule with a thin film that is superconducting at 20 K. Such non-epitaxial film growth on non-planar substrates is challenging. Here, we study Mg vapor annealing and solid-phase reactive inter-diffusion methods to form superconducting magnesium diboride (MgB 2 ) films on different planar and spherical carbon substrates, evaluating glassy carbon, polycrystalline diamond made by chemical vapor deposition, and carbon deposited by magnetron sputtering of graphite targets. Thin films of B and Mg are produced by magnetron sputtering onto stationary planar or rolling spherical substrates and annealed at either 600 or 850 °C in Mg vapor. The films are characterized by a combination of high-energy ion scattering, electron microscopy, and magnetometry. Results show that the critical superconducting temperature of resultant films depends on film microstructure and oxygen impurity content. The formation of MgB 2 films is also strongly substrate dependent, even in this case of non-epitaxial growth. Important factors to consider are oxygen outgassing of the substrate during thermal processing, substrate surface roughness, and the matching of thermal expansion coefficients of different layers in the multilayer structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A high-throughput framework for lattice dynamics

We develop an automated high-throughput workflow for calculating lattice dynamical properties from first principles including those dictated by anharmonicity. The pipeline automatically computes interatomic force constants (IFCs) up to 4th order from perturbed training supercells, and uses the IFCs to calculate lattice thermal conductivity, coefficient of thermal expansion, and vibrational free energy and entropy. It performs phonon renormalization for dynamically unstable compounds to obtain real effective phonon spectra at finite temperatures and calculates the associated free energy corrections. The methods and parameters are chosen to balance computational efficiency and result accuracy, assessed through convergence testing and comparisons with experimental measurements. Deployment of this workflow at a large scale would facilitate materials discovery efforts toward functionalities including thermoelectrics, contact materials, ferroelectrics, aerospace components, as well as general phase diagram construction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Nonperturbative Guiding Center Model for Magnetized Plasmas

Perturbative guiding center theory adequately describes the slow drift motion of charged particles in the strongly magnetized regime characteristic of thermal particle populations in various magnetic fusion devices. However, it breaks down for particles with large-enough energy. Here, we report on a data-driven method for learning a nonperturbative guiding center model from full-orbit particle simulation data. We show the data-driven model significantly outperforms traditional asymptotic theory in magnetization regimes appropriate for fusion-born α particles in stellarators, thus opening the door to nonperturbative guiding center calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fabrication of silicon-tipped fiber-optic temperature sensors using aerogel-assisted glass soldering with precise laser heating

We demonstrate the fabrication of fiber-optic Fabry–Perot interferometer (FPI) temperature sensors by bonding a small silicon diaphragm to the tip of an optical fiber using low melting point glass powders heated by a 980 nm laser on an aerogel substrate. The heating laser is delivered to the silicon FPI using an optical fiber, while the silicon temperature is being monitored using a 1550 nm white-light system, providing localized heating with precise temperature control. The use of an aerogel substrate greatly improves the heating efficiency by reducing the thermal loss of the bonding parts to the ambient environment. A desirable temperature for bonding can be achieved with relatively small heating laser power. The bonding process is carried out in an open space at room temperature for convenient optical alignment. The precise temperature control ensures minimum perturbation to the optical alignment and no induced thermal damage to the optical parts during the bonding process. For demonstration, we fabricated a low-finesse and high-finesse silicon FPI sensor and characterized their measurement resolution and temperature capability. The results show that the fabrication method has a good potential for high-precision fabrication of fiber-optic sensors.

Sheng, Qiwen↗
DOE OSTI2024

An extensible perturbed parameter ensemble for the Community Atmosphere Model version 6

This paper documents the methodology and preliminary results from a perturbed parameter ensemble (PPE) technique, where multiple parameters are varied simultaneously and the parameter values are determined with Latin hypercube sampling. This is done with the Community Atmosphere Model version 6 (CAM6), the atmospheric component of the Community Earth System Model version 2 (CESM2). We apply the PPE method to CESM2–CAM6 to understand climate sensitivity to atmospheric physics parameters. The initial simulations vary 45 parameters in the microphysics, convection, turbulence and aerosol schemes with 263 ensemble members. These atmospheric parameters are typically the most uncertain in many climate models. Control simulations and targeted simulations to understand climate forcing due to aerosols and fast climate feedbacks are analyzed. The use of various emulators is explored in the multi-dimensional space mapping input parameters to output metrics. Parameter impacts on various model outputs, such as radiation, cloud and aerosol properties, are evaluated. Machine learning is also used to probe optimal parameter values against observations. Our findings show that PPE is a valuable tool for climate uncertainty analysis. Furthermore, by varying many parameters simultaneously, we find that many different combinations of parameter values can produce results consistent with observations, and thus careful analysis of tuning is important. The CESM2–CAM6 PPE is publicly available and extensible to other configurations to address questions of other model processes in the atmosphere and other model components (e.g., coupling to the land surface).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Integrating field observations and process-based modeling to predict watershed water quality under environmental perturbations

Watersheds play a critical role in supplying water resources needed for human use and ecosystem health. Understanding and predicting how, when, and where changes in the quantity and quality of water resources occur under different environmental stresses including extreme events is crucial for sustainable management of water resources under a changing environment. However, few studies have attempted to quantify or identify the factors and process interactions controlling the impact of extreme events across water-shed systems. Only few large-scale studies include coordinated monitoring and modeling efforts, which limits our ability to assess the large-scale impact of extreme events on water supply and quality. Methods are lacking to propagate uncertainty in process understanding through an integrated hydro-biogeochemical model framework and evaluate its importance, thus failing to take full advantage of the information potentially available through transformative advances in characterization technologies from high-resolution mass spectrometry to airborne and satellite-based remote sensing. There are consequent risks to our nations water security and to human and ecosystem health that may become exacerbated with the increasing frequency of extreme events that is projected for the coming decades. This paper reviews the current status of watershed science for both water quantity and quality and identifies critical gaps in our current knowledge and modeling capability in addressing the emergent needs in predicting watershed hydrologic and biogeochemical responses (i.e., water quantity and quality) under natural and anthropogenic perturbations. We highlight the need to (1) understand how environmental perturbations including extreme events like floods and droughts propagate through watershed systems and assess their short- and long-term impacts on watershed biogeochemistry, water quality and their recovery pathways; (2) develop and improve a watershed water quality model that reflects the state of scientific understanding gained from observations; and (3) construct a data-model fusion system for watershed characterization, process identification, and mechanistic model parameterization. A large base of modeling, monitoring and data capabilities have been built by various federal government agencies given the relevance of water to their critical missions. An emerging need is to build an integrated national capability for watershed water availability and quality that can address water-related missions across multiple federal agencies.

13 HYDRO ENERGY↗
DOE OSTI2025

A semi-analytic estimate for the effective sound speed counterterm in the EFTofLSS

The Effective Field Theory of Large Scale Structure (EFTofLSS) has found tremendous success as a perturbative framework for the evolution of large scale structure, and it is now routinely used to compare theoretical predictions against cosmological observations. The model for the total matter field includes one nuisance parameter at 1-loop order, the effective sound speed, which can be extracted by matching the EFT to full N-body simulations. In this work we first leverage the Layzer-Irvine cosmic energy equation to show that the equation of state can be exactly computed with knowledge of the fully nonlinear power spectrum. When augmented with separate universe methods, we show one can estimate the effective sound speed. This estimate is in good agreement with simulation results, with errors at the few tens of percent level. Here, we apply our method to investigate the cosmology dependence of the effective sound speed and to shed light on what cosmic structures shape its value.

Cosmological perturbation theory in GR and beyond↗
DOE OSTI2020

A novel CMB component separation method: hierarchical generalized morphological component analysis

ABSTRACT We present a novel technique for cosmic microwave background (CMB) foreground subtraction based on the framework of blind source separation. Inspired by previous work incorporating local variation to generalized morphological component analysis (GMCA), we introduce hierarchical GMCA (HGMCA), a Bayesian hierarchical graphical model for source separation. We test our method on Nside = 256 simulated sky maps that include dust, synchrotron, free–free, and anomalous microwave emission, and show that HGMCA reduces foreground contamination by $25{{\ \rm per\ cent}}$ over GMCA in both the regions included and excluded by the Planck UT78 mask, decreases the error in the measurement of the CMB temperature power spectrum to the 0.02–0.03 per cent level at ℓ > 200 (and $\lt 0.26{{\ \rm per\ cent}}$ for all ℓ), and reduces correlation to all the foregrounds. We find equivalent or improved performance when compared to state-of-the-art internal linear combination type algorithms on these simulations, suggesting that HGMCA may be a competitive alternative to foreground separation techniques previously applied to observed CMB data. Additionally, we show that our performance does not suffer when we perturb model parameters or alter the CMB realization, which suggests that our algorithm generalizes well beyond our simplified simulations. Our results open a new avenue for constructing CMB maps through Bayesian hierarchical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Highest fusion performance without harmful edge energy bursts in tokamak

The path of tokamak fusion and International thermonuclear experimental reactor (ITER) is maintaining high-performance plasma to produce sufficient fusion power. This effort is hindered by the transient energy burst arising from the instabilities at the boundary of plasmas. Conventional 3D magnetic perturbations used to suppress these instabilities often degrade fusion performance and increase the risk of other instabilities. This study presents an innovative 3D field optimization approach that leverages machine learning and real-time adaptability to overcome these challenges. Implemented in the DIII-D and KSTAR tokamaks, this method has consistently achieved reactor-relevant core confinement and the highest fusion performance without triggering damaging bursts. This is enabled by advances in the physics understanding of self-organized transport in the plasma edge and machine learning techniques to optimize the 3D field spectrum. The success of automated, real-time adaptive control of such complex systems paves the way for maximizing fusion efficiency in ITER and beyond while minimizing damage to device components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Transformer-based operator learning framework for self-energy in strongly correlated systems

We introduce Σ-Attention, a transformer-based operator-learning framework for approximating the self-energy operator of strongly correlated electronic systems. By creating a batched dataset that combines results from three complementary approaches, i.e., many-body perturbation theory, strong-coupling expansion, and exact diagonalization, each effective in specific parameter regimes, Σ-Attention is applied to learn an accurate approximation for the self-energy operator that is valid across a wide range of parameter regimes. This hybrid strategy leverages the strengths of existing methods while relying on the transformer's ability to generalize beyond individual limitations. More importantly, the scalability of the transformer architecture allows the learned self-energy to be extended to systems with larger sizes, leading to much improved computational scaling. Using the one-dimensional Hubbard model, we demonstrate that Σ-Attention can accurately predict the Matsubara Green's function of large systems with a wide range of coupling strength. Our framework offers a promising and scalable pathway for studying strongly correlated systems with many possible generalizations.

Zhu, Yuanran↗
DOE OSTI2025

Object detection with deep learning for rare event search in the GADGET II TPC

In the pursuit of identifying rare two-particle events within the GADGET II Time Projection Chamber (TPC), this paper presents a comprehensive approach for leveraging Convolutional Neural Networks (CNNs) and various data processing methods. To address the inherent complexities of 3D TPC track reconstructions, the data is expressed in 2D projections and 1D quantities. This approach capitalizes on the diverse data modalities of the TPC, allowing for the efficient representation of the distinct features of the 3D events, with no loss in topology uniqueness. Additionally, it leverages the computational efficiency of 2D CNNs and benefits from the extensive availability of pre-trained models. Given the scarcity of real training data for the rare events of interest, simulated events are used to train the models to detect real events. To account for potential distribution shifts when predominantly depending on simulations, significant perturbations are embedded within the simulations. This produces a broad parameter space that works to account for potential physics parameter and detector response variations and uncertainties. These parameter-varied simulations are used to train sensitive 2D CNN object detectors. When combined with 1D histogram peak detection algorithms, this multi-modal detection framework is highly adept at identifying rare, two-particle events in data taken during experiment 21072 at the Facility for Rare Isotope Beams (FRIB), demonstrating a 100% recall for events of interest. Here, we present the methods and outcomes of our investigation and discuss the potential future applications of these techniques.

Convolutional neural network↗
DOE OSTI2022

Impact of Increased Latent Generations on Sensitivity Calculations with SCALE

Analyses of cross section sensitivity data from systems with fissile material allow analysts to associate an importance for each material, nuclide, reaction, and neutron energy by simulating real world criticality scenarios. Although criticality safety validation efforts can be guided by the cross-section sensitivity and uncertainty data generated for a particular system, these calculations can often be computationally expensive and sometimes cumbersome without proper guidance. The TSUNAMI suite within the SCALE code package has several methods for generating sensitivity data, including multigroup and continuous energy (CE) capabilities. The release of SCALE 6.3 has three different CE methods for generating cross section sensitivity data: (1) the Iterated Fission Probability (IFP) method with the KENO Monte Carlo transport solver, (2) the IFP method with the Shift Monte Carlo transport solver, and (3) the Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) method with the KENO Monte Carlo transport solver. Although the CLUTCH method has additional parameters for generating sensitivity data files relative to the IFP method, all three methods use latent generations, which are the generations between an event (i.e., fission) and the assessment of importance based on the asymptotic population of progeny neutrons. Increasing the number of latent generations in a calculation leads to increased discrimination of the sensitivity coefficients but at the cost of the increased uncertainty associated with those generated values. Analysts must balance the accuracy of the sensitivity calculations and its uncertainty with the associated computational cost involved in generating the values. This paper discusses the impact of adjusting the latent generation parameter for a range of sensitivity values and how these changes compare with the direct perturbation values obtained from a change of ±0.5% Δ k in both benchmark and safety application models. Two benchmarks from the International Handbook of Evaluated Criticality Safety Benchmark Experiments and the MPC-32 dual purpose canister for spent nuclear fuel are used for analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Heterotic instantons for monad and extension bundles

We consider non-perturbative superpotentials from world-sheet instantons wrapped on holomorphic genus zero curves in heterotic string theory. These superpotential contributions feature prominently in moduli stabilization and large field axion inflation, which makes their presence or absence, as well as their functional dependence on moduli, an important issue. We develop geometric methods to compute the instanton superpotentials for heterotic string theory with monad and extension bundles. Using our methods, we find a variety of examples with a non-vanishing superpotential. In view of standard vanishing theorems, we speculate that these results are likely to be attributed to the non-compactness of the instanton moduli space. We test this proposal, for the case of monad bundles, by considering gauged linear sigma models where compactness of the instanton moduli space can be explicitly checked. In all such cases, we find that the geometric results are consistent with the vanishing theorems. Surprisingly, linearly dependent Pfaffians even arise for cases with a non-compact instanton moduli space. This suggests some gauged linear sigma models with a non-compact instanton moduli space may still have a vanishing instanton superpotential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗