Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Expression systems”

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 1,027 records · Page 57

Achieving performance portability in Gaussian basis set density functional theory on accelerator based architectures in NWChemEx

The numerical integration of the exchange–correlation (XC) potential is one of the primary computational bottlenecks in Gaussian basis set Kohn–Sham density functional theory (KS-DFT). To achieve optimal performance and accuracy, care must be taken in this numerical integration to preserve local sparsity as to allow for near linear weak scaling with system size. This leads to an integration scheme with several performance critical kernels which must be hand optimized for each architecture of interest. As the set of available accelerator hardware goes more diverse, a key challenge for developers of KS-DFT software is to maintain performance portability across a wide range of computational architectures. In this article, we examine a modular software design pattern which decouples the implementation details of performance critical kernels from the expression of high-level algorithmic workflows in a device-agnostic language such as C++; thus allowing for developers to target existing and emerging accelerator hardware within a single code base. We consider the efficacy of such a design pattern in the numerical integration of the XC potential by demonstrating its ability to achieve performance portability across a set of accelerator architectures which are representative of those on current and future U.S. Department of Energy Leadership Computing Facilities.

97 MATHEMATICS AND COMPUTING↗

Double Inozemtsev limits of the quantum DELL system

In this letter we study various Inozemtsev-type limits of the quantum double elliptic (DELL) system when both elliptic parameters are sent to zero at different rates, while the coupling constant is sent to infinity, such that a certain combination of the three parameters is kept fixed. We find a regime in which such double Inozemtsev limit of DELL produces the elliptic Ruijsenaars-Schneider (eRS) Hamiltonians albeit in an unconventional normalization. We discuss other double scaling limits and anisotropic scaling of coordinates and momenta. In addition we provide a formal expression for the eigenvalues of the eRS Hamiltonians solely in terms of their eigenfunctions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parallel Runtime Interface for Fortran (PRIF) Specification (Rev. 0.3)

This document specifies an interface to support the parallel features of Fortran, named the Parallel Runtime Interface for Fortran (PRIF). PRIF is a proposed solution in which the runtime library is responsible for coarray allocation, deallocation and accesses, image synchronization, atomic operations, events, and teams. In this interface, the compiler is responsible for transforming the invocation of Fortran-level parallel features into procedure calls to the necessary PRIF procedures. The interface is designed for portability across shared- and distributed-memory machines, different operating systems, and multiple architectures. Implementations of this interface are intended as an augmentation for the compiler's own runtime library. With an implementation-agnostic interface, alternative parallel runtime libraries may be developed that support the same interface. One benefit of this approach is the ability to vary the communication substrate. A central aim of this document is to define a parallel runtime interface in standard Fortran syntax, which enables us to leverage Fortran to succinctly express various properties of the procedure interfaces, including argument attributes.

97 MATHEMATICS AND COMPUTING↗

Parallel Runtime Interface for Fortran (PRIF) Specification (Rev. 0.5)

This document specifies an interface to support the parallel features of Fortran, named the Parallel Runtime Interface for Fortran (PRIF). PRIF is a proposed solution in which the runtime library is primarily responsible for implementing coarray allocation, deallocation and accesses, image synchronization, atomic operations, events, teams and collective subroutines. In this interface, the compiler is responsible for transforming the invocation of Fortran-level parallel features into procedure calls to the necessary PRIF subroutines. The interface is designed for portability across shared- and distributed-memory machines, different operating systems, and multiple architectures. Implementations of this interface are intended as an augmentation for the compiler's own runtime library. With an implementation-agnostic interface, alternative parallel runtime libraries may be developed that support the same interface. One benefit of this approach is the ability to vary the communication substrate. A central aim of this document is to define a parallel runtime interface in standard Fortran syntax, which enables us to leverage Fortran to succinctly express various properties of the procedure interfaces, including argument attributes.

97 MATHEMATICS AND COMPUTING↗

Parallel Runtime Interface for Fortran (PRIF) Design Document (Rev. 0.2)

