Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computing 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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Loop-by-loop differential equations for dual (elliptic) Feynman integrals

We present a loop-by-loop method for computing the differential equations of Feynman integrals using the recently developed dual form formalism. We give explicit prescriptions for the loop-by-loop fibration of multi-loop dual forms. Then, we test our formalism on a simple, but non-trivial, example: the two-loop three-mass elliptic sunrise family of integrals. We obtain an ε-form differential equation within the correct function space in a sequence of relatively simple algebraic steps. In particular, none of these steps relies on the analysis of q-series. Then, we discuss interesting properties satisfied by our dual basis as well as its simple relation to the known ε-form basis of Feynman integrands. The underlying K3-geometry of the three-loop four-mass sunrise integral is also discussed. Finally, we speculate on how to construct a “good” loop-by-loop basis at three-loop.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Functional methods for Heavy Quark Effective Theory

We use functional methods to compute one-loop effects in Heavy Quark Effective Theory. The covariant derivative expansion technique facilitates the efficient extraction of matching coefficients and renormalization group evolution equations. This paper pro- vides the first demonstration that such calculations can be performed through the algebraic evaluation of the path integral for the class of effective field theories that are (i) constructed using a non-trivial one-to-many mode decomposition of the UV theory, and (ii) valid for non-relativistic kinematics. We discuss the interplay between operators that appear at intermediate steps and the constraints imposed by the residual Lorentz symmetry that is encoded as reparameterization invariance within the effective description. The tools presented here provide a systematic approach for computing corrections to higher order in the heavy mass expansion; precision applications include predictions for experimental data and connections to theoretical tests via lattice QCD. A set of pedagogical appendices comprehensively reviews modern approaches to performing functional calculations algebraically, and derives contributions from a term with open covariant derivatives for the first time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Normalization of ZZ instanton amplitudes in type 0B minimal superstring theory

Abstract We study ZZ instanton corrections in the (2,4k)$$ \mathcal{N} $$ N = 1 minimal superstring theory with the type 0B GSO projection, which becomes the type 0B$$ \mathcal{N} $$ N = 1 super-JT gravity in thek→ ∞ limit. Each member of the (2,4k) family of theories has two phases distinguished by the sign of the Liouville bulk cosmological constant. The worldsheet method for computing the one-loop normalization constant multiplying the instanton corrections gives an ill-defined answer in both phases. We fix these divergences using insights from string field theory and find finite, unambiguous results. Each member of the (2,4k) family of theories is dual to a double-scaled one-matrix integral, where the double-scaling limit can be obtained starting either from a unitary matrix integral with a leading one-cut saddle point, or from a hermitian matrix integral with a leading two-cut saddle point. The matrix integral exhibits a gap-closing transition, which is the same as the double-scaled Gross-Witten-Wadia transition whenk= 1. We also compute instanton corrections in the double-scaled matrix integral for allkand in both phases, and find perfect agreement with the string theory results.

Physics↗
DOE OSTI2020

BTZ one-loop determinants via the Selberg zeta function for general spin

We relate the heat kernel and quasinormal mode methods of computing the 1-loop partition function of arbitrary spin fields on a rotating (Euclidean) BTZ background using the Selberg zeta function associated with H 3 /Z, extending (arXiv:1811.08433) [1]. Previously, Perry and Williams [2] showed for a scalar field that the zeros of the Selberg zeta function coincide with the poles of the associated scattering operator upon a relabeling of integers. We extend the integer relabeling to the case of general spin, and discuss its relationship to the removal of non-square-integrable Euclidean zero modes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Effect of disorder on thermodynamic instability of binary Rare-earth – Nickel – Palladium compounds

