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 361 records · Page 20

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

UMass 2-Body WEC Techno-Economic Assessment

The University of Massachusetts (UMass) is developing a 2-body wave energy converter (WEC) device that is converting mechanical power into electricity using a mechanical motion rectifier that allows the system to couple to a flywheel. UMass has completed numerical modeling, wave tank testing, and PTO sub-system testing and needed assistance in developing a techno-economic model to enable optimization of their topology, comparison to a generic heaving point absorber topology, and guide the next steps in their development efforts. The core objective was to develop a techno-economic approach and modeling tool that allows benchmarking of the two topologies across a wide range of scales to evaluate their respective competitiveness in different application spaces. This data includes the final report as well as a supporting spreadsheet containing the data produced for this report.

16 TIDAL AND WAVE POWER↗

HYBRID Modeling Validation and Verification Status Matrix

The HYBRID modeling repository is a premier resource for integrated energy systems modeling. HYBRID models have been developed since 2015 to describe the physical operation of tightly coupled thermal systems including power generators, thermal transport systems, thermal storage, thermal-to-electric conversion systems, and other thermal applications. Due to the increased size of the repository, a concise summary matrix of available models is desired. This matrix will consolidate not only the list of available models but also indicate original information sources, publications that have model examples, and level of validation and verification that exists for the models. Validation and verification (V&V) levels begin from simplified algebraic relationship and advance to dynamic data validation. Moving forward, as new models are contributed to HYBRID, their V&V level will be included, updating this matrix.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation and Verification of TEDS Facility HYBRID Modeling

The HYBRID modeling repository is an in-house developed library of models for selected integrated energy systems (IES) modelling. HYBRID models have been developed since 2015 to describe the physical operation of tightly coupled thermal systems including power generators, thermal transport systems, thermal storage, thermal-to-electric conversion systems, and other thermal applications. Here, validation and verification (V&V) capabilities are demonstrated using the Thermal Energy Distribution System (TEDS) at INL. Building upon prior work, the TEDS model has been updated and verified so that it better represents the installed system configuration and the operating control system. The model control system was changed to allow replication of actual experimental procedures. Experimental operations focusing primarily on thermocline tank performance were devised and performed. Several anomalies were found in the operation data of the experiment facility. V&V activities calibrating a selected input parameter are demonstrated on a single component as well as with a single parameter within the thermocline. Calibrating is then demonstrated on multiple components and a multi-parameter metric for the entire system. The validation methodology is successfully applied to validate the model with experimental data. It is also used to confirm a hypothesis behind one of the anomalies in experimental performance.

25 ENERGY STORAGE↗

Anisotropic superconductivity mediated by ferroelectric fluctuations in cubic systems with spin-orbit coupling

Motivated by the experimental observation that superconductivity in bulk doped SrTiO 3 is enhanced as a putative ferroelectric quantum critical point (FE-QCP) is approached, we study the pairing instability of a cubic system in which electrons exchange low-energy ferroelectric fluctuations. Instead of the gradient coupling to the lattice distortion associated with ferroelectricity, we consider direct coupling between the electrons and the bosonic ferroelectric field that appears in the presence of spin-orbit coupling. Working in the weak-coupling regime, we find that the pairing interaction is dominated by the soft transverse optical (TO) mode, resulting in a Tc enhancement upon approaching the FE-QCP. Focusing on even-parity states, we find that although the s-wave state always wins, states with higher Cooper-pair angular momentum become close competitors as the TO mode softens. We show that the cubic anisotropy of the FE fluctuations mixes the s-wave and g-wave states, resulting in a characteristic anisotropy of the gap function. Here, the gap anisotropy behaves nonmonotonically as the FE-QCP is approached: upon decreasing the TO mode frequency, the gap anisotropy first changes sign and then increases in magnitude. We discuss the possible applications of our results to the superconducting state of SrTiO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Topological flat bands in a kagome lattice multiorbital system

Flat bands and dispersive Dirac bands are known to coexist in the electronic bands in a two-dimensional kagome lattice. Including the relativistic spin-orbit coupling, such systems often exhibit nontrivial band topology, allowing for gapless edge modes between flat bands at several locations in the band structure, and dispersive bands or at the Dirac band crossing. Here, we theoretically demonstrate that a multiorbital system on a kagome lattice is a versatile platform to explore the interplay between nontrivial band topology and electronic interaction. Specifically, here we report that the multiorbital kagome model with the atomic spin–orbit coupling naturally supports topological bands characterized by nonzero Chern numbers $\mathscr{C}$, including a flat band with |$\mathscr{C}$|=1. When such a flat band is 1/3 filled, the non-local repulsive interactions induce a fractional Chern insulating state. We also discuss the possible realization of our findings in real kagome materials.

36 MATERIALS SCIENCE↗

Coupled topological flat and wide bands: Quasiparticle formation and destruction