This design document proposes an interface to support the parallel features of Fortran, named the Parallel Runtime Interface for Fortran (PRIF). PRIF is a proposed solution in which the runtime library is responsible for coarray allocation, deallocation and accesses, image synchronization, atomic operations, events, and teams. In this interface, the compiler is responsible for transforming the invocation of Fortran-level parallel features into procedure calls to the necessary PRIF procedures. The interface is designed for portability across shared- and distributed-memory machines, different operating systems, and multiple architectures. Implementations of this interface are intended as an augmentation for the compiler's own runtime library. With an implementation-agnostic interface, alternative parallel runtime libraries may be developed that support the same interface. One benefit of this approach is the ability to vary the communication substrate. A central aim of this document is to define a parallel runtime interface in standard Fortran syntax, which enables us to leverage Fortran to succinctly express various properties of the procedure interfaces, including argument attributes.

97 MATHEMATICS AND COMPUTING↗

Parallel Runtime Interface for Fortran (PRIF) Specification (Rev. 0.6)

This document specifies an interface to support the multi-image parallelism features of Fortran, named the Parallel Runtime Interface for Fortran (PRIF). PRIF is a solution in which a runtime library is primarily responsible for implementing coarray allocation, deallocation and accesses, image synchronization, atomic operations, events, teams and collective subroutines. The Fortran compiler is responsible for transforming the invocation of Fortran-level multi-image parallelism features into procedure calls to the necessary PRIF subroutines. The interface is designed for portability across shared- and distributed-memory machines, different operating systems, and multiple architectures. Implementations of this interface are intended as an augmentation for the compiler's own runtime library. With an implementation-agnostic interface, alternative parallel runtime libraries may be developed that support the same interface. One benefit of this approach is the ability to vary the communication substrate. A central aim of this document is to define a parallel runtime interface in standard Fortran syntax, which enables us to leverage Fortran to succinctly express various properties of the procedure interfaces, including argument attributes.

97 MATHEMATICS AND COMPUTING↗

Parallel Runtime Interface for Fortran (PRIF) Specification (Rev. 0.4)

This document specifies an interface to support the parallel features of Fortran, named the Parallel Runtime Interface for Fortran (PRIF). PRIF is a proposed solution in which the runtime library is responsible for coarray allocation, deallocation and accesses, image synchronization, atomic operations, events, and teams. In this interface, the compiler is responsible for transforming the invocation of Fortran-level parallel features into procedure calls to the necessary PRIF procedures. The interface is designed for portability across shared- and distributed-memory machines, different operating systems, and multiple architectures. Implementations of this interface are intended as an augmentation for the compiler's own runtime library. With an implementation-agnostic interface, alternative parallel runtime libraries may be developed that support the same interface. One benefit of this approach is the ability to vary the communication substrate. A central aim of this document is to define a parallel runtime interface in standard Fortran syntax, which enables us to leverage Fortran to succinctly express various properties of the procedure interfaces, including argument attributes.

97 MATHEMATICS AND COMPUTING↗

Biphasic response of human iPSC-derived neural network activity following exposure to a sarin-surrogate nerve agent

Organophosphorus nerve agents (OPNA) are hazardous environmental exposures to the civilian population and have been historically weaponized as chemical warfare agents (CWA). OPNA exposure can lead to several neurological, sensory, and motor symptoms that can manifest into chronic neurological illnesses later in life. There is still a large need for technological advancement to better understand changes in brain function following OPNA exposure. The human-relevant in vitro multi-electrode array (MEA) system, which combines the MEA technology with human stem cell technology, has the potential to monitor the acute, sub-chronic, and chronic consequences of OPNA exposure on brain activity. However, the application of this system to assess OPNA hazards and risks to human brain function remains to be investigated. In a concentration-response study, we have employed a human-relevant MEA system to monitor and detect changes in the electrical activity of engineered neural networks to increasing concentrations of the sarin surrogate 4-nitrophenyl isopropyl methylphosphonate (NIMP). We report a biphasic response in the spiking (but not bursting) activity of neurons exposed to low (i.e., 0.4 and 4 μM) versus high concentrations (i.e., 40 and 100 μM) of NIMP, which was monitored during the exposure period and up to 6 days post-exposure. Regardless of the NIMP concentration, at a network level, communication or coordination of neuronal activity decreased as early as 60 min and persisted at 24 h of NIMP exposure. Once NIMP was removed, coordinated activity was no different than control (0 μM of NIMP). Interestingly, only in the high concentration of NIMP did coordination of activity at a network level begin to decrease again at 2 days post-exposure and persisted on day 6 post-exposure. Notably, cell viability was not affected during or after NIMP exposure. Also, while the catalytic activity of AChE decreased during NIMP exposure, its activity recovered once NIMP was removed. Gene expression analysis suggests that human iPSC-derived neurons and primary human astrocytes resulted in altered genes related to the cell’s interaction with the extracellular environment, its intracellular calcium signaling pathways, and inflammation, which could have contributed to how neurons communicated at a network level.

