Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “operator theory”

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

High Performance SOFCs with a Superior Stability for Reliable and Durable Power Systems

Next generation of fuel cells, electrolyzers, and batteries requires higher power, faster kinetics, and larger energy density, which necessitate the use of compositionally complex oxides to achieve multifunctionalities and activity. These compositionally complex oxides may change their phases and structures during an electrochemical process – a so-called “electrochemically driven phase transformation”. The origin for such a phase change has remained obscure. More importantly, there is a need to develop high performance solid oxide fuel cells with an enhanced stability. In this work, the La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3 (LSCF) cathode surface is modified by infiltration of Pr 6 O 11 and the power density at 0.8V and 750 °C is improved by 21%. In addition, by replacing the traditional barrier layer Gd 0.2 Ce 0.8 O 1.9 with mixed conducting Pr 0.1 Gd 0.1 Ce 0.8 O 1.9 , the power density increases by 38%. The different mechanism of promotions was investigated by electrochemical impedance spectroscopy. The ohmic resistance is dramatically reduced by applying the PGCO interlayer, and the distribution of relaxation time was used to analyze the mechanism for which the polarization resistance was decreased attributing to the mixed conduction nature in PrO x . An increase of power density at 0.8 V of 0.358 W/cm 2 (71%) is achieved with the implementation of both surface modification and buffer layer engineering. An experimental study and a theoretical analysis were then carried out on phase evolution in praseodymium nickelates. Nickelate-based electrodes show up to 60× greater phase transformation during operation when compared to thermally annealed ones. Theoretical analysis suggests that the presence of a reduced oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is the origin for the phase change in an oxygen electrode. Guided by the theory, an addition of the electronic conduction in the interface layer leads to the significant suppression of phase change, while improving cell performance and performance stability. When an oxygen electrode is under polarization, the oxygen partial pressure at the interface between the oxygen electrode and the electrolyte is lower than that of incoming oxidant. Under a high polarization, the environment at the aforementioned interface may lead to phase transformation of the oxygen electrode. The local oxygen partial pressure is determined by the transport properties at the interfaces. An addition of the electronic conduction in the interface layer, for instance using (Pr,Gd)-doped ceria to replace Gd-doped ceria, results in improved cell performance and performance stability, while the phase transformation is significantly suppressed. This work provides a fundamental understanding of the origin for phase transformation in oxygen electrodes during operation and use this knowledge to develop a high-performance electrode that exhibits improved performance stability.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Efficient machine-learning model for fast assessment of elastic properties of high-entropy alloys

We combined descriptor-based analytical models for stiffness-matrix and elastic-moduli with mean-field methods to accelerate assessment of technologically useful properties of high-entropy alloys, such as strength and ductility. Model training for elastic properties uses Sure-Independence Screening (SIS) and Sparsifying Operator (SO) method yielding an optimal analytical model, constructed with meaningful atomic features to predict target properties. Computationally inexpensive analytical descriptors were trained using a database of elastic properties determined from density functional theory for binary and ternary subsets of Nb-Mo-Ta-W-V refractory alloys. The optimal Elastic-SISSO models, extracted from an exponentially large feature space, give an extremely accurate prediction of target properties, similar to or better than other models, with some verified from existing experiments. Here we also show that electronegativity variance and elastic-moduli can directly predict trends in ductility and yield strength of refractory HEAs, and reveals promising alloy concentration regions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Field-assisted suppression of coke in the methane steam reforming reaction

Electric fields were applied to an unsupported nickel catalyst during methane steam reforming with a low steam-to-methane ratio of 2:1. As expected, catalysts operated in standard mode (without an applied field) experienced high levels of coking. Coke formation was almost completely suppressed by the application of a positive electric field. This is explained using density functional theory computations which indicate that positive electric fields (1) increase the affinity for water adsorption on the surface, thereby increasing the amount of oxygen available to remove carbon, (2) stabilize individual carbon fragments, thus reducing their tendency to polymerize, and (3) increase the oxygen storage capabilities of the catalyst by inducing sub-surface oxygen formation and bulk oxidation.

Gray, Jake T.↗
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 OSTI2024

Temperature-dependent iron motion in extremophile rubredoxins – no need for ‘corresponding states’

