Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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

A novel formulation of the unintegrated gluon distribution for DIS

We provide a semi-classical description of the inclusive gluon mediated Deep Inelastic Scattering cross section in a way that accounts for the leading powers in both the Regge and Bjorken limits. Our approach thus allows a systematic matching of small and moderate xBj regimes of gluon proton structure functions. We find a new unintegrated gluon distribution with an explicit dependence on the longitudinal momentum fraction x which entirely spans both the dipole operator and the gluonic Parton Distribution Function. Computing this gauge invariant gluon operator on the lattice could allow to probe the energy dependence of the saturation scale from first principles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Liquid Structure of CO 2 –Reactive Aprotic Heterocyclic Anion Ionic Liquids from X-ray Scattering and Molecular Dynamics

A combination of X-ray scattering experiments and molecular dynamics simulations were conducted to investigate the structure of ionic liquids (ILs) which chemically bind CO 2 . The structure functions were measured and computed for four different ILs consisting of two different phosphonium cations, triethyloctylphosphonium ([P 2228 ] + ) and trihexyltetradecylphosphonium ([P 66614 ] + ), paired with two different aprotic heterocyclic anions which chemically react with CO 2 , 2-cyanopyrrolide, and 1,2,4-triazolide. Simulations were able to reproduce the experimental structure functions, and by deconstructing the simulated structure functions, further information on the liquid structure was obtained. All structure functions of the ILs studied had three primary features which have been seen before in other ILs: a prepeak near 0.3–0.4 Å–1 corresponding to polar/nonpolar domain alternation, a charge alternation peak near 0.8 Å–1, and a peak near 1.5 Å–1 due to interactions of adjacent molecules. The liquid structure functions were only mildly sensitive to the specific anion and whether or not they were reacted with CO 2 . Upon reacting with CO 2 , small changes were observed in the structure functions of the [P 2228 ] + ILs, whereas virtually no change was observed upon reacting with CO 2 in the corresponding [P 66614 ] + ILs. When the [P 2228 ] + cation was replaced with the [P 66614 ] + cation, there was a significant increase in the intensities of the prepeak and adjacency interaction peak. While many of the liquid structure functions are similar, the actual liquid structures differ as demonstrated by computed spatial distribution functions.

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

Energetics and Structure of Ag–Water Clusters Formed in Mordenite

Zeolite-supported Ag 0 clusters have broad applications from catalysis to medicine, necessitating a mechanistic understanding of the formation of Ag 0 clusters in situ. Density functional theory (DFT) simulations have been performed on silver, water, and silver–water clusters in silica mordenite (Si-MOR), to identify the role of the confinement on the structure and energetics of Ag 0 cluster formation. The most favorable binding energy in the 12-membered ring (MR) pore of the Si-MOR is a 10–15-atom Ag 0 cluster. Computational pair distribution function (PDF) data indicates that the Ag 0 and Ag 0 –H 2 O clusters formed in vacuum versus in Si-MOR exhibit structural differences. Additionally, when the Ag 0 cluster is confined, the density decreases and the surface area increases, hypothesized to be due to the limiting geometry of the 12-MR main channel. An energetic drive toward formation of larger Ag 0 clusters was also identified, with hydrated silver atoms generating higher energy structures. In conclusion, this work identifies mechanistic and structural insight into the role of nanoconfinement on formation of Ag 0 clusters in mordenite.

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

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

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

Consideration of the Aluminum Distribution in Zeolites in Theoretical and Experimental Catalysis Research

Research efforts in zeolite catalysis have become increasingly cognizant of the diversity in structure and function resulting from the distribution of framework aluminum atoms, through emerging reports of catalytic phenomena that fall outside those recognizable as the shape-selective ones emblematic of its earlier history. Molecular-level descriptions of how active-site distributions affect catalysis are an aspirational goal articulated frequently in experimental and theoretical research, yet they are limited by imprecise knowledge of the structure and behavior of the zeolite materials under interrogation. In experimental research, higher precision can result from more reliable control of structure during synthesis and from more robust and quantitative structural and kinetic characterization probes. In theoretical research, construction of models with specific aluminum locations and distributions seldom capture the heterogeneity inherent to the materials studied by experiment. In this Perspective, we discuss research findings that appropriately frame the challenges in developing more predictive synthesis-structure-function relations for zeolites, highlighting studies on ZSM-5 zeolites that are among the most structurally complex molecular sieve frameworks and the most widely studied because of their versatility in commercial applications. We discuss research directions to address these challenges and forge stronger connections between zeolite structure, composition, and active sites to catalytic function. Such connections promise to aid in bridging the findings of theoretical and experimental catalysis research, and transforming zeolite active site design from an empirical endeavor into a more predictable science founded on validated models.

09 BIOMASS FUELS↗
DOE OSTI2016