59 BASIC BIOLOGICAL SCIENCES↗

Entropy and de Haas–van Alphen oscillations of a three-dimensional marginal Fermi liquid

Here we study de Haas-van Alphen oscillations in a marginal Fermi liquid resulting from a three-dimensional metal tuned to a quantum-critical point (QCP). We show that the conventional approach based on extensions of the Lifshitz-Kosevich formula for the oscillation amplitudes becomes inapplicable when the correlation length exceeds the cyclotron radius. This breakdown is due to (i) non-analytic finite-temperature contributions to the fermion self-energy (ii) an enhancement of the oscillatory part of the self-energy by quantum fluctuations, and (iii) non-trivial dynamical scaling laws associated with the quantum critical point. We properly incorporate these effects within the Luttinger-Ward-Eliashberg framework for the thermodynamic potential by treating the fermionic and bosonic contributions on equal footing. As a result, we obtain the modified expressions for the oscillations of entropy and magnetization that remain valid in the non-Fermi liquid regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantifying the quantum heat contribution from a driven superconducting circuit

Heat flow management at the nanoscale is of great importance for emergent quantum technologies. For instance, a thermal sink that can be activated on-demand is a highly desirable tool that may accommodate the need to evacuate excess heat at chosen times, e.g., to maintain cryogenic temperatures or reset a quantum system to ground, and the possibility of controlled unitary evolution otherwise. Here we propose a design of such heat switch based on a single coherently driven qubit. We show that the heat flow provided by a hot source to the qubit can be switched on and off by varying external parameters, the frequency and the intensity of the driving. The complete suppression of the heat flow is a quantum effect occurring for specific driving parameters that we express and we analyze the role of the coherences in the free-qubit energy eigenbasis. We finally study the feasibility of this quantum heat switch in a circuit QED setup involving a charge qubit coupled to thermal resistances. In this work, we demonstrate robustness to experimental imperfections such as additional decoherence, paving the road towards experimental verification of this effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Statistical Performance of Forced Oscillation Detectors in the Presence of Missing Measurements

In bulk power systems, measurement-based monitoring for large oscillations can help maintain system reliability. One of the challenges encountered in a recent field demonstration was the unavailability of measurements due to underlying measurement quality or communication problems. During the demonstration, the oscillation detector ignored a measurement location if even 10 seconds of data was missing. To extend the detector's ability to operate in these conditions, this paper evaluates the impact of three methods for addressing missing data. The strengths and weaknesses of each approach are evaluated using theoretical expressions for the probability of detection along with results from simulated data and publicly available field measurements. Based on these results, a suitable approach is identified that can extend the oscillation detector's performance when large segments of data are missing.

Follum, James D.↗

Nuclear Data Libraries Sensitivity Studies for ORSA Using SCALE

Subcritical assemblies offer valuable training capabilities in nuclear criticality safety (NCS) for individuals handling fissile material. At Oak Ridge National Laboratory(ORNL), the Oak Ridge Subcritical Assembly (ORSA), a new experimental facility, is being established to provide hands-on training for the Nuclear Criticality Safety Program (NCSP).It is essential to accurately determine the neutron multiplication factor (keff) to ensure that ORSA remains subcritical and safe during operations. This study investigated the sensitivity of keff to variations across Evaluated Nuclear Data File (ENDF/B) libraries, consisting of ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1. The analysis was conducted using the CSAS6 sequence in the SCALE-6.3 code system. Individual isotopes in the ORSA model were replaced one at a time with ENDF/B-VII.1 as the base library and changing to ENDF/B-VIII.0or ENDF/B-VIII.1. The results demonstrated that the changes in keffof the nuclides associated with the ORSA model were mostly within the uncertainty of the base model (~24 pcm), except for primary nuclides like Uranium-235and H-poly with few other nuclides. The relative delta keff values of Uranium-235 and H-poly, expressed in pcm, were +280 and -212 in ENDF/B-VIII.0 and +332 and -275 in ENDF/B-VIII.1, respectively, which were notable changes in reactivity. These results demonstrated that ORSA was largely insensitive to variations across these nuclear data libraries.

