Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling method”

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

Elucidating the impact of the ionomer equivalent weight on a platinum group metal-free PEMFC cathode via oxygen limiting current

Leveraging the interactions between ionomer and catalyst can increase the performance of proton exchange membrane fuel cells. The impacts of the equivalent weight (EW) of perfluorosulfonic acid–based ionomers on the platinum group metal-free electrode structure and fuel cell performance have not been fully explored. Four membrane electrode assemblies (MEAs) were prepared by using a commercial Fe–N–C catalyst, two perfluorosulfonic acid ionomers with different EWs, that is, Aquivion 720 (A720) and Nafion 1100 (N1100), and two ionomer-to-catalyst (I/C) ratios. The four MEAs were characterized to understand the impact of the ionomer EW and content on the capacitance, proton conductivity, and mass transport on the cathode. The mass transport resistance was measured for the first time using a new oxygen reduction reaction limiting current method enabling to couple the effects of oxygen diffusion with liquid water generation. Low EW ionomer combined with a moderate I/C results in improved performance due to its enhanced proton conductivity. However, when used at high I/C, it can cause severe water flooding at high current density due to the enhanced liquid water uptake, especially at high relative humidity, resulting in lower catalyst utilization and higher mass transport resistance.

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

Towards a nonperturbative construction of the S-matrix

We present a nonperturbative recipe for directly computing the S-matrix in strongly-coupled QFTs. The method makes use of spectral data obtained in a Hamiltonian framework and can be applied to a wide range of theories, including potentially QCD. We demonstrate the utility of this prescription in the specific example of the 2+1d O(N) model at large N, using energy eigenstates computed with Hamiltonian truncation to reproduce the full 2 → 2 scattering amplitude for arbitrary (complex) center-of-mass energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Microstructure prediction for Ti-22Al-25Nb in laser powder bed fusion

This work presents a physics-informed framework for predicting solidification morphology and defect susceptibility in additively manufactured Ti–22Al–25Nb across a broad processing space. The framework integrates solidification microstructure selection (SMS) analysis with a single-track defect-based printability map to establish a unified methodology linking processing parameters to both interfacial morphology and manufacturability. Thermal gradients G and solidification rates R are first computed using the Thermo-Calc Additive Manufacturing (TC-AM) module, a finite-interface-dissipation (FID) phase-field (PF) model coupled with CALPHAD method is then employed to systematically distinguish planar and dendritic regimes as functions of $G$ and $R$. By superimposing the printability map onto the morphology projections, a comprehensive process–structure framework is obtained. Across most processing conditions, the predicted microstructure is predominantly dendritic, while planar growth emerges only under selected laser power $P$ and scan speed $v$ combinations. In addition to morphology classification, the framework quantifies the dendritic area fraction and introduces a width-based morphology descriptor to characterize the spatial extent of planar/dendritic regions within the melt pool. It provides mechanistic insight into the interplay between solidification physics and defect formation, offering practical guidance for parameter selection and microstructural control in Ti–22Al–25Nb additive manufacturing (AM).

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multi–length scale characterization of point defects in thermally oxidized, proton irradiated iron oxides