Flat bands amplify correlation effects and are of extensive current interest. They provide a platform to explore both topology in correlated settings and correlation physics enriched by topology. Recent experiments in correlated kagome metals have found evidence for strange-metal behavior. A major theoretical challenge is to study the effect of local Coulomb repulsion when the band topology obstructs a real-space description. In a variant to the kagome lattice, we identify an orbital-selective Mott transition in any system of coupled topological flat and wide bands. This was made possible by the construction of exponentially localized and Kramers-doublet Wannier functions, which, in turn, leads to an effective Kondo-lattice description. Our findings show how quasiparticles are formed in such coupled topological flat-wide band systems and, equally important, how they are destroyed. Our work provides a conceptual framework for the understanding of the existing and emerging strange-metal properties in kagome metals and beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

CNRS Molten Salt Reactor Benchmark Analysis Using Griffin-Pronghorn Coupled Multi-physics Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time, since the fuel salt is circulating through the whole primary loop. This movement of fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and the corresponding redistribution of the DNPs in the active region of the core. To capture this effect, neutronics and thermal-hydraulics (T/H) codes need to be able to handle the movement of DNPs and their decay [1, 2]. Idaho National Laboratory (INL) is actively working on developing the neutronics and T/H codes to model the MSRs with flowing fuel. As one of the efforts, the neutronics code Griffin [3] and the T/H code Pronghorn [4] that are built upon the MOOSE framework [5] have been extended to handle the flowing fuel involving the drift of DNPs, and this capability is being tested and verified with MSR benchmark problems including the CNRS benchmark problem [6]. Additionally, a multi-physics analysis of the molten salt reactor experiment (MSRE) is being conducted to verify and validate the codes against available experimental data [7]. This work presents the analysis of the CNRS benchmark problem using the Griffin-Pronghorn coupled multi-physics code system of INL. Both steady-state and transient problems are analyzed, and the results are compared with those of other participants presented in [6]. The following section provides a description of the benchmark, along with input parameters and observables for each step. Then, the results for each step are presented in comparison with other participants’ results, followed by summary and conclusion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ORMATEX

The Oak Ridge Matrix Exponential (ORMATEX) software library contains methods to compute the matrix exponential and the action of the matrix exponential on a vector. Additionally, this package contains the related methods for the phi-functions which commonly appear in a wide class of exponential time integration methods. Krylov methods are provided to evaluate the matrix exponential-vector and phi-vector products for cases where the matrix is large and sparse. Utilizing these methods, ORMATEX implements performant exponential integrators for large systems of coupled ordinary differential equations (ODEs). The exponential time integration routines in ORMATEX are particularly suitable to large, stiff systems of equations. These routines may be utilized as a competitive alternative to classical implicit and explicit time integration schemes for certain classes of differential equations where the problem stiffness can be predominately explained by the linear terms.