Hong, Evan [North Carolina State University]↗

Assessing parallel path cooling tower performance via artificial neural networks

Real-time monitoring of a research nuclear reactor, a system in which all generated power is dissipated to the environment, can be performed via analysis of the heat rejection from the cooling system. Given an inlet water temperature and flow rate, the reactor power can be well-approximated from the outlet water temperature; however, the instrumentation to measure outlet conditions may not be robust or accurate. If we know how a cooling tower performs from historical data, but cannot measure the outlet temperature, a mathematical representation of the system can be inverted to obtain the outlet water temperature that describes the cooling capacity. Unfortunately, model inversion processes are computationally expensive. To address this, an artificial neural network (ANN) is implemented to assess the performance of a multi-cell cooling tower for a nuclear reactor. This approach leverages the Merkel model to obtain an extensive data set describing performance of the cooling tower cells throughout a wide array of potential operating conditions. The Merkel model is expressed as a function of four parameters: the inlet and outlet water temperatures, inlet air wet bulb temperature, and ratio of liquid-to-gas mass flow rates (L/G), which together provide a non-dimensional number indicative of cooling tower performance, called the Merkel integral. Computing a 4-dimensional data structure that describes finite combinations of the Merkel integral, an inverse model is then generated using an ANN to determine the cell outlet water temperature from the other three model parameters along with the computed Merkel integral. Compared to traditional model inversion methods, the ANN reduces the computational time by approximately 4 orders of magnitude, with effectively no sacrifice to solution accuracy, and could be applied for different cooling towers in the event the performance curve is known. Finally, three use cases of the ANN are then reviewed: (1) determining the cell outlet water temperatures when gas flow at rated conditions (GFRC) is known, (2) performing the prior case without knowledge of the GRFC, and (3) assessing performance differences between the individual tower cells.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Peridynamic neural operators: A data-driven nonlocal constitutive model for complex material responses

Neural operators, which can act as implicit solution operators of hidden governing equations, have recently become popular tools for learning the responses of complex real-world physical systems. Nevertheless, most neural operator applications have thus far been data-driven and neglect the intrinsic preservation of fundamental physical laws in data. Here, in this work, we introduce a novel integral neural operator architecture called the Peridynamic Neural Operator (PNO) that learns a nonlocal constitutive law from data. This neural operator provides a forward model in the form of state-based peridynamics, with objectivity and momentum balance laws automatically guaranteed. As applications, we demonstrate the expressivity and efficacy of our model in learning complex material behaviors from both synthetic and experimental data sets. We also compare the performances with baseline models that use predefined constitutive laws. We show that, owing to its ability to capture complex responses, our learned neural operator achieves improved accuracy and efficiency. Moreover, by preserving the essential physical laws within the neural network architecture, the PNO is robust in treating noisy data. The method shows generalizability to different domain configurations, external loadings, and discretizations.

42 ENGINEERING↗

Normal modes, rotational inertia, and thermal fluctuations of trapped ion crystals

The normal modes of a trapped ion crystal are derived using an approach based on the Hermitian properties of the system’s dynamical matrix. This method is equivalent to the standard Bogoliubov method, but for classical systems it is arguably simpler and more general in that canonical coordinates are not necessary. The theory is developed for stable, unstable, and neutrally-stable systems. The method is then applied to ion crystals in a Penning trap. Reduced eigenvalue problems for the case of large applied magnetic field are developed, for which the spectrum breaks into E x B drift modes, axial modes, and cyclotron modes. Thermal fluctuation levels in these modes are analyzed and shown to be consistent with the Bohr-van-Leeuwen theorem, provided that neutrally-stable modes associated with crystal rotations are included in the analysis. An expression for the rotational inertia of the crystal is derived, and a magnetic contribution to this inertia, which dominates in large magnetic fields, is described. An unusual limit is discovered for the special case of spherically-symmetric confinement, in which the rotational inertia does not exist and changes in angular momentum leave the rotation frequency unaffected.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying multipoint ordering in alloys

