Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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 649 records · Page 36

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Managing Software Provenance to Enhance Reproducibility in Computational Research

Scientific processes rely on software as an important tool for data acquisition, analysis, and discovery. Over the years, sustainable software development practices have made progress in being considered as an integral component of research. However, management of computation-based scientific studies is often left to individual researchers who design their computational experiments based on personal preferences and the nature of the study. Here, we believe that the quality, efficiency, and reproducibility of computation-based scientific research can be improved by explicitly creating an execution environment that allows researchers to provide a clear record of traceability. This is particularly relevant to complex computational studies in high-performance computing (HPC) environments. In this article, we review the documentation required to maintain a comprehensive record of HPC computational experiments for reproducibility. We also provide an overview of tools and practices that we have developed to perform such studies around Flash-X, a multiphysics scientific software.

97 MATHEMATICS AND COMPUTING↗

Simulation of Proton Beams in the AGS Booster with Space Charge

As higher intensity accelerators are designed, built and operated, the resonances and instabilities caused by space charge can become dominant effects in limiting the intensity of the beam that can be delivered. Particle loss has many causes but space charge can be a major factor in reducing beam intensity and generating emittance growth, particularly at low energies and in the lower energy stages of accelerators such as the AGS Booster. Successful operation relies on a good understanding of basic particle dynamics and it becomes an important challenge to ascertain how the beam evolves when charged particles experience both external forces through the fields of the magnets combined with internal forces through Coulomb interactions. A possible approach to such a study is to include calculations of space charge in computer simulation codes that track the motions of charged particles in 6D phase space. Unfortunately, repeated space charge calculations in 6D phase space are computationally intensive even with the most modern computing facilities. Study with traditional models of space charge is then either limited to a short real-time scale, or sacrifices accuracy and often omits some of the physics. The study described in this report attempts to explore all aspects of the Booster from injection, through RF capture and acceleration to the extraction energy. The simulations are allied to real machine settings. They do not attempt to model effects that are well known and for which accurate experimental data are available, but they try in many cases to extend beyond the normal Booster parameter range so as to understand the limitations on operations and suggest areas for improvement.

43 PARTICLE ACCELERATORS↗

Hydroxyl radical‐driven transformations of bisphenol A and 2,4‐dinitroanisole: Experimental and computational analysis

Abstract This study used experimental and computational analysis to investigate the advanced oxidation of bisphenol A (BPA) and 2,4‐dinitroanisole (DNAN). The pseudo first‐order reaction rate constants depended on the molar peroxide ratio and were between 0.13 and 0.28 min −1 for BPA and between 0.018 and 0.032 min −1 for DNAN. The kinetic differences appear to be due in part to the energy requirements for oxidation, which depended on the reaction mechanism but were typically lower for BPA than they were for DNAN. Density functional theory (DFT) was used to develop transformation pathways that included experimentally‐detected byproducts. The most energetically favored pathway for BPA oxidation begins with the formation of hydroxylated derivatives, while for DNAN, the most energetically favorable degradation pathway begins with the substitution of the methoxy group. Overall, these findings demonstrate the power of combining experimental and computational tools to reveal transformation mechanisms during water treatment. Practitioner Points Advanced oxidation transformations for two emerging water pollutants, bisphenol A and dinitroanisole, was investigated. The observed reaction kinetics depended on molar peroxide ratio in a manner that is in keeping with previous findings. Density functional theory‐based analysis revealed reaction energy requirements and degradation pathways.

Engineering↗

Increasing CO 2 Capture Rate in Liquid-Solvent Direct-Air Carbon Capture via Additive Manufacturing

Carbon capture is essential to meeting climate change mitigation goals. One approach currently being commercialized utilizes liquid-based solvents to capture CO 2 directly from the atmosphere but is limited by slow absorption of CO 2 into the liquid. Improved air/solvent liquid mixing increases CO 2 absorption rate, and this increased CO 2 absorption efficiency allows for smaller carbon capture systems with lower capital costs and better economic viability. In this project, we study the use of passive micromixers fabricated by metal additive manufacturing. The micromixer’s small-scale surface geometric features perturb and mix the liquid film to enhance mass transfer and CO 2 absorption. In this project, we evaluated this hypothesis through computational and experimental studies. Computational investigations focused on developing capabilities to simulate thin film (~ 100μm) fluid flow on rough surfaces. Such thin films are in a surface-tension dominated regime and simulations in this regime are prone to instabilities. Improvements to the Nalu code completed in this project resulted in a 10x timestep stability improvement for these problems.

36 MATERIALS SCIENCE↗

Tuning Two-Dimensional Phthalocyanine Dual Site Metal–Organic Framework Catalysts for the Oxygen Reduction Reaction