Gurecky, William [Oak Ridge National Laboratory (O↗

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

Sub-system self-consistency in coupled cluster theory

Here, in this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dzyaloshinskii-Moriya interaction induced magnetoelectric coupling in a tetrahedral molecular spin-frustrated system

We have investigated magnetoelectric (ME) coupling in the single-molecule magnet Mn 4 Te 4 (PEt 3 ) 4 with tetrahedral spin frustration. Our density functional studies found that an electric dipole moment can emerge with various noncollinear spin orderings. The forms of spin-dependent dipole are determined and consistent with that in noncentrosymmetric magnets driven by the Dzyaloshinskii-Moriya interaction. Writing a parameterized spin Hamiltonian, after solving for eigenvalues and eigenstates, we quantified the ME coupling by calculating the thermal average of the electric and magnetic susceptibilities, which can be influenced by external magnetic and electric fields, respectively. The quadratic relations are expected to be observable in experiments.

36 MATERIALS SCIENCE↗

On the effects of quadrupolar relaxation in Earth’s field NMR spectra

There is growing interest in using low-field magnetic resonance experiments for routine chemical characterization. Earth’s field NMR is one such technique that can garner structural information and enable sample differentiation with low cost and highly portable designs. The resulting NMR spectra are primarily influenced by $J$-couplings, resulting in so-called $J$-coupled spectra (JCS). Many small molecules include atoms with NMR-active nuclei that are quadrupolar either at natural abundance or are often isotopically enriched (e.g., 2 H, 6 Li, 11 B, 14 N, 17 O, etc.) where the effects of quadrupolar J-couplings and relaxation on JCS of strongly- and weakly-coupled spin systems have not been explored to date. Herein, using a set of seven fluoropyridine samples with unique substitution and $J$-couplings, we demonstrate that the 14 N relaxation rates can induce drastic line-broadening in the JCS. This includes a previously unexplored unique line broadening mechanism enabled by strongly coupled spins at low-field. In conclusion, numerical simulations are used to model and refine the magnitudes and signs of $J$-couplings, as well as indirectly determine the 14 N relaxation rates in a single 1D experiment that has a higher fidelity than observed in high-field NMR experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lieb-Schultz-Mattis, Luttinger, and 't Hooft - anomaly matching in lattice systems

We analyze lattice Hamiltonian systems whose global symmetries have ’t Hooft anomalies. As is common in the study of anomalies, they are probed by coupling the system to classical background gauge fields. For flat fields (vanishing field strength), the nonzero spatial components of the gauge fields can be thought of as twisted boundary conditions, or equivalently, as topological defects. The symmetries of the twisted Hilbert space and their representations capture the anomalies. We demonstrate this approach with a number of examples. In some of them, the anomalous symmetries are internal symmetries of the lattice system, but they do not act on-site. (We clarify the notion of “on-site action.”) In other cases, the anomalous symmetries involve lattice translations. Using this approach we frame many known and new results in a unified fashion. In this work, we limit ourselves to 1+1d systems with a spatial lattice. In particular, we present a lattice system that flows to the c=1 c = 1 compact boson system with any radius (no BKT transition) with the full internal symmetry of the continuum theory, with its anomalies and its T-duality. As another application, we analyze various spin chain models and phrase their Lieb-Shultz-Mattis theorem as an ’t Hooft anomaly matching condition. We also show in what sense filling constraints like Luttinger theorem can and cannot be viewed as reflecting an anomaly. As a by-product, our understanding allows us to use information from the continuum theory to derive some exact results in lattice model of interest, such as the lattice momenta of the low-energy states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Polariton enhanced free charge carrier generation in donor–acceptor cavity systems by a second-hybridization mechanism

Cavity quantum electrodynamics has been studied as a potential approach to modify free charge carrier generation in donor–acceptor heterojunctions because of the delocalization and controllable energy level properties of hybridized light–matter states known as polaritons. However, in many experimental systems, cavity coupling decreases charge separation. In this work, we theoretically study the quantum dynamics of a coherent and dissipative donor–acceptor cavity system, to investigate the dynamical mechanism and further discover the conditions under which polaritons may enhance free charge carrier generation. We use open quantum system methods based on single-pulse pumping to find that polaritons have the potential to connect excitonic states and charge separated states, further enhancing free charge generation on an ultrafast timescale of several hundred femtoseconds. The mechanism involves polaritons with optimal energy levels that allow the exciton to overcome the high Coulomb barrier induced by electron–hole attraction. Moreover, we propose that a second-hybridization between a polariton state and dark states with similar energy enables the formation of the hybrid charge separated states that are optically active. These two mechanisms lead to a maximum of 50% enhancement of free charge carrier generation on a short timescale. However, our simulation reveals that on the longer timescale of picoseconds, internal conversion and cavity loss dominate and suppress free charge carrier generation, reproducing the experimental results. Thus, our work shows that polaritons can affect the charge separation mechanism and promote free charge carrier generation efficiency, but predominantly on a short timescale after photoexcitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗

Implementation of Chemistry in the Athena++ Code

Chemistry plays a key role in many aspects of astrophysical fluids. Atoms and molecules are agents for heating and cooling, determine the ionization fraction, serve as observational tracers, and build the molecular foundation of life. We present the implementation of a chemistry module in the publicly available magnetohydrodynamic code Athena++ . We implement several chemical networks and heating and cooling processes suitable for simulating the interstellar medium (ISM). A general chemical network framework in the KIDA format is also included, allowing users to easily implement their own chemistry. Radiation transfer and cosmic-ray ionization are coupled with chemistry and solved with the simple six-ray approximation. The chemical and thermal processes are evolved as a system of coupled ordinary differential equations with an implicit solver from the CVODE library. We perform and present a series of tests to ensure the numerical accuracy and convergence of the code. Many tests combine chemistry with gas dynamics, including comparisons with analytic solutions, 1D problems of the photodissociation regions and shocks, and realistic 3D simulations of the turbulent ISM. We release the code with the new public version of Athena++ , aiming to provide a robust and flexible code for the astrochemical simulation community.

79 ASTRONOMY AND ASTROPHYSICS↗

Mass-conserving implicit–explicit methods for coupled compressible Navier–Stokes equations

Earth system models are composed of coupled components that separately model systems such as the global atmosphere, ocean, and land surface. While these components are well developed, coupling them in a single system can be a significant challenge. Computational efficiency, accuracy, and stability are principal concerns. In this study we focus on these issues. In particular, implicit–explicit (IMEX) tight and loose coupling strategies are explored for handling different time scales. For a simplified model for the air–sea interaction problem, we consider coupled compressible Navier–Stokes equations with an interface condition. Under the rigid-lid assumption, horizontal momentum and heat flux are exchanged through the interface. Several numerical experiments are presented to demonstrate the stability of the coupling schemes. Here, we show both numerically and theoretically that our IMEX coupling methods are mass conservative for a coupled compressible Navier–Stokes system with the rigid-lid condition

42 ENGINEERING↗