A central problem in multicomponent lattice systems is to systematically quantify multipoint ordering. Ordering in such systems is often described in terms of pairs, even though this is not sufficient when three-point and higher-order interactions are included in the Hamiltonian. Current models and parameters for multipoint ordering are often only applicable for very specific cases or require approximating a subset of correlated occupational variables on a lattice as being uncorrelated. In this paper, cluster order parameters are introduced to systematically quantify arbitrary multipoint ordering motifs in substitutional systems through direct calculations of normalized cluster probabilities. These parameters can describe multipoint chemical ordering in crystal systems with multiple sublattices, multiple components, and systems with reduced symmetry. These are defined in this paper and applied to quantify four-point chemical ordering motifs in platinum/palladium alloy nanoparticles that are of practical interest to the synthesis of catalytic nanocages. Impacts of chemical ordering on nanocage stability are discussed. It is demonstrated that approximating four-point probabilities from superpositions of lower-order pair probabilities is not sufficient in cases where three- and four-body terms are included in the energy expression. Furthermore, conclusions about the formation mechanisms of nanocages may change significantly when using common pair approximations.

36 MATERIALS SCIENCE↗

Kondo Lattice Model of Magic-Angle Twisted-Bilayer Graphene: Hund’s Rule, Local-Moment Fluctuations, and Low-Energy Effective Theory

For this work, we apply a generalized Schrieffer-Wolff transformation to the extended Anderson-like topological heavy fermion (THF) model for the magic-angle (θ=1.05°) twisted bilayer graphene (MATBLG) [Phys. Rev. Lett. 129, 047601 (2022)], to obtain its Kondo lattice limit. In this limit localized $\mathcal{f}$ electrons on a triangular lattice interact with topological conduction $\mathcal{c}$ electrons. By solving the exact limit of the THF model, we show that the integer fillings ν=0, ±1, ±2 are controlled by the heavy $\mathcal{f}$ electrons, while ν=±3 is at the border of a phase transition between two $\mathcal{f}$-electron fillings. For ν=0, ±1, ±2, we then calculate the Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions between the $\mathcal{f}$ moments in the full model and analytically prove the SU(4) Hund’s rule for the ground state which maintains that two $\mathcal{f}$ electrons fill the same valley-spin flavor. Our (ferromagnetic interactions in the) spin model dramatically differ from the usual Heisenberg antiferromagnetic interactions expected at strong coupling. We show the ground state in some limits can be found exactly by employing a positive semidefinite “bond-operators” method. We then compute the excitation spectrum of the $\mathcal{f}$ moments in the ordered ground state, prove the stability of the ground state favored by RKKY interactions, and discuss the properties of the Goldstone modes, the (reason for the accidental) degeneracy of (some of) the excitation modes, and the physics of their phase stiffness. We develop a low-energy effective theory for the $\mathcal{f}$ moments and obtain analytic expressions for the dispersion of the collective modes. We discuss the relevance of our results to the spin-entropy experiments in TBG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding the Quadrupole Mass Filter and Testing a High-Resolution QMS RGA for ITER

A common type of residual gas analyzer is the quadrupole mass spectrometer. One of the main components within this instrument is a mass filter known as the quadrupole. It is responsible for the selective throughput of the ionized gas particles - by ascending mass number - prior to ion impacts on the analyzer (or detector) surface from which the ion current signal is generated for processing. However, the quadrupole is not fully described in relation to the electric field characteristics and the function as an ion mass separator. This paper describes the basic origins of the electrical design, the intricate assembly criteria, and performance of the quadrupole within the spectrometer. A specialized quadrupole mass spectrometer is part of a configuration for a diagnostic gas analyzer system planned for ITER, a fusion research machine. It has a verified capability, essential as a diagnostic criterion for this reactor project, to successfully deconvolute the mass signals of Helium-4 and deuterium (reactor fuel exhaust gases, separated by only 0.026 atomic mass units), down to a relative three-percent concentration of the former gas. The associated preliminary testing, performed at the Oak Ridge National Laboratory, is also addressed. Finally, one of the key parameters used to express gas concentration, the relative sensitivity factor, will be explained, including an evaluation of dependency on other variables.

Marcus, Chris [ORNL] (ORCID:0000000190139636)↗