Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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

Soft integrals and soft anomalous dimensions at $\mathrm{N^3LO}$ and beyond

We calculate soft phase-space and loop master integrals for the computation of color-singlet cross sections through N 3 LO in perturbative QCD. Our results are functions of homogeneous transcendental weight and include the first nine terms in the expansion in the dimensional regulator ϵ. We discuss the application of our results to the computation of deeply-inelastic scattering and e + e - annihilation processes. We use these results to compute the perturbative coefficient functions for the Drell-Yan and gluon-fusion Higgs boson production cross sections to higher orders in ϵ through N 3 LO in QCD in the limit where only soft partons are produced on top of the colorless final state. Furthermore, we extract the anomalous dimension of the inclusive threshold soft function and of the N-Jettiness beam and jet functions to N 4 LO in perturbative QCD.

higher-order perturbative calculations↗
DOE OSTI2022

A deep learning approach for semantic segmentation of unbalanced data in electron tomography of catalytic materials

In computed TEM tomography, image segmentation represents one of the most basic tasks with implications not only for 3D volume visualization, but more importantly for quantitative 3D analysis. In case of large and complex 3D data sets, segmentation can be an extremely difficult and laborious task, and thus has been one of the biggest hurdles for comprehensive 3D analysis. Heterogeneous catalysts have complex surface and bulk structures, and often sparse distribution of catalytic particles with relatively poor intrinsic contrast, which possess a unique challenge for image segmentation, including the current state-of-the-art deep learning methods. To tackle this problem, we apply a deep learning-based approach for the multi-class semantic segmentation of a γ-Alumina/Pt catalytic material in a class imbalance situation. Specifically, we used the weighted focal loss as a loss function and attached it to the U-Net’s fully convolutional network architecture. We assessed the accuracy of our results using Dice similarity coefficient (DSC), recall, precision, and Hausdorff distance (HD) metrics on the overlap between the ground-truth and predicted segmentations. Our adopted U-Net model with the weighted focal loss function achieved an average DSC score of 0.96 ± 0.003 in the γ-Alumina support material and 0.84 ± 0.03 in the Pt NPs segmentation tasks. We report an average boundary-overlap error of less than 2 nm at the 90th percentile of HD for γ-Alumina and Pt NPs segmentations. The complex surface morphology of γ-Alumina and its relation to the Pt NPs were visualized in 3D by the deep learning-assisted automatic segmentation of a large data set of high-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) tomography reconstructions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accuracy of ghost rotationally invariant slave-boson and dynamical mean field theory as a function of the impurity-model bath size

Here, we compare the accuracy of the ghost rotationally invariant slave-boson (g-RISB) theory and dynamical mean field theory (DMFT) on the single-band Hubbard model, as a function of the number of bath sites in the embedding impurity Hamiltonian. Our benchmark calculations confirm that the accuracy of g-RISB can be systematically improved by increasing the number of bath sites, similar to DMFT. With a few bath sites, we observe that g-RISB is systematically more accurate than DMFT for the ground-state observables. On the other hand, the relative accuracy of these methods is generally comparable for the quasiparticle weight and the spectral function. As expected, we observe that g-RISB satisfies the variational principle in infinite dimensions, as the total energy decreases monotonically towards the exact value as a function of the number of bath sites, suggesting that the g-RISB wave function may approach the exact ground state in infinite dimensions. Our results suggest that the g-RISB is a promising method for first-principles simulations of strongly correlated matter, which can capture the behavior of both static and dynamical observables, at a relatively low computational cost.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Bootstrapping outperforms community‐weighted approaches for estimating the shapes of phenotypic distributions

