Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite precision”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

we report a multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Quantifying Errors in Effective Cluster Interactions of Lattice Gas Cluster Expansions

The promise of lattice gas (LG) cluster expansions (CEs) is that they can describe a given system property to any level of accuracy since the orthogonal “cluster basis functions” span the complete space of available configurations. Such an approach can be constructed to an arbitrarily large surface of a finite number of distinct adsorption sites. Unfortunately, this is only true for the case of an ideal, fixed lattice decorated with components at precise lattice points (the lattice “sites”) with no distortions or relaxations subsequently allowed. Since most systems, and surfaces specifically, do not conform to such an ideal set of constraints, errors in LG CEs must be expected or CE convergence severely hampered. Beyond this, numerical errors in the provided data can complicate the proper construction of a truly predictive and/or physically significant CE. In this work, we show here how reliance on typical statistical tools like confidence intervals cannot be expected to provide an accurate representation of the uncertainty of effective cluster interactions (ECIs) in the CE due to the nature of the target ab initio data and the nature of CEs themselves. We develop a method for estimating these errors that does not rely on statistical assumptions about the model or data. We then use these ECI errors to quantify fundamental consequences on the uncertainty of ECIs in CEs built from O/Fe(100) data whose surface and adsorbates have been allowed to relax in a typical manner and from O/Fe(100) data whose surface and adsorbates are fixed in ideal lattice positions. We also quantify the effect of using a different density functional theory exchange–correlation functional, using these ECI errors to assess the significance in any deviations. In both cases, our method is shown to have remarkable utility in the quantification of errors in the ECIs of CEs. While we stick to the lattice gas convention in this work, the method is equally applicable to the Ising convention or, in principle, any linear model of sufficient complexity.

08 HYDROGEN↗
DOE OSTI2022

Hot Electron-Driven Photocatalysis Using Sub-5 nm Gap Plasmonic Nanofinger Arrays

Semiconductor photocatalysis has received increasing attention because of its potential to address problems related to the energy crisis and environmental issues. However, conventional semiconductor photocatalysts, such as TiO2 and ZnO, can only be activated by ultraviolet light due to their wide band gap. To extend the light absorption into the visible range, the localized surface plasmon resonance (LSPR) effect of noble metal nanoparticles (NPs) has been widely used. Noble metal NPs can couple incident visible light energy to strong LSPR, and the nonradiative decay of LSPR generates nonthermal hot carriers that can be injected into adjacent semiconductor material to enhance its photocatalytic activity. Here we demonstrate that nanoimprint-defined gap plasmonic nanofinger arrays can function as visible light-driven plasmonic photocatalysts. The sub-5 nm gaps between pairs of collapsed nanofingers can support ultra-strong plasmon resonance and thus boost the population of hot carriers. The semiconductor material is exactly placed at the hot spots, providing an efficient pathway for hot carrier injection from plasmonic metal to catalytic materials. This nanostructure thus exhibits high plasmon-enhanced photocatalytic activity under visible light. The hot carrier injection mechanism of this platform was systematically investigated. The plasmonic enhancement factor was calculated using the finite-difference time-domain (FDTD) method and was consistent with the measured improvement of the photocatalytic activity. This platform, benefiting from the precise controllable geometry, provides a deeper understanding of the mechanism of plasmonic photocatalysis.

42 ENGINEERING↗
DOE OSTI2021

Multi-species collisions for delta-f gyrokinetic simulations: Implementation and verification with GENE

A multi-species linearized collision operator based on the model developed by Sugama et al. has been implemented in the nonlinear gyrokinetic code, GENE. Such a model conserves particles, momentum, and energy to machine precision, and is shown to have negative definite free energy dissipation characteristics, satisfying Boltzmann’s H-theorem, including for realistic mass ratio. Finite Larmor Radius (FLR) effects have also been implemented into the local version of the code. For the global version of the code, the collision operator has been developed to allow for block-structured velocity space grids, allowing for computationally tractable collisional global simulations. The validity of the collision operator has been demonstrated by relaxation and conservation tests, as well as appropriate benchmarks. The newly implemented operator shall be used in future simulations to study magnetically confined fusion plasma turbulence and transport in more extreme regions with higher collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Fast GPU 3D diffeomorphic image registration