Metal-organic frameworks (MOFs) offer an interesting opportunity for catalysis, particularly for metal-nitrogen-carbon (M-N-C) motifs by providing an organized porous structural pattern and well-defined active sites for the oxygen reduction reaction (ORR), a key need for hydrogen fuel cells and related sustainable energy technologies. Here, in this work, we leverage electrochemical testing with computational models to study the electronic and structural properties in these systems and their relationship to ORR activity and stability based on dual transitional metal centers. These consists of two M1 metals with amine nodes coordinated to a single M2 metal with a phthalocyanine linker, where M1/M2 = Co, Ni, or Cu. Co-based metal centers, in particular Ni-Co, demonstrate the highest overall activity of all nine tested MOFs. Computationally, we identify the dominance of Co-sites, relative higher importance of the M2 site, and the role of layer M1 interactions on the ORR activity. Selectivity measurements indicate that M1 sites of MOFs, particularly Co, exhibits lowest (< 4%), and Ni demonstrates highest (>46%) two-electron selectivity, in good agreement with computational studies. Direct in-situ stability characterization, measuring dissolved metal ions, and calculations, using an alkaline stability metric, confirm that Co is the most stable metal in the MOF, while Cu exhibits notable instability at the M1. Overall, this study reveals how atomistic coupling of electronic and structural properties affects the ORR performance of dual site MOF catalysts and opens new avenues for tunable design and future development of these systems for practical electrochemical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automated, reliable, and efficient continental-scale replication of 7.3 petabytes of computational simulation data: A case study

We report on our experiences replicating 7.3 petabytes (PB) of Earth System Grid Federation (ESGF) computational simulation data from Lawrence Livermore National Laboratory (LLNL) in California to Argonne National Laboratory (ANL) in Illinois and Oak Ridge National Laboratory (ORNL) in Tennessee—a task motivated by a need for increased reliability, capacity, and performance. This task presented significant challenges: the need to move 29 million files twice under time pressure from aging storage hardware; a source file system bottleneck limiting throughput to 1.5 GB/s; frequent site maintenance windows; and the need for complete reliability at scale. We addressed these challenges using a simple replication tool that invoked Globus to transfer large bundles of files while tracking progress in a database, dynamically rerouting transfers to work around maintenance periods and file system limitations. Under the covers, Globus organized transfers to make efficient use of the high-speed Energy Sciences network (ESnet) and the data transfer nodes deployed at participating sites, and also addressed security, integrity checking, and recovery from a variety of transient failures. This success demonstrates the considerable benefits that can accrue from the adoption of performant data replication infrastructure. The replication tool is available at https://github.com/esgf2-us/data-replication-tools.

Globus↗

Exploring cycling induced crystallographic change in NMC with X-ray diffraction computed tomography

This study presents the application of X-ray diffraction computed tomography for the first time to analyze the crystal dimensions of LiNi 0.33 Mn 0.33 Co 0.33 O 2 electrodes cycled to 4.2 and 4.7 V in full cells with graphite as negative electrodes at 1 μm spatial resolution to determine the change in unit cell dimensions as a result of electrochemical cycling. The nature of the technique permits the spatial localization of the diffraction information in 3D and mapping of heterogeneities from the electrode to the particle level. An overall decrease of 0.4% and 0.6% was observed for the unit cell volume after 100 cycles for the electrodes cycled to 4.2 and 4.7 V. Additionally, focused ion beam-scanning electron microscope cross-sections indicate extensive particle cracking as a function of upper cut-off voltage, further confirming that severe cycling stresses exacerbate degradation. Finally, the technique facilitates the detection of parts of the electrode that have inhomogeneous lattice parameters that deviate from the bulk of the sample, further highlighting the effectiveness of the technique as a diagnostic tool, bridging the gap between crystal structure and electrochemical performance.

25 ENERGY STORAGE↗

S -matrix positivity without Lorentz invariance: a case study

We investigate the analytic structure of scattering amplitudes in theories in which Lorentz invariance is spontaneously broken. We do so by computing and studying the S-matrix for a simple example: a superfluid described by a complex scalar with quartic interactions. The computation is confined to tree-level, for there are no absolutely stable single-particle states, though the lifetime can be made long by lowering the chemical potential. For the 2 → 2 amplitude in center-of-mass configurations, not only is crossing symmetry violated, there appears a tree level branch cut for unphysical kinematics. Its appearance is a consequence of non-analyticity in the dispersion relation. The branch point defines a new scale in the problem, which scales inversely with the chemical potential. In this example, even derivatives of the forward amplitude are positive while odd derivatives are negative. This pattern can be understood in a general way in the limit of a small chemical potential, or weak Lorentz breaking.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Test Points for Online Monitoring of Quantum Circuits

Noisy Intermediate-Scale Quantum (NISQ) computers consisting of tens of inherently noisy quantum bits (qubits) suffer from reliability problems. Qubits and their gates are susceptible to various types of errors. Due to limited numbers of qubits and high error rates, quantum error correction cannot be applied. Physical constraints of quantum hardware including the error rates are used to guide the design and the layout of quantum circuits. The error rates determine the selection of qubits and their operations. The resulting circuit is executed on the quantum computer. This study explores the risk of unexpected changes in the error rates of NISQ computers post-calibration. We show that unexpected changes in error rates can alter the output state of a quantum circuit. To detect these changes, we propose the insertion of test points into the quantum circuit to enable online monitoring of the physical qubit behavior. We utilize classical, superposition, and uncompute test points. Furthermore, we use a gate error coverage metric to assess the quality of the tests. We verify the effectiveness of the proposed scheme on different IBM quantum computers (IBM Q), in addition to a noisy simulation that shows the scalability of the proposed approach.