In this work, we have investigated the thermodynamic stability of disordered rare-earth phases SmX 2 and Sm 10 X 21 (X=Ni, Pd) using machine-learning based analytical descriptor and first-principles density functional theory methods. The absence of Laves phase compounds in R-Pd binary systems is a longstanding problem of rare earth science: even though Ni and Pd belong to the same group of the periodic table and have similar electronic structure, the Pd compound crystallizes in a monoclinic (C2/m) phase with 10:21 stoichiometry, i.e., Sm 10 Pd 21 , while the Ni compound adopts a cubic Laves phase (MgCu 2 ) structure. To understand this contrasting phase stability, we performed thermodynamic convex hull analysis of Sm x Ni 1-x and Sm x Pd 1-x binary systems, which is experimentally validated using powder X-ray diffraction (PXRD) analyzes of polycrystalline Sm(Ni x Pd 1-x ) 2 samples with x=0, 0.5, and 1. A detailed electronic-structure (band-structure, charge density, and Fermi-surface) analysis of the differences between SmNi 2 /SmPd 2 and Sm 10 Ni 21 /Sm 10 Pd 21 compounds provides the quantum mechanical origin of the unfavorable mixing of Pd with Ni in cubic Laves phase. We show that the stability of Sm-Pd in 10:21 stoichiometry arises from improved intra-/inter-layer 5d-4d bonding compared to the 1:2 stoichiometry. Our work emphasizes the importance of ab-initio methods and computationally inexpensive analytical descriptors for the detailed analysis of thermodynamic and electronic properties of hard-to-prepare rare-earth compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of applied potential on metal surfaces: Surface energy, Wulff shape and charge distribution

Here we use grand canonical density functional theory to predict the surface energies, Wulff shapes, charge distributions and catalytically active sites of different metal surfaces under electrochemical conditions. We propose a method for computing surface energies from grand canonical density functional theory (GC-DFT) calculations of periodic slab models and use it to compute the surface energies of the facets of Pt, Cu, and Ag crystals to predict their Wulff shapes under electrochemical conditions. GC-DFT predicts that, for the pure metals studied, solvation only slightly affects the Wulff shape while applied potentials considerably affect the surface energies and corresponding Wulff shapes. We used Bader charge analysis of GC-DFT computed electron densities to investigate the effect of applied potential on the distribution of electron density over the atoms of the surfaces of Pt, Cu, Ag, and the 75–25 Ag-Pt and Au-Ni alloys. This analysis shows that, under an applied potential, the electron density is unevenly distributed over the surface atoms and that the charges of atoms more exposed to solvent are more sensitive to bias. Our results show that the most sensitive atom to bias can be used to identify the most favorable adsorption site and thus, the active sites of electrochemical reactions, which is computationally less demanding than calculating the adsorption energies on all possible adsorption sites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Investigation of countercurrent flow profile and liquid holdup in random packed column with local CFD data

Liquid holdup and mass transfer area are critical parameters used in packed column design and CO 2 capture efficiency prediction. Compared to experimental methods, 3D computational fluid dynamics (CFD) can provide high-resolution data for studying countercurrent flow hydrodynamics in a packed column. In this paper, a framework was established for modeling the liquid–gas countercurrent flow in a random packed column with pall rings. Besides the column-averaged information, the radial pall ring distribution, velocity, and liquid holdup profiles are obtained in the simulations. The radial two-phase velocity and liquid holdup profiles show significant fluctuations caused by the non-uniform pall ring distribution. The entrance effect and wall influence are also observed in the packed column. By using local countercurrent flow information, the packing specific area $a_p$ and liquid velocity $u_L$ range are significantly expanded for liquid holdup hL prediction. This enables us to develop the liquid holdup correlation with one packing size and six CFD runs only. The proposed power law of the liquid holdup correlation $h_L~u_L^{0.44}$ indicates the random packed column falls in a viscous to turbulent transition regime with a Reynolds Number range of [6.7–40.2]. The derived liquid holdup correlation is in good agreement with existing correlations that are developed using the column-averaged experimental data.

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

Assessing Mg–Sc–(rare earth) ternary phase stability via constituent binary cluster expansions

The disordered Mg–Sc body-centered cubic (bcc) phase is both lightweight and strong; however, the system is impractical for general industrial use due to the high cost of scandium. Here we propose a computationally efficient metric that assesses ternary rare earth element additions that may stabilize the bcc phase at lower Sc concentrations. We find that the bcc phase is stabilized by the ternary addition of Y or Er, but not by La, Ce, or Nd, and we validate these predictions by experimental production and characterization of Mg–Sc–(Y,Er,Nd) alloys. The results suggest a computationally efficient method to anticipate integration of ternary elements into binary systems using cluster expansions of constituent binaries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Challenges for the statistical gravitational-wave method to measure the Hubble constant

