Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled system”

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

Hybrid Energy System to H2 to Green Steel/Ammonia

Green hydrogen is a vital pathway to decarbonize industrial sectors responsible for the largest portion of greenhouse gas emissions, including ammonia/fertilizer and steel. It takes approximately a GW-scale renewables hybrid-H2 system to fully convert one average U.S. steel plant to a green steel plant, with an annual production capacity of 100 million metric tons per year. The optimized cost-effective operation of renewables-H2 plants are highly dependent on location and the specific end use designs. This project focused on the analysis of GW-scale off-grid, tightly-coupled, co-located systems to reduce overall cost compared with steam methane reform with and without carbon capture and grid connected systems (see Figure 1).

ammonia↗
DOE OSTI2023

Superconductivity in the Li-B-C system at 100 GPa

Layer Li-B-C compounds have been shown to have feasible superconductivity. Using the adaptive genetic algorithm, we predict the structures of the Li-B-C system at 100 GPa. Here we identify several low-enthalpy metallic phases with stoichiometries of Li B 2 C, Li B 3 C, Li 2 BC 2 , Li 3 B 2 C 3 , Li 3 BC, and Li 5 BC. Using a fast screening method of electron-phonon interaction, we find that Li B 3 C is a promising candidate for superconductivity. The consecutive calculations using the full Brillouin zone confirm the existence of the strong electron-phonon coupling (EPC) in this system. The anharmonic B-C phonon modes near the zone center provide the major contribution to the EPC. The EPC constant is 1.40, and the estimated critical temperature is 22 K. In this paper, we indicate that superconductivity can also happen without a layered structural motif in the Li-B-C system. We also demonstrate an effective strategy for crystal structure prediction of superconducting materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Polynomial Scaling Localized Active Space Unitary Selective Coupled Cluster Singles and Doubles