Abstract Estimating phenotypic distributions of populations and communities is central to many questions in ecology and evolution. These distributions can be characterized by their moments (mean, variance, skewness and kurtosis) or diversity metrics (e.g. functional richness). Typically, such moments and metrics are calculated using community‐weighted approaches (e.g. abundance‐weighted mean). We propose an alternative bootstrapping approach that allows flexibility in trait sampling and explicit incorporation of intraspecific variation, and show that this approach significantly improves estimation while allowing us to quantify uncertainty. We assess the performance of different approaches for estimating the moments of trait distributions across various sampling scenarios, taxa and datasets by comparing estimates derived from simulated samples with the true values calculated from full datasets. Simulations differ in sampling intensity (individuals per species), sampling biases (abundance, size), trait data source (local vs. global) and estimation method (two types of community‐weighting, two types of bootstrapping). We introduce the traitstrap R package, which contains a modular and extensible set of bootstrapping and weighted‐averaging functions that use community composition and trait data to estimate the moments of community trait distributions with their uncertainty. Importantly, the first function in the workflow, trait_fill , allows the user to specify hierarchical structures (e.g. plot within site, experiment vs. control, species within genus) to assign trait values to each taxon in each community sample. Across all taxa, simulations and metrics, bootstrapping approaches were more accurate and less biased than community‐weighted approaches. With bootstrapping, a sample size of 9 or more measurements per species per trait generally included the true mean within the 95% CI. It reduced average percent errors by 26%–74% relative to community‐weighting. Random sampling across all species outperformed both size‐ and abundance‐biased sampling. Our results suggest randomly sampling ~9 individuals per sampling unit and species, covering all species in the community and analysing the data using nonparametric bootstrapping generally enable reliable inference on trait distributions, including the central moments, of communities. By providing better estimates of community trait distributions, bootstrapping approaches can improve our ability to link traits to both the processes that generate them and their effects on ecosystems.

Maitner, Brian S.↗
DOE OSTI2021

Tuning dispersity of linear polymers and polymeric brushes grown from nanoparticles by atom transfer radical polymerization

Molecular weight distribution imposes considerable influence on the properties of polymers, making it an important parameter, impacting morphology and structural behavior of polymeric materials. Atom transfer radical polymerization (ATRP) has established itself as a powerful tool to prepare polymers with predetermined molecular weight, preserved chain-end functionality, and low dispersity. More recently, ATRP has also been shown to provide a means to deliberately broaden molecular weight distributions, and, via retaining living chain-ends, to enable the formation of block copolymers with designed block dispersity, featuring new microstructures and potentially attractive properties. Here, similar methodologies have been developed to facilitate tuning of the dispersity of polymeric brushes grown from nanoparticles thus resulting in hybrid materials with enhanced fracture toughness and high inorganic content. Recent advances have given access to brush architectures comprised of uni- and bimodal block copolymers with unique morphologies along with interesting mechanical, thermal, and optical properties.

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

Ion Mobilities, Transference Numbers, and Inverse Haven Ratios of Polymeric Ionic Liquids

We probe the ion mobilities, transference numbers, and inverse Haven ratio of ionic liquids and polymerized ionic liquids as a function of their molecular weight using a combination of atomistic equilibrium and nonequilibrium molecular dynamics simulations. In contrast to expectations, we demonstrate that the inverse Haven ratio increases with increasing degree of polymerization (N) and then decreases at larger N. For a fixed center of mass reference frame, we demonstrate that such results arise as a consequence of the strong cation–cation correlated motions, which exceed (in magnitude) the selfdiffusivity of cations. Together, our findings challenge the premise underlying the pursuit of pure polymeric ionic liquids as high transference number, single-ion conducting electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Robotic High-Throughput Grid-Search Platform for Mapping Phase Behavior in Triblock Copolymer–Homopolymer Blends

We present a high-throughput experimental investigation of the phase behavior in triblock copolymers (PS-b-PB-b-PS and PS-b-PI-b-PS) and polystyrene (PS) homopolymer blends as a function of homopolymer molecular weight (MW) and blend ratio. Using a robotic thin-film processing platform (NOVA) integrated with Grazing Incidence Small-Angle X-ray Scattering (GISAXS) and Atomic Force Microscopy (AFM), we systematically mapped the order–disorder transition (ODT) boundaries and domain spacing evolution across a broad MW range (4.0–101.3 kDa) with varying homopolymer loadings (10% to 90%). The results reveal three distinct regimes: low-MW homopolymers, corresponding to the wet-brush regime produced only gradual domain swelling before disordering at high blend ratios (weight fraction); medium-MW homopolymers, corresponding to thedry-brush regime induced significant domain spacing increase up to 80% followed by earlier disordering, while high-MW homopolymers led to macrophase separation with minimal changes in domain spacing. Additionally, coarse-grained molecular dynamics simulations confirmed our experimental finding that in the low-MW region, the PS homopolymer uniformly distributed in the PS domain. These findings demonstrate that homopolymer molecular weight critically governs both the extent of domain swelling and the onset of disorder in triblock copolymer systems. This high-throughput platform enables the rapid mapping of composition–morphology relationships and can be integrated with AI/ML tools for designing next-generation nanostructured polymers.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A variational framework for residual-based adaptivity in neural PDE solvers and operator learning