A key for the success of safe nuclear power generation system is to consider structural materials that are economical, meet mechanical property needs, possess good corrosion resistance, and are radiation tolerant. Nevertheless, fundamental insights that elucidate the details of radiation damage on materials corrosion performance are lacking. This includes the behavior of surface oxides which often regulate corrosion. For example, it is unclear how non-equilibrium point defects, oxide structure, mass transport in oxides, and subsequent oxidation behavior are altered by the radiation. Here, in this work, some of the effects of proton irradiation on the attributes of point defects, iron oxide microstructures, and the physical nature of the oxidation product were correlated with corrosion behavior. Iron oxides, fabricated by thermal oxidation in air at 400°C and 800°C for 1 h, were subjected to 200 keV, 0.03 dpa (displacements per atom) of proton irradiation, and subjected to corrosion reactivity assessment using AC and DC electrochemical methods. Experimental methods that target materials information at different length scales, such as positron annihilation spectroscopy (atomistic), transmission electron microscopy (mesoscopic), and electrochemical methods (macroscopic), were coupled to shed light on the impact of radiation-induced defect modifications and structural changes in oxides on corrosion reactivity which ultimately affects durability in harsh environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via 13 C-metabolic flux analysis ( 13 C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated 13 C labeling data were used with 13 C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Theory of Plasmonic Hot-Carrier Generation and Relaxation

Hot-carrier (HC) generation from (localized) surface plasmon decay has recently attracted much attention due to its promising applications in physical, chemical, materials, and energy science. However, the detailed mechanisms of plasmonic HC generation, relaxation, and trapping are less studied. In this work, we developed and applied a quantum-mechanical model and coupled master equation method to study the generation of HCs from plasmon decay and their following relaxation processes with different mechanisms treated on equal footing. First, a quantum-mechanical model for HC generation is developed. Its connection to existing semiclassical models and time-dependent density functional theory (TDDFT) is discussed. Second, the relaxation and lifetimes of HCs are investigated in the presence of electron–electron and electron–phonon interactions. A GW-like approximation is introduced to account for the electron–electron scattering. The numerical simulations on the Jellium nanoparticles with a size up to 1.6 nm demonstrate the electron–electron scattering and electron–phonon scattering dominate different time scale in the relaxation dynamics. We also generalize the model to study the extraction of HCs to attached molecules. The quantum yield of extracting HCs for other applications is found to be size-dependent. In general, the smaller size of NP improves the quantum yield, which is in agreement with recent experimental measurements. Even though we demonstrate this newly developed theoretical formalism with Jellium model, the theory applies to any other atomistic models.

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

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Spiers Memorial Lecture: activating metal sites for biological electron transfer

Metal sites in biology often exhibit unique spectroscopic features that reflect novel geometric and electronic structures imposed by the protein that are key to reactivity. The blue copper active site involved in long range, rapid biological electron transfer is a classic example. This review presents an overview of both traditional and synchrotron based spectroscopic methods and their coupling to electronic structure calculations to understand the unique features of the blue copper active site, their contributions to function and the role of the protein in determining the geometric and electronic structure of the active site (called the “entatic state”). The relation of this active site to other biological electron transfer sites is further developed. Specifically, ultrafast XFEL spectroscopy is used to evaluate the methionine–S–Fe bond in cytochrome c, and its entatic control by the protein in determining function (electron transfer vs. apoptosis).

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

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

Here, the exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron–ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron–nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

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

Binding and excitations in Si x H y molecular systems using quantum Monte Carlo

In this work, we present high-accuracy correlated calculations of small SixHy molecular systems in both the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within ≈0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between 1% and 3.5% with absolute values at most ≈0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations.

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

Self-consistent description of the halo nature of 31 Ne with continuum and pairing correlations

Here a relativistic structure model has previously been used to predict a halo structure for 31 Ne, consistent with halo signatures from measured reaction cross sections of Ne isotopes bombarding Carbon targets. However, previous attempts to calculate those cross sections with reaction models were missing contributions from resonances and pairing correlations in their structure input. Use a reaction model with our relativistic fully microscopic structure model input to predict these cross sections and momentum distributions and analyze for possible halo signatures. Structure input for exotic Ne isotopes were obtained via the analytical continuation of the coupling constant (ACCC) method based on the relativistic mean field (RMF) theory with Bardeen - Cooper - Schrieffer (BCS) pairing approximation, the RAB approach. Total reaction cross sections, one-neutron removal cross sections, and momentum distributions of breakup reaction products were calculated with a Glauber model using our relativistic structure input. Our predictions of total reaction and one-neutron removal cross sections of 31 Ne on a Carbon target were significantly enhanced compared with those of neighboring Neon isotopes, agreeing with measurements at 240 MeV/nucleon and consistent with a single neutron halo. Furthermore, our calculations of the inclusive longitudinal momentum distribution of the 30 Ne and valence neutron residues from the 31 Ne breakup reaction indicate a dilute density distribution in coordinate space, another halo signature. We give a full description of the halo nature of 31 Ne that includes a self-consistent use of pairing and continuum contributions that makes predictions consistent with reaction cross section measurements. This approach can be utilized to determine the halo structure of other exotic nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Origin of multifractality in solar wind turbulence: the role of current sheets

ABSTRACT In this work, a multifractal framework is proposed to investigate the effects of current sheets in solar wind turbulence. By using multifractal detrended fluctuation analysis coupled with surrogate methods and volatility, two solar wind magnetic field time series are investigated, one with current sheets and one without current sheets. Despite the lack of extreme-events intermittent bursts in the current sheet-free series, both series are shown to be strongly multifractal, although the current sheet-free series displays an almost linear behaviour for the scaling exponent of structure functions. Long-range correlations are shown to be the main source of multifractality for the series without current sheets, while a combination of heavy-tail distribution and non-linear correlations is responsible for multifractality in the series with current sheets. The multifractality in both time series is formally shown to be associated with an energy-cascade process using the p-model.

Gomes, Leonardo F.↗
DOE OSTI2023

Solving the nuclear pairing model with neural network quantum states

In this work, we present a variational Monte Carlo method that solves the nuclear many-body problem in the occupation number formalism exploiting an artificial neural network representation of the groundstate wave function. A memory-efficient version of the stochastic reconfiguration algorithm is developed to train the network by minimizing the expectation value of the Hamiltonian. We benchmark this approach against widely used nuclear many-body methods by solving a model used to describe pairing in nuclei for different types of interaction and different values of the interaction strength. Despite its polynomial computational cost, our method outperforms coupled-cluster and provides energies that are in excellent agreement with the numerically-exact full configuration interaction values.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Magnetic excitations from the hexagonal spin clusters in the 𝑆 = $\frac{1}{2}$ distorted honeycomb lattice antiferromagnet Cu 2 ⁢(pymca)⁢ 3 ⁢(ClO 4 )

Cu 2 ⁢(pymca) ⁢3 (ClO 4 ) (pymca: pyrimidine-2-carboxylate) consists of a slightly distorted honeycomb lattice of Cu 2+ spins, which shows no long-range magnetic order down to 0.6 K. A magnetization study revealed 1/3 and 2/3 plateau phases [A. Okutani et al., J. Phys. Soc. Jpn. 88, 013703 (2019)], which is not expected for regular honeycomb antiferromagnets. Inelastic neutron scattering experiments were performed using a powder sample to investigate the exchange interactions of this material. The spin excitations from the singlet ground state to the first three triplet states, predicted from the antiferromagnetic hexagonal spin cluster interacting with 3.9 meV, were observed. Using the exact diagonalization methods, the intercluster coupling was estimated from the excitation peak width to be about 20% of the intracluster interaction, which is consistent with the previously reported value. Finally, our exchange path model explains the anisotropic exchange interactions in the distorted honeycomb plane.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Grid Parameters and Voltage Estimation Approach Integrating Data-Driven Converter Model

With Measurements of grid voltage and current are essential for the optimal operation of the grid protection and control (P&C) systems. Grid parameters vary through time during the faults and especially in the converter interfaced resources (CIRs) rich power grid, and thus accurate estimation is critical to avoid the mis-operation of the P&C systems. In this paper, a moving horizon estimation (MHE) as an observer is devised and applied to estimate the grid line parameters and grid voltages for protection enhancement. Due to the proprietary and confidentiality of CIRs, the proposed approach uses the black-box model to represent their dynamics. Leveraging the easily accessible measurements of output current from the black-box model of CIR and voltage at the point of common coupling, the proposed method estimates the grid impedance and grid voltage during normal and faulty operating conditions. The performance shows that the optimization-based observer was able to closely observe the accurate states and parameters, which can be utilized by the P&C systems.

Subedi, Sunil↗
DOE OSTI2022

Verification of the LOTUS code with C5G7 benchmark

In this study, the verification of the neutron transport code LOTUS through the well-known C5G7 benchmark is performed. LOTUS uses the current coupled collision probability method with the expansion of the flux by orthogonal polynomials (CCCPO) for the solution of the neutron transport equation. The expansion of the flux by orthogonal polynomials allows one to avoid discretization of the calculation regions and significantly decreases the simulation time. The results of the LOTUS calculations are compared with the results of the reference OpenMC Monte Carlo simulations. The results of the comparison demonstrate almost perfect agreement with Monte-Carlo for the second order of the flux expansion. Even though the results obtained in the presented study agree well with the results of reference Monte Carlo calculations, further investigations are necessary for a better understanding of the stability and limitations of the code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