3D image registration is one of the most fundamental and computationally expensive operations in medical image analysis. Here, we present a mixed-precision, Gauss–Newton–Krylov solver for diffeomorphic registration of two images. Our work extends the publicly available CLAIRE library to GPU architectures. Despite the importance of image registration, only a few implementations of large deformation diffeomorphic registration packages support GPUs. Our contributions are new algorithms to significantly reduce the run time of the two main computational kernels in CLAIRE: calculation of derivatives and scattered-data interpolation. Additionally, we deploy (i) highly-optimized, mixed-precision GPU-kernels for the evaluation of scattered-data interpolation, (ii) replace Fast-Fourier-Transform (FFT)-based first-order derivatives with optimized 8th-order finite differences, and (iii) compare with state-of-the-art CPU and GPU implementations. As a highlight, we demonstrate that we can register clinical images in less than 6 s on a single NVIDIA Tesla V100. This amounts to over 20 speed-up over the current version of CLAIRE and over 30 speed-up over existing GPU implementations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Universal relation between doping content and normal-state resistance in gate voltage tuned ultrathin Bi 2 Sr 2 CaCu 2 O 8 + x flakes

We report gate voltage tunable ultra-thin high-T c cuprates supply a unique platform to investigate the electronic phase diagram and superconductor-insulator transition. One of the challenges in this field is the precise determination of the doping content in the underdoped non-superconducting region. Here we report the discovery of a universal relation between the doping content p and the normal-state resistance at a fixed temperature R(T f ), p = α + β ln[1/R(T f )], in the ultra-thin Bi 2 Sr 2 CaCu 2 O 8+x flakes. The in-depth analysis shows that the evolution of carrier scattering probability with doping content and the change of effective mass caused by superconductor-insulator transition are two key factors leading to this logarithmic relation. Based on our finding, the more precise electronic phase diagram can be established. In addition, the superconductor-insulator transition is verified to be a quantum phase transition using a finite size scaling analysis. The scaling exponent z v is found to have a close correlation with the disorder levels. The present result provides an important foundation to investigate the fascinating electronic states in the ultra-thin cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effect of time and spatial domains on monolayer 2D material interface thermal conductance measurement using ns ET-Raman

Raman-based approaches to characterizing thermal properties of nanoscale and microscale materials have proven to be powerful experimental methods. Here, this work explores the potential of energy-transport state-resolved Raman (ET-Raman) for the precise determination of interfacial thermal conductance (G) between monolayer 2D materials and their substrates. ET-Raman employs a continuous wave laser and nanosecond (ns) pulsed laser to induce both steady state and transient heating while simultaneously collecting 2D material temperature information from the Raman signal. A 3D numerical model built using the finite volume method simulates the two differentiated thermal transport processes. The simulated temperature rises are compared with experimental Raman temperature information to extract G. In this study, we investigate the temporal effects of the pulse frequency and pulse width of the ns laser on precise characterization of G. By running multiple simulations for a range of pulsed laser frequencies and pulse widths, we determine that ET-Raman is most sensitive to precise interfacial thermal conductance measurements with short pulse widths (~below 20 ns) and with pulsing frequencies that ensure sufficient cooling time between consecutive pulses (~below 500 kHz). We also conclude that substrate thickness and laser spot size play critical roles in heat distribution and subsequent sensitivity to meaningful G measurements. As both laser spot size and SiO 2 substrate thickness decrease, detection of interfacial thermal transport becomes more difficult. These results offer essential guidance in the design of future Raman-based approaches to thermal measurements using pulsed lasers.

2D materials↗
DOE OSTI2024

Broad sterile neutrinos & the reactor/gallium tension

Significant evidence exists for the apparent disappearance of electron-type neutrinos in radioactive source experiments. Yet, interpreted within the standard ‘3+1 sterile neutrino scenario’, precision short-baseline measurements of electron antineutrinos from nuclear reactors strongly disagree with these results. Recently, it has been demonstrated that allowing for a finite wavepacket size for the reactor neutrinos can ameliorate such a tension, however the smallness of the required wavepackets is a subject of intense debate. In this work, we demonstrate that a ‘broad’ sterile neutrino may relax this tension in much the same way. Such a phenomenological possibility can arise in plausible hidden sector scenarios, such as a clockwork-style sector, for which we provide a concrete microscopic model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Normalization of ZZ instanton amplitudes in minimal string theory