Residual-based adaptive strategies are widely used in scientific machine learning yet remain largely heuristic. We introduce a variational framework that formalizes these methods through convex transformations of the residual, where different transformations correspond to distinct objective functionals. For instance, exponential weights target uniform error minimization, while linear weights recover quadratic error minimization. This perspective reveals adaptive weighting as a means of selecting sampling distributions that optimize a primal objective, directly linking discretization choices to error metrics. This principled approach yields three key benefits: it enables systematic design of adaptive schemes, reduces discretization error by lowering estimator variance, and enhances learning dynamics by improving gradient signal-to-noise ratio. Extending the framework to operator learning, we demonstrate substantial performance gains across diverse optimizers and architectures. Our results provide a theoretical perspective for residual-based adaptivity and establish a foundation for principled discretization and training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Model of Hydrogen Solubility in Palladium-Silver Alloys

Fusion fuel cycle designs depend on palladium and palladium-silver alloys. The selective permeability of these is used to separate hydrogen isotopes from other elements. Predicting the solubility of hydrogen isotopes in palladium-silver alloys is important for the design of key unit operations for tritium processing. This paper presents a model of hydrogen isotope solubility in palladium and palladium-silver alloys based on Guggenheim’s quasichemical lattice theory. The model parameters were determined as functions of isotope molecular weight and the silver content of the alloy. The model predicts reasonable values for the solubility of protium, deuterium, and tritium over a wide range of temperatures, pressures, and alloy compositions.

08 HYDROGEN↗
DOE OSTI2021

Redshift space power spectrum beyond Einstein-de Sitter kernels

Abstract We develop a framework to compute the redshift space power spectrum (PS), with kernels beyond Einstein-de Sitter (EdS), that can be applied to a wide variety of generalized cosmologies. We build upon a formalism that was recently employed for standard cosmology in Chen, Vlah & White (2020), and utilize an expansion of the density-weighted velocity moment generating function that explicitly separates the magnitude of the k-modes and their angle to the line-of-sight direction dependencies. We compute the PS for matter and biased tracers to 1-loop Perturbation Theory (PT) and show that the expansion has a correct infrared and ultraviolet behavior, free of unwanted divergences. We also add Effective Field Theory (EFT) counterterms, necessary to account for small-scale contributions to PT, and employ an IR-resummation prescription to properly model the smearing of the BAO due to large scale bulk flows within Standard-PT. To demonstrate the applicability of our formalism, we apply it on the ΛCDM and the Hu-Sawicki f(R) models, and compare our numerical results against the elephant suite of N-body simulations, finding very good agreement up to k = 0.27 Mpc -1 h at z = 0.5 for the first three non-vanishing Legendre multipoles of the PS. To our knowledge, the model presented in this work is the most accurate theoretical EFT-PT for modified gravity to date, being the only one that accounts for beyond linear local biasing in redshift-space. Hence, we argue our RSD modeling is a promising tool to construct theoretical templates in order to test deviations from ΛCDM using real data obtained from the next stage of cosmological surveys such as DESI and LSST.

Astronomy & Astrophysics↗
DOE OSTI2022

A Q-Learning-Based Approach for Distributed Beam Scheduling in mmWave Networks

We consider the problem of distributed downlink beam scheduling and power allocation for millimeter-Wave (mmWave) cellular networks where multiple base stations (BSs) belonging to different service operators share the same unlicensed spectrum with no centralized control or explicit coordination among them. Our goal is to design efficient distributed beam scheduling and power allocation schemes such that the network-level payoff, defined as the weighted sum of the total throughput and a power penalization term, can be maximized. To this end, we propose a distributed scheduling approach to power allocation and adaptation for efficient interference management over the shared spectrum by modeling each BS as an independent Q-learning agent. As a baseline, we compare the proposed approach to the state-of-the-art non-cooperative game-based approach which was previously developed for the same problem. We conduct extensive experiments under various scenarios to verify the effect of multiple factors on the performance of both approaches. Experiment results show that the proposed approach adapts well to different interference conditions by learning from experience and thus can achieve higher payoff than the game-based approach. The proposed approach can also be integrated into our previously developed Lyapunov stochastic optimization framework for the purpose of network utility maximization with optimality guarantee. As a result, the weights in the payoff function can be automatically and optimally determined by the virtual queue values from the sub-problems derived from the Lyapunov optimization framework.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