We present a polynomial-scaling algorithm for the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. In this approach, cluster excitations are selected based on a threshold ϵ determined by the absolute gradients of the LAS-UCCSD energy with respect to cluster amplitudes. Using the generalized Wick’s theorem for multireference wave functions, we derive the gradient expression as a polynomial function of one-, two-, and three-body reduced density matrices and 1- and 2-electron integrals, valid for any multireference wave function. The resulting gradient implementation exhibits a memory scaling of 𝒪(N 6 ), with N spin orbitals in the combined active space of all fragments. The variational quantum eigensolver is used to optimize the selected cluster excitations on a quantum simulator. Furthermore, by plotting the energy error, defined as the difference between the LAS-USCCSD and corresponding CASCI energies, against the inverse cluster amplitude selection threshold (ϵ –1 ) for polyene chains containing 2 to 5 π-bond units, we establish a relationship between the energy error and the threshold. To further validate the accuracy of LAS-USCCSD, we computed the cis–trans isomerization energy of stilbene (a 20-qubit system) and the magnetic coupling constant of the tris-hydroxo-bridged chromium dimer [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ (evaluated as both 12- and 20-qubit systems) using the Qiskit-Qulacs simulator. Assessing such examples is important to determine the practical feasibility of quantum simulations for chemically realistic systems. Toward this goal, with the LAS-USCCSD algorithm we estimated the quantum resources required for simulating an active space of (30e,22o) in [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ , a size that remains beyond the reach of current quantum simulators for accurate treatment.

Algorithms↗
DOE OSTI2025

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Barriers to carriers: faults and recombination in non-stoichiometric perovskite scintillators

Abstract Tuning the efficiency and speed of charge carrier recombination in inorganic scintillators can potentially improve their performance in diverse applications. Recent work suggests that this maybe be achieved via a two-phase scintillator AB that naturally phase separates into A -rich and B -rich domains. In addition, a favorable electronic structure and band-edge alignment such that the charge carriers are confined or are thermodynamically driven to preferentially accumulate in one of the two domains, might lead to an improved radiative recombination rate. Here, we use density functional theory computations and ab initio molecular dynamics (AIMD), including non-adiabatic molecular dynamics (NAMD) simulations, to examine an alternative phase structure and its potential impact on recombination. Using a model perovskite SrTiO $$_3$$ 3 system with one-, two- and three-dimensional Ruddlesden–Popper (RP) phases, we demonstrate that RP faults induce band structure changes in the material that can act as barriers to carrier transport. Our AIMD/NAMD simulations indicate competing effects of a lower mean free path (potentially enhancing the desired radiative recombination and overall scintillating efficiency) and faster non-radiative recombination (undesired) due to enhanced electron–phonon coupling in the faulted system. Full exploitation of such a rational design approach would require tuning of the effective scintillation efficiency by varying the perovskite chemistry using appropriate arrangements of RP faults in the bulk material. Finally, other effects, such as the tendency of point defects to segregate at the interface, that might affect the overall performance, are briefly discussed. We expect the basic results found here to apply to other nanostructured scintillators. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dynamics and transport in the boundary-driven dissipative Klein-Gordon chain

Motivated by experiments on chains of superconducting qubits, we consider the dynamics of a classical Klein-Gordon chain coupled to coherent driving and subject to dissipation solely at its boundaries. As the strength of the boundary driving is increased, this minimal classical model recovers the main features of the "dissipative phase transition" seen experimentally. Between the transmitting and non-transmitting regimes on either side of this transition (which support ballistic and diffusive energy transport respectively), we observe additional dynamical regimes of interest. These include a regime of superdiffusive energy transport at weaker driving strengths, together with a "dissipative breather crystal" regime at stronger driving strengths, which is characterized by emergent translation symmetry, ballistic energy transport, and coherent oscillations of a nonlinear normal mode. Here, we propose a non-local Lyapunov exponent as an experimentally measurable diagnostic of many-body chaos in this system, and more generally in open systems that are only coupled to an environment at their boundaries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Positronium: an illustration of nonperturbative renormalization in a basis light-front approach

We calculate the mass spectrum and the structure of the positronium system at a strong coupling in a basis light-front approach. We start from the light-front QED Hamiltonian and retain one dynamical photon in our basis. We perform the fermion mass renormalization associated with the nonperturbative fermion self-energy correction. We present the resulting mass spectrum and wave functions for the selected low-lying states. Next, we apply this approach to QCD and calculate the heavy meson system with one dynamical gluon retained. We illustrate the obtained mass spectrum and wave functions for the selected low-lying states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Enriched immersed boundary method (EIBM) for interface-coupled multi-physics and applications to convective conjugate heat transfer

It is challenging to develop numerical methods that simultaneously maintain accuracy of resolving boundary conditions and mesh flexibility to handle the interface in interface-coupled multi-physics problems involving large property discontinuity. On the one hand, boundary-fitted methods possess high accuracy in capturing the interfacial phenomenon but involve complicated volumetric mesh generation, mesh-motion, and even re-meshing procedures. On the other hand, immersed boundary methods (IBM) provide mesh flexibility but sometimes suffer from inferior interface representations, leading to poor enforcement of boundary conditions. This paper presents an enriched immersed boundary method (EIBM) to overcome this challenge and demonstrates its efficacy in conjugate heat transfer, a representative example in interface-coupled multi-physics systems that have implications for many industrial processes. Further, the core technique of the method is to enhance IBM’s accuracy of the fluid–solid interface by enriching the degrees of freedom of the cut elements to enforce temperature and flux compatibilities and resolve all the physical unknowns on the background mesh to simplify the volumetric mesh generation. We implement the EIBM under the framework of a variational multiscale formulation for coupled Navier–Stokes and thermodynamics equations. The enriched DoFs enable better enforcement of temperature and flux compatibilities with large conductivity ratios across the fluid–solid interface. At the same time, the immersed nature of the proposed method still attains mesh flexibility. The EIBM’s accuracy is thoroughly evaluated through a set of examples, ranging from benchmark problems with analytical solutions to real-world cooling processes of a moving metallic structure with complex geometry.

42 ENGINEERING↗
DOE OSTI2024

Innovative SCR Materials and System for Low Temperature - CRADA 350 (Abstract)

The aim of this CRADA is focused on providing a new enabling SCR catalyst system that will function at very high efficiency to attain the most demanding emissions regulations and thereby facilitate the market introduction of advanced powertrains that will support domestic energy independence and security. Future powertrains, that will be significantly more efficient than currently available technologies, will be needed by automotive manufacturers to meet rapidly increasing CAFE and GHG standards. However, these powertrains cannot enter into the US light duty vehicle market unless they are coupled with an aftertreatment system that will sufficiently remediate tailpipe emissions to meet EPA Tier III and California SULEV emissions standards. The low temperature exhaust associated with these powertrains is especially challenging for any current aftertreatment technology to meet these standards. The key focus of this CRADA is to further develop newly invented materials for the selective catalytic reduction (SCR) of NOx by ammonia (NH3) that show promise for significantly reducing ‘light-off’ temperatures compared to current commercial catalysts. Specifically, the goal of the proposed work is to achieve ‘light-off’ of NH3 SCR at 150 ºC in order to realize conversion efficiencies of 90% at these low temperatures. This will enable deployment of lean combustion powertrains with significantly increased fuel efficiencies but lower exhaust temperatures. To accomplish this overall goal, it will be essential to also identify an appropriate NH3 supply strategy for the SCR aftertreatment device that can controllably deliver NH3 at these low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2024

Integrating subsystem embedding subalgebras and coupled cluster Green’s function: a theoretical foundation for quantum embedding in excitation manifold

Here, in this study, we introduce a novel approach to coupled-cluster Green's function (CCGF) embedding by seamlessly integrating conventional CCGF theory with the state-of-the-art sub-system embedding sub-algebras coupled cluster (SES-CC) formalism. This integration focuses primarily on delineating the characteristics of the sub-system and the corresponding segments of the Green's function, defined explicitly by active orbitals. Crucially, our work involves the adaptation of the SES-CC paradigm, addressing the left eigenvalue problem through a distinct form of Hamiltonian similarity transformation. This advancement not only facilitates a comprehensive representation of the interaction between the embedded sub-system and its surrounding environment but also paves the way for the quantum mechanical description of multiple embedded domains, particularly by employing the emergent quantum flow algorithms. Our theoretical underpinnings further set the stage for a generalization to multiple embedded sub-systems. This expansion holds significant promise for the exploration and application of non-equilibrium quantum systems, enhancing the understanding of system–environment interactions. In doing so, the research underscores the potential of SES-CC embedding within the realm of quantum computations and multi-scale simulations, promising a good balance between accuracy and computational efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

gdess: A framework for evaluating simulated atmospheric CO 2 in Earth System Models

Atmospheric carbon dioxide (CO 2 ) plays a key role in the global carbon cycle and global warming. Climate-carbon feedbacks are often studied and estimated using Earth System Models (ESMs), which couple together multiple model components—including the atmosphere, ocean, terrestrial biosphere, and cryosphere—to jointly simulate mass and energy exchanges within and between these components. Despite tremendous advances, model intercomparisons and benchmarking are aspects of ESMs that warrant further improvement (Fer et al., 2021; Smith et al., 2014). Such benchmarking is critical because comparing the value of state variables in these simulations against observed values provides evidence for appropriately refining model components; moreover, researchers can learn much about Earth system dynamics in the process (Randall et al., 2019). We introduce `gdess` (a.k.a., Greenhouse gas Diagnostics for Earth System Simulations), which parses observational datasets and ESM simulation output, combines them to be in a consistent structure, computes statistical metrics, and generates diagnostic visualizations. In its current incarnation, `gdess` facilitates evaluating a model's ability to reproduce observed temporal and spatial variations of atmospheric CO 2 . The diagnostics implemented modularly in `gdess` support more rapid assessment and improvement of model-simulated global CO 2 sources and sinks associated with land and ocean ecosystem processes. We intend for this set of automated diagnostics to form an extensible, open source framework for future comparisons of simulated and observed concentrations of various greenhouse gases across Earth system models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗
DOE OSTI2022

Surface fluctuating hydrodynamics methods for the drift-diffusion dynamics of particles and microstructures within curved fluid interfaces

Here we introduce fluctuating hydrodynamics approaches on surfaces for capturing the drift-diffusion dynamics of particles and microstructures immersed within curved fluid interfaces of spherical shape. We take into account the interfacial hydrodynamic coupling, traction coupling with the surrounding bulk fluid, and thermal fluctuations. For fluid-structure interactions, we introduce Immersed Boundary Methods (IBM) and related Stochastic Eulerian-Lagrangian Methods (SELM) for curved surfaces. We use these approaches to investigate the statistics of surface fluctuating hydrodynamics and microstructures. For velocity autocorrelations, we find characteristic power-law scalings $τ^{-1}, τ^{-2}$, and plateaus can emerge. This depends on the physical regime associated with the geometry, surface viscosity, and bulk viscosity. This differs from the characteristic $τ^{-3/2}$ scaling for bulk three dimensional fluids. We develop theory explaining these observed power-laws associated with time-scales for dissipation within the fluid interface and coupling to the surrounding fluid. We then use our introduced methods to investigate a few example systems and roles of hydrodynamic coupling and thermal fluctuations including for the kinetics of passive particles and active microswimmers in curved fluid interfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Nonlinear nanomechanical resonators approaching the quantum ground state

It is an open question whether mechanical resonators can be made nonlinear with vibrations approaching the quantum ground state. This requires the engineering of a mechanical nonlinearity far beyond what has been realized so far. Here we discover a mechanism to boost the Duffing nonlinearity by coupling the vibrations of a nanotube resonator to single-electron tunnelling and by operating the system in the ultrastrong-coupling regime. We find that thermal vibrations become highly nonlinear when lowering the temperature. The average vibration amplitude at the lowest temperature is 13 times the zero-point motion, with approximately 42% of the thermal energy stored in the anharmonic part of the potential. Our work may enable the realization of mechanical Schrödinger cat states, mechanical qubits and quantum simulators emulating the electron–phonon coupling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microscopic Imprints of Learned Solutions in Tunable Networks

In physical networks trained using supervised learning, physical parameters are adjusted to produce desired responses to inputs. An example is an electrical contrastive local learning network of nodes connected by edges that adjust their conductances during training. When an edge conductance changes, it upsets the current balance of every node. In response, physics adjusts the node voltages to minimize the dissipated power. Learning in these systems is therefore a coupled double-optimization process, in which the network descends both a cost landscape in the high-dimensional space of edge conductances and a physical landscape—the power dissipation—in the high-dimensional space of node voltages. Because of this coupling, the physical landscape of a trained network contains information about the learned task. Here, we derive a structure-function relation for trained tunable networks and demonstrate that all the physical information relevant to the trained input-output relation can be captured by a tuning susceptibility, an experimentally measurable quantity. We supplement our theoretical results with simulations to show that the tuning susceptibility is correlated with functional importance and that we can extract physical insight into how the system performs the task from the conductances of highly susceptible edges. Our analysis is general and can be applied directly to mechanical networks, such as networks trained for protein-inspired function such as allostery.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Quest for Stable Potassium‐Ion Battery Chemistry

Abstract Potassium‐ion batteries (KIBs) have attracted wide interest for energy storage because of the abundance of the electrode materials involved; however, their electrochemical performances are far behind what can be achieved from lithium‐ion batteries (LIBs) or sodium‐ion batteries (SIBs). Herein, key promising electrode and electrolyte materials for potassium‐ion batteries are identified, the coupled electrochemical reactions in the cell are investigated, and the compatibility between different materials is demonstrated to play the most important role. K 2 Mn[Fe(CN) 6 ] cathode can deliver a high capacity of ≈125 mAh g −1 and exceptional cycling stability over 61 000 cycles (≈9 months) if the side reactions from the anode can be prevented. Graphite is a good anode material but is subjected to degradation in traditional carbonate electrolytes. New concentrated electrolytes are developed and evaluated. A stable KIB system is demonstrated by coupling a stable K 2 Mn[Fe(CN) 6 ] cathode, a prepotassiated graphite anode with a concentrated electrolyte to achieve a high energy density of ≈260 Wh kg −1 (based on the active mass of cathode and anode) and good cycling of over 1000 cycles.

Wu, Xianyong↗