We use insights from string field theory to analyze and cure the divergences in the cylinder diagram in minimal string theory with both boundaries lying on a ZZ brane. We focus on theories with worldsheet matter consisting of the (2, p) minimal model plus Liouville theory, with total central charge 26, together with the usual bc-ghosts. The string field theory procedure gives a finite, purely imaginary normalization constant for non-perturbative effects in minimal string theory, or doubly non-perturbative effects in JT gravity. We find precise agreement with the prediction from the dual double-scaled one-matrix integral. We also make a few remarks about the extension of this result to the more general (p', p) minimal string.

2D gravity↗
DOE OSTI2025

Cabibbo-Kobayashi-Maskawa unitarity deficit reduction via finite nuclear size

We revisit the extraction of the |𝑉 𝑢⁢𝑑 | Cabibbo-Kobayashi-Maskawa (CKM) matrix element from the superallowed transition decay rate of 26⁢𝑚 Al → 26 Mg, focusing on finite nuclear size effects. The decay rate dependence on the 26⁢𝑚 Al charge radius is found to be four times higher than previously believed, necessitating precise determination. However, for a short-lived isotope of an odd 𝑍 element such as 26⁢𝑚 Al , radius extraction relies on challenging many-body atomic calculations. We performed the needed calculations, finding an excellent agreement with previous ones, which used a different methodology. This sets a new standard for the reliability of isotope shift factor calculations in many-electron systems. The ℱ⁡𝑡 value obtained from our analysis is lower by 2.2𝜎 than the corresponding value in the previous critical survey, resulting in an increase in |𝑉 𝑢⁢𝑑 | 2 by 0.9𝜎. Adopting |𝑉 𝑢⁢𝑑 | from this decay alone reduces the CKM unitarity deficit by one standard deviation, irrespective of the choice of |𝑉 𝑢⁢𝑠 |.

beta-decay↗
DOE OSTI2022

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots

The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.

divertor physics↗
DOE OSTI2020

Crossing symmetry, transcendentality and the Regge behaviour of 1d CFTs

We develop the technology for Polyakov-Mellin (PM) bootstrap in one- dimensional conformal field theories (CFT 1 ). By adding appropriate contact terms, we bootstrap various effective field theories in AdS 2 and analytically compute the CFT data to one loop. The computation can be extended to higher orders in perturbation theory, if we ignore mixing, for any external dimension. We develop PM bootstrap for O(N) theories and derive the necessary contact terms for such theories (which also involves a new higher gradient contact term absent for N = 1). We perform cross-checks which include considering the diagonal limit of the 2d Ising model in terms of the 1d PM blocks. As an independent check of the validity of the results obtained with PM bootstrap, we propose a suitable basis of transcendental functions, which allows to fix the four-point correlators of identical scalar primaries completely, up to a finite number of ambiguities related to the number of contact terms in the PM basis. We perform this analysis both at tree level (with and without exchanges) and at one loop. We also derive expressions for the corresponding CFT data in terms of harmonic sums. Finally, we consider the Regge limit of one-dimensional correlators and derive a precise connection between the latter and the large-twist limit of CFT data. Exploiting this result, we study the crossing equation in the three OPE limits and derive some universal constraints for the large-twist limit of CFT data in Regge-bounded theories with a finite number of exchanges.

1/N expansion↗
DOE OSTI2023

Elastic nucleon-pion scattering at m π = 200 MeV from lattice QCD

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π = 200 MeV. The s-wave scattering lengths with both total isospins I = 1/2 and I = 3/2 are inferred from the finite-volume spectrum below the inelastic threshold together with the I = 3/2 p-wave containing the Δ(1232) resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels, enabling a determination of the I = 3/2 scattering length with statistical errors below 5%, while the I = 1/2 scattering length is somewhat less precisely evaluated. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a step toward future computations down to physical light quark masses with multiple lattice spacings and volumes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Limitations to the ‘basic’ HOD model and beyond

ABSTRACT We use the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based halo occupation distribution (HOD) approach to modelling the galaxy–halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (basic) HOD model, we find that, on average, the ‘basic’ HOD model underpredicts the real-space correlation function in the TNG300 box by ∼15 per cent on scales of $1 \,\,\lt\,\, r \,\,\lt\,\, 20 \ {\rm Mpc}\, h^{-1}$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ‘secondary’ parameters with which to augment the ‘basic’ HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, β, and the product of the half-mass radius and the velocity dispersion, σ2Rhalfmass, are the three most effective measures of assembly bias that help reconcile the ‘basic’ HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch, and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy–halo connection can be constructed if we combine both mass and local environment information about the halo.

42 ENGINEERING↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Optomechanical self-organization in a mesoscopic atom array

Increasing the number of particles in a system often leads to qualitative changes in its properties, such as breaking of symmetries and the appearance of phase transitions. This renders a macroscopic system fundamentally different from its individual microscopic constituents. Lying between these extremes, mesoscopic systems exhibit microscopic fluctuations that influence behaviour on longer length scales, leading to critical phenomena and dynamics. Therefore, tracing the properties of well-controlled mesoscopic systems can help bridge the gap between an exact description of few-body microscopic systems and the emergent description of many-body systems. Here we explore the mesoscopic signatures of an optomechanical self-organization phase transition using arrays of cold atoms inside an optical cavity. By precisely engineering atom–cavity interactions, we reveal how critical behaviour depends on the atom number, identify characteristic dynamical behaviours in the self-organized regime and observe a finite optomechanical susceptibility at the critical point. These findings advance our understanding of particle-number- and time-resolved properties of phase transitions in mesoscopic systems.

Atomic and molecular interactions with photons↗