Accelerated Aging of Fast Pyrolysis Bio-Oil Using Carbonyl Titration: Laboratory Analytical Procedure (LAP)

This laboratory analytical procedure covers the accelerated aging of fast pyrolysis bio-oils. Bio-oils undergo reactions that result in physical and chemical changes over time. These changes typically result in an increase in molecular weight, decrease in some functional groups such as carbonyls, and an increase in viscosity; additionally, the aging process often leads to phase separation. Studies have shown that accelerated aging of bio-oils using this method closely mimics room temperature aging for long periods of time (over 3 years). This procedure has been developed for the accelerated aging of fast pyrolysis bio-oils only. Fast pyrolysis bio-oils are more reactive than other bio-oils (e.g., catalytic fast pyrolysis) and undergo aging much more rapidly. Therefore, this procedure should only be used for fast pyrolysis bio-oil samples.

09 BIOMASS FUELS↗
DOE OSTI2024

The mechanical properties of Kel-F 800 (FK-800) as a function of crystallinity

Kel-F 800 is a copolymer of chlorotrifluoroethylene PTFE (75 wt. %) and vinylidene fluoride PVDF (25 wt. %). It has previously been used as a PBX binder for insensitive explosives such as PBX 9502 and LX-17. 3M started production of Kel-F 800 in 1957 and small-scale batches continued to be made until 2002 when production ceased due to environmental concerns regarding one of the emulsifiers used during production. Around 2000 the Kel-F 800 name was changed to FK-800 to avoid trademark concerns because rights to produce another polymer with a similar tradename (Kel-F 81) had been sold to another manufacturer. The Kel designation came from the original manufacturer of PCTFE (Kel-F 81), the Kellog company. To avoid confusion this document will only refer to Kel-F 800. In 2006, production of small-scale batches of Kel-F 800 was started again by 3M in response to customer enquiries. This new material, the first blended batch is referred to as LOT 1, was produced with a different emulsifier than used previously. Because Kel-F 800 is made in a small batch reactor, considerable variation in crystallinity can be expected from lot to lot and year to year. In many ways, this is not significant since the material is dissolved in a solvent (often MEK, ethylmethyl ketone or ethyl acetate) for PBX production purposes. This destroys the as received crystallinity and the resulting crystallinity in the processed material is a function of polymer molecular weight and thermal history. Producing large billets of Kel-F 800 from solvent extraction is not practical and so a compression molding technique has been used above the melting temperature. This method also removes residual crystallinity from the supplied granules. The molecular weight of a polymer can be estimated by several techniques, the most common being gel permittivity chromatography (GPC), size exclusion chromatography (SEC) and shear rheometry measurements of polymer/solvent solutions. Changes in molecular weight will affect the crystallization rate and the maximum crystallinity reached for a specific thermal history. Both references agree that the new LOT 1 material molecular weight falls within the deviation found from averaging previous historical lots of Kel-F 800.

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

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
DOE OSTI2020

Spectroscopic signatures of next-nearest-neighbor hopping in the charge and spin dynamics of doped one-dimensional antiferromagnets

In this work, we study the impact of next-nearest-neighbor (nnn) hopping on the low-energy collective excitations of strongly correlated doped antiferromagnetic cuprate spin chains. Specifically, we use exact diagonalization and the density matrix renormalization group method to study the single-particle spectral function, the dynamical spin and charge structure factors, and the Cu L-edge resonant inelastic x-ray scattering (RIXS) intensity of the doped t-t'-J model for a set of t' values. We find evidence for the breakdown of spin-charge separation as |t'| increases and identify its fingerprints in the dynamical response functions. The inclusion of nnn hopping couples the spinon and holon excitations, resulting in the formation of a spin-polaron, where a ferromagnetic spin polarization cloud dresses the doped carrier. The spin-polaron manifests itself as additional spectral weight in the dynamical correlation functions, which appear simultaneously in the spin- and charge-sensitive channels. We also demonstrate that RIXS can provide a unique view of the spin-polaron, due to its sensitivity to both the spin and charge degrees of freedom.

1D spin chains↗
DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