Extremophile organisms are known that can metabolize at temperatures down to – 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins – rubredoxins – from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of ‘corresponding states’ posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although ‘corresponding states’ would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗
DOE OSTI2026

Generalized Hall conductivities in local commuting projector models: Generalized symmetries and protected surface modes

Hall conductivities are important characterizations of phases of matter. It is known that nonzero Hall conductivities are difficult to realize in local commuting projector lattice models due to no-go theorems in (2+1)⁢D. In this work we construct local commuting projector models in (2+1)⁢D and (3+1)⁢D with nonzero generalized Hall conductivities for ordinary and higher-form continuous symmetries on tensor product Hilbert space of finite local dimension. The model is given by a standard ℤ 𝑁 toric code, but the symmetries do not admit expression in terms of on-site charge operators. The symmetry do not have local charges or currents on the lattice in the absence of boundaries, but there is still a notion of Hall conductivities that coincide with the continuum field theories. We construct protected gapless boundaries of the lattice models using modified Villain formalism. The generalized Hall conductivities are computed by surface currents as well as bulk flux insertion and many-body Chern number.

Anomalies↗
DOE OSTI2024

Scale and conformal invariance on (A)dS spacetimes

We examine the question of scale versus conformal invariance on maximally symmetric curved backgrounds and study general two-derivative conformally invariant free theories of vectors and tensors. For spacetime dimension D > 4 , these conformal theories can be diagonalized into standard massive fields in which unbroken conformal symmetry nontrivially mixes components of different spins. In D = 4 , the tensor case becomes a conformal theory mixing a partially massless spin-2 field with a massless spin-1 field. For massless linearized gravity in D = 4 , we confirm through direct calculation that correlation functions of gauge-invariant operators take the conformally invariant form, despite the absence of standard conformal symmetry at the level of the action. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2022

POLCA8 - modelling of cross section variations inside hexagonal assemblies

This paper presents the POLCA8 approach for modelling non-constant cross section distributions inside hexagonal fuel assemblies. The multigroup diffusion equation is modified to account for intranodal cross section variations. The obtained equation is solved in a node-wise manner based on the Fourier expansion method. As a result of varying cross sections, the solution includes a particular part additionally to the homogeneous one. A method for obtaining the particular solution is derived. Numerical tests on a VVER-1000 core are presented showing the impact of cross-section variations to some key parameters for reactor operation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Demonstrate Capability of NEAMS Tools to Generate Reactor Kinetics Parameters for Pebble-Bed HTGRs Transient Modeling

The system analysis of anticipated operating occurrences and design basis accident for pebblebed reactor systems requires knowledge of the neutron kinetic parameters. These parameters take into account various types of feedback from the pebble-bed core and are typically integrated into system analysis tools such as the System Analysis Module (SAM) in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tool package. These parameters are typically generated by a higher fidelity full-core coupled neutronics/thermal fluids analysis, such as using Griffin and Pronghorn, in 2D or 3D. The delayed neutron fractions (beta) and the neutron lifetimes (lambda) are typically generated by applying the adjoint solution of the neutron flux. Temperature reactivity coefficients are the most important feedback for all anticipated operating occurrences and design basis accident, then depending on the transient other effects might be important as rector control system worth, s-curves, and Xenon generation. This work package will establish a methodology and process for calculating the reactor kinetics parameters for high-temperature gas-cooled reactor and provide them in a suitable form for system level analyses of typical anticipated operating occurrences and design basis accident.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Entanglement and the double copy

We construct entangled states of gluons that scatter exactly as if they were gravitons. Operationally, these objects implement the double copy at the level of the wave function. Our analysis begins with a general ansatz for a wave function characterizing gluons in two copies of SU(N ) gauge theory. Given relatively minimal assumptions following from permutation invariance and dimensional analysis, the three- and four-particle wave functions generate scattering amplitudes that automatically coincide exactly with gravity, modulo normalization. For five-particle scattering the match is not automatic but imposing certain known selection rules on the amplitude is sufficient to uniquely reproduce gravity. The resulting amplitudes exhibit a color-dressed and permutation-invariant form of the usual double copy relations. We compute the entanglement entropy between the two gauge theory copies and learn that these states are maximally-entangled at large N. Moreover, this approach extends immediately to effective field theories, where Born-Infeld photons and Galileons can be similarly recast as entangled gluons and pions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Enhanced multi-carbon alcohol electroproduction from CO via modulated hydrogen adsorption