Gravitational waves (GW) can be employed as standard sirens that will soon measure the Hubble constant with sufficient precision to weigh in on the ~5σ Hubble tension. Most GW sources will have no identified electromagnetic counterpart, leading to uncertainty in the redshift of the source, and in turn a degeneracy between host galaxy distance, redshift, and H 0 . In the case where no electromagnetic counterparts are identified, it has been proposed that a statistical canvassing of candidate GW hosts, found in a large galaxy survey for example, can be used to accurately constrain the Hubble constant. We study and simulate this "galaxy voting" method to compute H 0 . Here, we find that the Hubble constant posterior is in general biased relative to the true value even when making optimistic assumptions about the statistical properties of the sample. Using the MICECAT light-cone catalog, we find that the bias in the H0 posteriors depends on the realization of the underlying galaxy sample and the precision of the GW source distance measurement.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Power allocation by load aggregator with heterogeneous loads using weighted projection

With the recently increasing adoption of transactive energy markets, the role of a load aggregator is becoming increasingly important. Among the many tasks of an aggregator, the task of managing the electricity consumption among a group of grid-interactive efficient buildings according to the demand on the grid is taken up. Specifically, here, the task of allocating a certain amount of power among the participating customers while respecting their preferences and comfort is considered. Modeling this as an aggregate flexibility, the virtual battery conditions that are employed to verify the feasibility of a candidate power profile are derived. A feasible power profile, when allocated, satisfies specified comfort bands. To determine the ideal power requirement of the loads, the model-free control method is employed. The feasible power profile is allocated by embedding the customers’ preferences through weights, posing the power constraint as a hyperplane, and employing a weighted projection from the ideal power requirements onto the constraint plane to solve the considered task. The proposed method is computationally efficient and scalable to any number of heterogeneous thermostatically controlled loads in buildings.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Interfaces in all solid state Li-metal batteries: A review on instabilities, stabilization strategies, and scalability

With technological advancements in electrochemical energy storage systems increasing at a spectacular rate, batteries equipped with a lithium anode hold the key towards unlocking high energy densities. While lithium-ion batteries with layered anodes (e.g. graphite) and liquid organic electrolytes have been ubiquitous in portable electronics, electric vehicles, and grid applications, all solid-state batteries that use the combination of a lithium anode and a solid-state electrolyte (SSE) will further advance the present technology. The underlying challenge that limits the successful development of all solid-state batteries (ASSBs) is dictated largely by the highly reactive interfaces at the anode/SSE and the cathode/SSE interface. In this comprehensive review, we present an overview of the following: (i) characterization of the electrode/SSE interface via multimodal characterization, which include X-ray-, electron-, neutron-, optical-, and computation-based methods: (ii) parasitic reactions that occur from chemical, mechanical or electrochemical instabilities and interfacial engineering strategies for improving the stability of these interfaces classified by the class of inorganic SSEs (sulfide-, NASICON-, and garnet-type SSEs); (iii) laboratory-to-industry scale processing perspectives of SSEs; (iv) scalability and manufacturing aspects of current interfacial strategies; and (v) the prospects of all ASSBs within the context of extreme fast charging capability. Here, this review seeks to highlight key efforts in the field of ASSBs, by focusing particularly on stabilizing the electrode/SSE interfaces, which will help to bridge fundamental studies to technological relevance.

25 ENERGY STORAGE↗
DOE OSTI2020

An accelerated approach to inline non-LTE modeling

A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Feasibility study for test rig assessments of fish passage conditions in a Kaplan turbine