97 MATHEMATICS AND COMPUTING↗

Curifactory: A research experiment manager

Curifactory is a command line tool and framework for organizing Python experiment code, configuration parameters, and results. It is an opinionated and lightweight approach to workflow management infrastructure and is primarily intended to support researchers conducting experiments on one machine. This software was developed to support the reproducibility of results for several data science projects in the Nuclear Nonproliferation Division at Oak Ridge National Laboratory. Curifactory is intended to be a general framework and is not specific to machine learning or data science. It can aid in any field in which experiments are primarily computation-based studies and can be implemented in Python (e.g., high-energy physics, astronomy, computational chemistry). Here, the design emphasizes the automated caching of intermediate data analysis artifacts to speed up development involving computationally intensive tasks. It also allows for data provenance and experiment reproduction. Individual experiment runs are tracked through logs and their output reports, and entire copies of a run with all cached data and metadata can be exported for others to run using Curifactory on another machine. Curifactory experiments can either be integrated into a project from the beginning or can be written on top of an existing codebase without needing significant modification. A few important views of the Curifactory library can be seen in Figure 1.

97 MATHEMATICS AND COMPUTING↗

Scattering using real-time path integrals

Background: Path integrals are a powerful tool for solving problems in quantum theory that are not amenable to a treatment by perturbation theory. Most path integral computations require an analytic continuation to imaginary time. While imaginary time treatments of scattering are possible, imaginary time is not a natural framework for treating scattering problems. More importantly, quantum algorithms for calculating path integrals require real-time evolution. Purpose: Here, we test a recently introduced method for performing direct calculations of scattering observables using real-time path integrals in order to understand the challenges facing real-time path integral calculations of scattering observables. Method: The computations are based on a new interpretation of the path integral as the expectation value of a potential functional on cylinder sets of continuous paths with respect to a complex probability distribution. The method can in principle be applied to arbitrary short-range potentials. Results: The method is applied to compute matrix elements of Møller wave operators applied to narrow wave packets. These are used to calculate half-shell sharp-momentum transition matrix elements for one-dimensional potential scattering. The calculations for half-shell transition operator matrix elements converge to the numerical solution of the Lippmann-Schwinger equation. Conclusions: This work presents a proof in principle that scattering observables can be computed using real-time Feynman path integrals. While the computational method is not efficient, it can be improved. It provides a laboratory for studying quantum computational algorithms that are applicable to scattering problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Significance of Aggregation Methods in Functional Group Modeling

The growth of forests and the feedbacks between forests and environmental changes are central issues in the planetary carbon cycle, global climate change, and basic plant ecology. A challenge to understanding both growth and feedbacks from local to global scales is that many critical metabolic processes vary among species. An innovation in solving this challenge is the recognition that species can be lumped into “functional groups” based on metabolic similarity, and these functional groups can then be studied in computational models that simulate ecosystem function. Despite the vast resources devoted to functional group studies and the progress made by them, an important logical and biological question has not been formally addressed, “How do the groupings alter the results of modeling studies?” To what extent do modeling results depend on the choices made in aggregating taxa into functional groups. Here, we consider the effects of using different aggregation strategies in simulating the carbon dynamics of a deciduous forest. Understanding the impacts that aggregation strategy has on efforts to simulate regional-to-global-scale forest dynamics offers insights into both ecosystem regulation and model function and addresses this central problem in the study of carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

Dimension Reduction and Redundancy Removal through Successive Schmidt Decompositions

Quantum computers are believed to have the ability to process huge data sizes, which can be seen in machine learning applications. In these applications, the data, in general, are classical. Therefore, to process them on a quantum computer, there is a need for efficient methods that can be used to map classical data on quantum states in a concise manner. On the other hand, to verify the results of quantum computers and study quantum algorithms, we need to be able to approximate quantum operations into forms that are easier to simulate on classical computers with some errors. Motivated by these needs, in this paper, we study the approximation of matrices and vectors by using their tensor products obtained through successive Schmidt decompositions. We show that data with distributions such as uniform, Poisson, exponential, or similar to these distributions can be approximated by using only a few terms, which can be easily mapped onto quantum circuits. The examples include random data with different distributions, the Gram matrices of iris flower, handwritten digits, 20newsgroup, and labeled faces in the wild. Similarly, some quantum operations, such as quantum Fourier transform and variational quantum circuits with a small depth, may also be approximated with a few terms that are easier to simulate on classical computers. Furthermore, we show how the method can be used to simplify quantum Hamiltonians: In particular, we show the application to randomly generated transverse field Ising model Hamiltonians. The reduced Hamiltonians can be mapped into quantum circuits easily and, therefore, can be simulated more efficiently.

97 MATHEMATICS AND COMPUTING↗