Multi-carbon alcohols such as ethanol are valued as fuels in view of their high energy density and ready transport. Unfortunately, the selectivity toward alcohols in CO 2 /CO electroreduction is diminished by ethylene production, especially when operating at high current densities (>100 mA cm –2 ). Here we report a metal doping approach to tune the adsorption of hydrogen at the copper surface and thereby promote alcohol production. Using density functional theory calculations, we screen a suite of transition metal dopants and find that incorporating Pd in Cu moderates hydrogen adsorption and assists the hydrogenation of C 2 intermediates, providing a means to favour alcohol production and suppress ethylene. We synthesize a Pd-doped Cu catalyst that achieves a Faradaic efficiency of 40% toward alcohols and a partial current density of 277 mA cm –2 from CO electroreduction. The activity exceeds that of prior reports by a factor of 2.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Entanglement phase transitions in random stabilizer tensor networks

Here we explore a class of random tensor network models with “stabilizer” local tensors which we name random stabilizer tensor networks (RSTNs). For RSTNs defined on a two-dimensional square lattice, we perform extensive numerical studies of entanglement phase transitions between volume-law and area-law entangled phases of the one-dimensional boundary states. These transitions occur when either (a) the bond dimension D of the constituent tensors is varied or (b) the tensor network is subject to random breaking of bulk bonds, implemented by forced measurements. In the absence of broken bonds, we find that the RSTN supports a volume-law entangled boundary state with bond dimension D ≥ 3 where D is a prime number, and an area-law entangled boundary state for D = 2. Upon breaking bonds at random in the bulk with probability p, there exists a critical measurement rate p c for each D ≥ 3 above which the boundary state becomes area-law entangled. To explore the conformal invariance at these entanglement transitions for different prime D, we consider tensor networks on a finite rectangular geometry with a variety of boundary conditions, and extract universal operator scaling dimensions via extensive numerical calculations of the entanglement entropy, mutual information, and mutual negativity at their respective critical points. Our results at large D approach known universal data of percolation conformal field theory, while showing clear discrepancies at smaller D, suggesting a distinct entanglement transition universality class for each prime D. We further study universal entanglement properties in the volume-law phase and demonstrate quantitative agreement with the recently proposed description in terms of a directed polymer in a random environment.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Precision Reconstruction of Rational Conformal Field Theory from Exact Fixed-Point Tensor Network

The novel concept of entanglement renormalization and its corresponding tensor network renormalization technique have been highly successful in developing a controlled real-space renormalization group (RG) scheme. Numerically approximate fixed-point (FP) tensors are widely used to extract the conformal data of the underlying conformal field theory (CFT) describing critical phenomena. In this paper, we present an explicit analytical construction of the FP tensor for 2D rational CFT. We define it as a correlation function between the “boundary-changing operators” (BCO) on triangles. Our construction fully captures all the real-space RG conditions. We also provide concrete examples, such as Ising, Yang-Lee, and tricritical Ising models, to compute the scaling dimensions explicitly based on the corresponding FP tensor. The BCO descendants turn out to be an optimal basis such that truncation in bond dimensions naturally produces comparable accuracies with the leading existing FP algorithms. Interestingly, our construction of FP tensors is closely related to a strange correlator, where the holographic picture naturally emerges. Our results also open a new door toward understanding CFT in higher dimensions. Published by the American Physical Society 2025

Cheng, Gong (ORCID:0009000891587404)↗
DOE OSTI2025

Nonresonant two-photon x-ray absorption in Cu

We present a real-space Green's function theory and calculations of two-photon x-ray absorption (TPA). Our focus is on nonresonant 𝐾-shell TPA in metallic Cu, which has been observed experimentally at intense x-ray free electron laser (XFEL) sources. The theory is based on an independent particle Green's function treatment of the Kramers-Heisenberg equation and an approximation for the sum over nonresonant intermediate states in terms of a static quadrupole transition operator. XFEL effects are modeled by a partially depleted 𝑑 band. This approach is shown to give results for 𝐾-shell TPA in quantitative agreement with XFEL experiment and with a Bethe-Salpeter equation approach. Furthermore, we also briefly discuss many-body corrections and TPA sum rules.

Approximation methods for many-body systems↗