Efficient photon coupling from a diamond nitrogen vacancy center by integration with silica fiber

A central goal in quantum information science is to efficiently interface photons with single optical modes for quantum networking and distributed quantum computing. Here, we introduce and experimentally demonstrate a compact and efficient method for the low-loss coupling of a solid-state qubit, the nitrogen vacancy (NV) center in diamond, with a single-mode optical fiber. In this approach, single-mode tapered diamond waveguides containing exactly one high quality NV memory are selected and integrated on tapered silica fibers. Numerical optimization of an adiabatic coupler indicates that near-unity-efficiency photon transfer is possible between the two modes. Experimentally, we find an overall collection efficiency between 16% and 37% and estimate a single photon count rate at saturation above 700 kHz. This integrated system enables robust, alignment-free, and efficient interfacing of single-mode optical fibers with single photon emitters and quantum memories in solids.

Optics↗
DOE OSTI2023

Deducing subnanometer cluster size and shape distributions of heterogeneous supported catalysts

Abstract Infrared (IR) spectra of adsorbate vibrational modes are sensitive to adsorbate/metal interactions, accurate, and easily obtainable in-situ or operando. While they are the gold standards for characterizing single-crystals and large nanoparticles, analogous spectra for highly dispersed heterogeneous catalysts consisting of single-atoms and ultra-small clusters are lacking. Here, we combine data-based approaches with physics-driven surrogate models to generate synthetic IR spectra from first-principles. We bypass the vast combinatorial space of clusters by determining viable, low-energy structures using machine-learned Hamiltonians, genetic algorithm optimization, and grand canonical Monte Carlo calculations. We obtain first-principles vibrations on this tractable ensemble and generate single-cluster primary spectra analogous to pure component gas-phase IR spectra. With such spectra as standards, we predict cluster size distributions from computational and experimental data, demonstrated in the case of CO adsorption on Pd/CeO 2 (111) catalysts, and quantify uncertainty using Bayesian Inference. We discuss extensions for characterizing complex materials towards closing the materials gap.

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

Geospatial mapping of distribution grid with machine learning and publicly-accessible multi-modal data

Abstract Detailed and location-aware distribution grid information is a prerequisite for various power system applications such as renewable energy integration, wildfire risk assessment, and infrastructure planning. However, a generalizable and scalable approach to obtain such information is still lacking. In this work, we develop a machine-learning-based framework to map both overhead and underground distribution grids using widely-available multi-modal data including street view images, road networks, and building maps. Benchmarked against the utility-owned distribution grid map in California, our framework achieves > 80% precision and recall on average in the geospatial mapping of grids. The framework developed with the California data can be transferred to Sub-Saharan Africa and maintain the same level of precision without fine-tuning, demonstrating its generalizability. Furthermore, our framework achieves a R 2 of 0.63 in measuring the fraction of underground power lines at the aggregate level for estimating grid exposure to wildfires. We offer the framework as an open tool for mapping and analyzing distribution grids solely based on publicly-accessible data to support the construction and maintenance of reliable and clean energy systems around the world.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Note: Interionic Potentials of Mean Force for Ca2+-Cl- in Polarizable Water

We describe our molecular dynamics simulation studies of the inter-ionic potentials of mean force (PMF) for the Ca{sup 2+}-Cl{sup -} ion pair in water. Ion-water and ion-ion polarizable potential models were constructed to reproduce experimental ion-hydration enthalpies as well as the electronic-structure calculations of the minimum energy for the Ca{sup 2+}-Cl{sup -} interactions. Our study predicted a PMF with no contact ion-pair state and a shallow solvent-separated ion-pair state. In addition, our computed radial distribution function, g{sub Ca-Cl}, for a 1-M CaCl{sub 2} electrolyte solution indicated that there is no contact-ion pair but, instead, showed a maximum near 5 {angstrom}. This observation is consistent with our PMF calculation. These new results are discussed and compared to the recent classical and ab initio PMF calculations on the same system reported by Timko, De Castro, and Kuyucak (J. Chem. Phys. 134, 204510 (2011)).

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

On a theory of stability for nonlinear stochastic chemical reaction networks

We present elements of a stability theory for small, stochastic, nonlinear chemical reaction networks. Steady state probability distributions are computed with zero-information (ZI) closure, a closure algorithm that solves chemical master equations of small arbitrary nonlinear reactions. Stochastic models can be linearized around the steady state with ZI-closure, and the eigenvalues of the Jacobian matrix can be readily computed. Eigenvalues govern the relaxation of fluctuation autocorrelation functions at steady state. Autocorrelation functions reveal the time scales of phenomena underlying the dynamics of nonlinear reaction networks. In accord with the fluctuation-dissipation theorem, these functions are found to be congruent to response functions to small perturbations. Significant differences are observed in the stability of nonlinear reacting systems between deterministic and stochastic modeling formalisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Structure and thermodynamics of core-softened models for alcohols