The assessment of fish passage conditions in hydroelectric turbines consists of identifying and quantifying physical magnitudes leading to increased risks of injury of fish passing through turbines in operation. Such assessments are usually carried out either with the use of computer-based methods during design or with field testing of live fish and sensors passing through prototypes. A method in between consists of test rig experimentation, which is critical for testing fish-focused design concepts and offers the opportunity for implementing the most effective design measures for improved fish survivability. However, fish-related assessments in test rigs are not sufficiently documented for industrial applications. This work presents the main findings of an experimental campaign to quantify fish-related hydraulic magnitudes in a physical model of a Kaplan turbine in a commercial test rig. Two operating conditions were tested by releasing miniaturized autonomous sensor devices (termed Sensor Fish Mini) at the turbine intake flow, passing them through the runner in motion and recovering them at the draft tube exit. During passage, time series of acceleration, absolute pressure and rotational velocity were recorded. The recordings were then interpreted to determine the magnitude and likely location of hydraulic stressors hazardous to fish. The statistical tests on the reported measurements indicated that low pressure, collisions on the runner and rotations in the draft tube were not different between the two tested operating points. On the other hand, pressure drop and collision rates on the distributor differed considerably as a function of net head. The outcomes of this investigation showed that test rig evaluations of fish-related properties with Sensor Fish Mini can contribute to an evidence-based development of turbine geometries designed for providing safer passage conditions. Future work will investigate the scaling of test rig measurements to hydraulically equivalent magnitudes in the prototype and their biological consequences.

13 HYDRO ENERGY↗
DOE OSTI2022

Efficient high-fidelity TRISO statistical failure analysis using Bison: Applications to AGR-2 irradiation testing

The ability of tri-structural isotropic (TRISO) fuel to contain fission products is largely dictated by the quality of the manufacturing process, since most of the fission product release is expected to occur due to coating layer failure in a small number of particles containing defects. The Bison fuel performance code has capabilities to predict failure in individual particles, accounting for the presence of defects, and to apply statistical analysis methods to compute the probability of failure in a set of fuel particles. Bison has recently undergone significant development both to improve its physical representations of fuel particle behavior and to improve the efficiency of its statistical failure calculations. Physical model improvements include new capabilities to account for the pressure generated by fission gases on inner pyrolytic carbon (IPyC) crack surfaces and to use local material coordinate orientation to accurately incorporate the anisotropy in the material properties in aspherical particles. To improve statistical modeling efficiency, a direct integration approach which involves directly integrating the failure probability function associated with statistically varying parameters has been developed. The direct integration approach is much more efficient than the Monte Carlo (MC) schemes commonly employed, and allows Bison to directly run high-dimensional fuel performance models, which improves the accuracy of failure probability calculations. Finally, a set of benchmark problems is considered here to compare the MC and direct integration approaches, and a statistical failure analysis of compacts in the Advanced Gas Reactor (AGR)-2 experiments is performed using the direct integration approach.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗
DOE OSTI2024

The carbon challenge: Design, synthesis, and chemisorption behavior of solid sorbents in direct air capture of carbon dioxide

Direct air capture (DAC) of CO 2 is a promising solution for reducing the carbon footprint through "negative emission" technology. However, the low CO 2 concentration (~400 ppm) and the dynamic nature of DAC processes present challenges in designing effective sorbent systems. Recent advancements in material design and structural engineering have led to the development of high-performance solid sorbents, offering a more stable, safe, and energy-efficient alternative to traditional liquid CO 2 capture methods. This review highlights progress in solid sorbent-based DAC, focusing on amine-modified materials, hydrogen-bonded frameworks, and ionic liquid-engineered scaffolds. The discussion covers design principles, synthesis methodologies, and their impact on CO 2 chemisorption, comparing the advantages and limitations of each approach. Characterization techniques, especially operando methods and computational tools, are reviewed to understand sorbent behavior during CO 2 integration and release. The reaction pathways and interaction mechanisms of these sorbents with CO 2 are analyzed to guide future design. Additionally, the CO 2 chemisorption behaviors, including capacity, sorption kinetics, recyclability, and durability in the presence of gaseous impurities and under humid conditions will be evaluated and compared. Further, the review offers unique insights into the physical properties, chemical structures, and surface engineering effects of these sorbents, based on comprehensive characterization and evaluation techniques.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