The phase behavior and the fluid structure of coarse-grain models for alcohols are studied by means of reference interaction site model (RISM) theory and Monte Carlo simulations. Specifically, we model ethanol and 1-propanol as linear rigid chains constituted by three (trimers) and four (tetramers) partially fused spheres, respectively. Thermodynamic properties of these models are examined in the RISM context, by employing closed formulæ for the calculation of free energy and pressure. Gas-liquid coexistence curves for trimers and tetramers are reported and compared with already existing data for a dimer model of methanol. Critical temperatures slightly increase with the number of CH{sub 2} groups in the chain, while critical pressures and densities decrease. Such a behavior qualitatively reproduces the trend observed in experiments on methanol, ethanol, and 1-propanol and suggests that our coarse-grain models, despite their simplicity, can reproduce the essential features of the phase behavior of such alcohols. The fluid structure of these models is investigated by computing radial distribution function g{sub ij}(r) and static structure factor S{sub ij}(k); the latter shows the presence of a low−k peak at intermediate-high packing fractions and low temperatures, suggesting the presence of aggregates for both trimers and tetramers.

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

Polarization response of clathrate hydrates capsulated with guest molecules

Clathrate hydrates are characterized by their water cages encapsulating various guest atoms or molecules. The polarization effect of these guest-cage complexes was studied with combined density functional theory and finite-field calculations. An addition rule was noted for these systems whose total polarizability is approximately equal to the polarizability sum of the guest and the cage. However, their distributional polarizability computed with Hirshfeld partitioning scheme indicates that the guest–cage interaction has considerable influence on their polarization response. The polarization of encapsulated guest is reduced while the polarization of water cage is enhanced. The counteraction of these two opposite effects leads to the almost unchanged total polarizability. Further analysis reveals that the reduced polarizability of encapsulated guest results from the shielding effect of water cage against the external field and the enhanced polarizability of water cage from the enhanced bonding of hydrogen bonds among water molecules. Although the charge transfer through the hydrogen bonds is rather small in the water cage, the polarization response of clathrate hydrates is sensitive to the changes of hydrogen bonding strength. The guest encapsulation strengthens the hydrogen bonding network and leads to enhanced polarizability.

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

Equilibrium reconstruction with 3D eddy currents in the Lithium Tokamak eXperiment

Axisymmetric free-boundary equilibrium reconstructions of tokamak plasmas in the Lithium Tokamak eXperiment (LTX) are performed using the PSI-Tri equilibrium code. Reconstructions in LTX are complicated by the presence of long-lived non-axisymmetric eddy currents generated by a vacuum vessel and first wall structures. To account for this effect, reconstructions are performed with additional toroidal current sources in these conducting regions. The eddy current sources are fixed in their poloidal distributions, but their magnitude is adjusted as part of the full reconstruction. Eddy distributions are computed by toroidally averaging currents, generated by coupling to vacuum field coils, from a simplified 3D filament model of important conducting structures. The full 3D eddy current fields are also used to enable the inclusion of local magnetic field measurements, which have strong 3D eddy current pick-up, as reconstruction constraints. Using this method, equilibrium reconstruction yields good agreement with all available diagnostic signals. Here, an accompanying field perturbation produced by 3D eddy currents on the plasma surface with a primarily n = 2, m = 1 character is also predicted for these equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

MATCHING THE SPECTRAL ENERGY DISTRIBUTION AND p-MODE OSCILLATION FREQUENCIES OF THE RAPIDLY ROTATING DELTA SCUTI STAR {alpha} OPHIUCHI WITH A TWO-DIMENSIONAL ROTATING STELLAR MODEL

Spectral energy distributions are computed using two-dimensional (2D) rotating stellar models and non-LTE plane-parallel model atmospheres. A rotating, 2D stellar model has been found that matches the observed ultraviolet and visible spectrum of {alpha} Oph. The SED match occurs for the interferometrically deduced surface shape and inclination, and is different from the SED produced by spherical models. The p-mode oscillation frequencies in which the latitudinal variation is modeled by a linear combination of eight Legendre polynomials were computed for this model. The five highest and seven of the nine highest amplitude modes show agreement between computed axisymmetric, equatorially symmetric mode frequencies and the mode frequencies observed by the Microvariability and Oscillations of Stars satellite (MOST) to within the observational error. Including nonaxisymmetric modes up through |m| = 2, and allowing for the possibility that the eight lowest amplitude modes could be produced by modes that are not equatorially symmetric, matches for 24 out of the 35 MOST modes to within the observational error and another eight modes to within twice the observational error. The remaining three observed modes can be fitted within 4.2 times the observational error, but even these may be fitted to within the observational error if the criteria for computed modes are expanded.

79 ASTRONOMY AND ASTROPHYSICS↗