Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “future of computing”

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 721 records · Page 40

The Significance of the 'Insignificant': Non-covalent Interactions in CO 2 Reduction Reactions with 3C-TM (TM=Sc-Zn) Single-Atom Catalysts

With energy shortages and excessive CO 2 emissions driving climate change, converting CO 2 into high-value-added products offers a promising solution for carbon recycling. We investigate CO 2 reduction reactions (CO2RR) catalyzed by 10 single-atom catalysts (SACs), incorporating weak non-covalent interactions, specifically lone pair-π and H-π interactions. The SACs, consisting of transition metals coordinated by three carbon atoms in a defective graphene substrate (3C-TM, TM=Sc-Zn), leverage these interactions to influence the energy fluctuations of intermediates and the limiting potentials of CO 2 RR, without altering the overall reaction pathway. Further, our findings show that SACs based on early transition metals (Sc, Ti, V, Cr) can serve as catalysts for C 1 products, including HCOOH, HCHO, CH 3 OH, and CH 4 , while those based on Fe and Co are suitable for CO formation. Driving force analysis helps bridge theoretical results with experimental observations and propose a modified approach for assessing hydrogen evolution reactions (HER) competition. SACs based on Ni and Cu exhibit moderate HER tolerance, while early transition metals excel in selective CO 2 reduction. We also identify a linear scaling relationship between the free energies of *COOH and *CO. This study offers valuable insights for future experimental studies and large-scale computational screenings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bayesian model aggregation for ensemble-based estimates of protein pKa values

This paper investigates an ensemble-based technique called Bayesian Model Averaging (BMA) to improve the performance of protein amino acid p$K_a$ predictions. Structure-based p$K_a$ calculations play an important role in the mechanistic interpretation of protein structure and are also used to determine a wide range of protein properties. A diverse set of methods currently exist for p$K_a$ prediction, ranging from empirical statistical models to {\it ab initio} quantum mechanical approaches. However, each of these methods are based on a set of assumptions that have inherent bias and sensitivities that can effect a model's accuracy and generalizability for p$K_a$ prediction in complicated biomolecular systems. We use BMA to combine eleven diverse prediction methods that each estimate pKa values of amino acids in staphylococcal nuclease. These methods are based on work conducted for the pKa Cooperative and the pKa measurements are based on experimental work conducted by the Garc{\'i}a-Moreno lab. Our study demonstrates that the aggregated estimate obtained from BMA outperforms all individual prediction methods in our cross-validation study with improvements from 40-70\% over other method classes. This work illustrates a new possible mechanism for improving the accuracy of p$K_a$ prediction and lays the foundation for future work on aggregate models that balance computational cost with prediction accuracy.

electrostatics↗

Performance analysis of data intensive cloud systems based on data management and replication: a survey

As we delve deeper into the ‘Digital Age’, we witness an explosive growth in the volume, velocity, and variety of the data available on the Internet. For example, in 2012 about 2.5 quintillion bytes of data was created on a daily basis that originated from myriad of sources and applications including mobiledevices, sensors, individual archives, social networks, Internet of Things, enterprises, cameras, software logs, etc. Such ‘Data Explosions’ has led to one of the most challenging research issues of the current Information and Communication Technology era: how to optimally manage (e.g., store, replicated, filter, and the like) such large amount of data and identify new ways to analyze large amounts of data for unlocking information. It is clear that such large data streams cannot be managed by setting up on-premises enterprise database systems as it leads to a large up-front cost in buying and administering the hardware and software systems. Therefore, next generation data management systems must be deployed on cloud. The cloud computing paradigm provides scalable and elastic resources, such as data and services accessible over the Internet Every Cloud Service Provider must assure that data is efficiently processed and distributed in a way that does not compromise end-users’ Quality of Service (QoS) in terms of data availability, data search delay, data analysis delay, and the like. In the aforementioned perspective, data replication is used in the cloud for improving the performance (e.g., read and write delay) of applications that access data. Through replication a data intensive application or system can achieve high availability, better fault tolerance, and data recovery. In this paper, we survey data management and replication approaches (from 2007 to 2011) that are developed by both industrial and research communities. The focus of the survey is to discuss and characterize the existing approaches of data replication and management that tackle the resource usage and QoS provisioning with different levels of efficiencies. Moreover, the breakdown of both influential expressions (data replication and management) to provide different QoS attributes is deliberated. Furthermore, the performance advantages and disadvantages of data replication and management approaches in the cloud computing environments are analyzed. Open issues and future challenges related to data consistency, scalability, load balancing, processing and placement are also reported.

Malik, Saif Ur Rehman↗

Equation-based and data-driven modeling: Open-source software current state and future directions

Here, a review of current trends in scientific computing reveals a broad shift to open-source and higher-level programming languages such as Python and growing career opportunities over the next decade. Open-source modeling tools accelerate innovation in equation-based and data-driven applications. Significant resources have been deployed to develop data-driven tools (PyTorch, TensorFlow, Scikit-learn) from tech companies that rely on machine learning services to meet business needs while keeping the foundational tools open. Open-source equation-based tools such as Pyomo, CasADi, Gekko, and JuMP are also gaining momentum according to user community and development pace metrics. Integration of data-driven and principles-based tools is emerging. New compute hardware, productivity software, and training resources have the potential to radically accelerate progress. However, long-term support mechanisms are still necessary to sustain the momentum and maintenance of critical foundational packages.

97 MATHEMATICS AND COMPUTING↗

Policy support and technology development trajectory for renewable natural gas in the U.S.

Renewable natural gas (RNG) is a clean alternative to fossil natural gas, which can be used as transportation fuel, among other applications. This study projects the development trajectory of RNG and evaluates its impacts on the future U.S. transportation market using a hybrid computable general equilibrium model. This analysis considers various factors and uncertainties affecting RNG production, such as technology development, market conditions, competition with other advanced biofuels, and national and state policies. In 2050, RNG production will grow to 2.7 billion gallons (10 billion liters), mostly from swine manure, under current policy provisions. This will lead to a reduction in greenhouse gas (GHG) emissions by 58.56 million metric tonne of CO 2e in 2050. Analysis of different technology cases finds RNG from animal manure to be predominant, while RNG from corn stover and cellulosic ethanol are less competitive. Furthermore, a high mandatory target of 1 billion gallons will drive RNG production higher by 8–18 %, while an extended 2 nd -generation biofuel production tax credit will mostly increase cellulosic ethanol production. The model also finds RNG production being affected by uncertainties in market conditions, such as GDP growth, fossil fuel prices, and oil and gas supply.

Biomethane↗

Deciphering Distinct Overpotential-Dependent Pathways for Electrochemical CO 2 Reduction Catalyzed by an Iron–Terpyridine Complex

[Fe(tpyPY2Me)] 2+ ([Fe] 2+ ) is a homogeneous electrocatalyst for converting CO 2 into CO featuring low overpotentials of <100 mV, near-unity selectivity, and high activity with turnover frequencies faster than 100 000 s -1 . To identify the origins of its exceptional performance and inform future catalyst design, we report a combined computational and experimental study that establishes two distinct mechanistic pathways for electrochemical CO 2 reduction catalyzed by [Fe] 2+ as a function of applied overpotential. Electrochemical data shows the formation of two catalytic regimes at low (η TOF/2 of 160 mV) and high (η TOF/2 of 590 mV) overpotential plateaus. We propose that at low overpotentials [Fe] 2+ undergoes a two-electron reduction, two-proton-transfer mechanism (electrochemical-electrochemical-chemical-chemical, EECC), where turnover occurs through the dicationic iron complex, [Fe] 2+ . Computational analysis supports the importance of the singlet ground-state electronic structure for CO 2 binding and that the rate-limiting step is the second protonation in this low-overpotential regime. When more negative potentials are applied, an additional electron-transfer event occurs through either a stepwise or proton-coupled electron-transfer (PCET) pathway, enabling catalytic turnover from the monocationic iron complex ([Fe] + ) via an electrochemical-chemical-electrochemical-chemical (ECEC) mechanism. Comparison of experimental kinetic data obtained from variable controlled potential electrolysis (CPE) experiments with direct product detection with calculated rates obtained from the energetic span model supports the PCET pathway as the most likely mechanism. Moreover, we build upon this mechanistic understanding to propose the design of an improved ligand framework that is predicted to stabilize the key transition states identified in our study and explore their electronic structures using an energy decomposition analysis. Finally, taken together, this work highlights the value of synergistic computational/experimental approaches to decipher mechanisms of new electrocatalysts and direct the rational design of improved platforms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perovskite neural trees

Trees are used by animals, humans and machines to classify information and make decisions. Natural tree structures displayed by synapses of the brain involves potentiation and depression capable of branching and is essential for survival and learning. Demonstration of such features in synthetic matter is challenging due to the need to host a complex energy landscape capable of learning, memory and electrical interrogation. We report experimental realization of tree-like conductance states at room temperature in strongly correlated perovskite nickelates by modulating proton distribution under high speed electric pulses. This demonstration represents physical realization of ultrametric trees, a concept from number theory applied to the study of spin glasses in physics that inspired early neural network theory dating almost forty years ago. We apply the tree-like memory features in spiking neural networks to demonstrate high fidelity object recognition, and in future can open new directions for neuromorphic computing and artificial intelligence.

36 MATERIALS SCIENCE↗

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

RADIAL DEPENDENCE OF PEAK PROTON AND IRON ION FLUXES IN SOLAR ENERGETIC PARTICLE EVENTS: APPLICATION OF THE PATH CODE

The radial dependence of particle peak fluxes in large solar energetic particle (SEP) events is important in determining the potential impact of space weather hazards on space missions. Using the Particle Acceleration and Transport in the Heliosphere code, we model the acceleration and transport of protons and iron ions at evolving coronal mass ejection shocks propagating throughout the inner heliosphere from about 0.1 to 2.5 AU. An example shock with a compression ratio of 3.9 and a speed of 1000 km s{sup -1} close to the Sun is modeled using a two-dimensional MHD ZEUS code. The compression ratio and shock speed weakened to 2.3 and 630 km s{sup -1}, correspondingly, at 2 AU. Shocks with 15 Degree-Sign , 45 Degree-Sign , and 75 Degree-Sign angles between the upstream magnetic field and the shock normal were studied. The shock angle was kept constant throughout the simulation. Both gradual and impulsive events are studied. Diffusive shock acceleration is assumed at the shock and we use a total diffusion coefficient that includes a parallel diffusion coefficient which takes into account the upstream wave amplification, and a perpendicular diffusion coefficient which is based on the NonLinear Guiding Center theory. The transport of particles escaping from the shock is modeled using a Monte Carlo approach. Time-intensity profiles for protons and iron ions are obtained. We analyzed the radial dependence of peak fluxes (J) for both protons and iron ions from 0.5 to 2 AU. We find that the functional dependence is softer than R {sup -3} and is about R {sup -2.9} to R {sup -1.8} in the energy range of 0.3-5 MeV nuc{sup -1}. Quasi-perpendicular shock showed a steeper radial dependence than a quasi-parallel shock. Mixed events show a softer radial dependence at energies above 500 keV nuc{sup -1} for iron ions and above 1 MeV for protons. The values of J(R) depend on seed particle composition, particle energy, shock obliquity, and the interplanetary turbulence level. Consequently, we advocate using SEP event-specific computer modeling rather than empirical formulae for future forecasting of the radiation environment throughout the heliosphere during large SEP events.

79 ASTRONOMY AND ASTROPHYSICS↗

Dissipative axial inflation

We analyze in detail the background cosmological evolution of a scalar field coupled to a massless abelian gauge field through an axial term φ/ f {sub γ} F ∼ F , such as in the case of an axion. Gauge fields in this case are known to experience tachyonic growth and therefore can backreact on the background as an effective dissipation into radiation energy density ρ{sub R}, which can lead to inflation without the need of a flat potential. We analyze the system, for momenta k smaller than the cutoff f {sub γ}, including the backreaction numerically. We consider the evolution from a given static initial condition and explicitly show that, if f {sub γ} is smaller than the field excursion φ{sub 0} by about a factor of at least O (20), there is a friction effect which turns on before the field can fall down and which can then lead to a very long stage of inflation with a generic potential. In addition we find superimposed oscillations, which would get imprinted on any kind of perturbations, scalars and tensors. Such oscillations have a period of 4–5 efolds and an amplitude which is typically less than a few percent and decreases linearly with f {sub γ}. We also stress that the curvature perturbation on uniform density slices should be sensitive to slow-roll parameters related to ρ{sub R} rather than φ-dot {sup 2}/2 and we discuss the existence of friction terms acting on the perturbations, although we postpone a calculation of the power spectrum and of non-gaussianity to future work and we simply define and compute suitable slow roll parameters. Finally we stress that this scenario may be realized in the axion case, if the coupling 1/ f {sub γ} to U(1) (photons) is much larger than the coupling 1/ f {sub G} to non-abelian gauge fields (gluons), since the latter sets the range of the potential and therefore the maximal allowed φ{sub 0∼} f {sub G}.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance of modern color decompositions for standard candle LHC tree amplitudes

In the last decade, developments of matrix element and phase space generators have focused on providing good efficiency and maximal flexibility and automation for a wide range of physical processes. However, as recent studies have shown, they are a major bottleneck in the established Monte Carlo event generator toolchains. With the advent of the HL-LHC and ever rising precision requirements, future developments will need to focus on computational performance, especially at intermediate to large jet multiplicities. We present the novel BlockGen family of fast matrix element algorithms that are amenable for GPU acceleration, making use of modern, minimal color decompositions. Moreover, we discuss the performance achieved for standard candle processes such as V +jets and tt̄+jets production.

Bothmann, E. [Gottingen U.]↗

Machine learning-enabled computer vision for plant phenotyping: a primer on AI/ML and a case study on stomatal patterning

Abstract Artificial intelligence and machine learning (AI/ML) can be used to automatically analyze large image datasets. One valuable application of this approach is estimation of plant trait data contained within images. Here we review 39 papers that describe the development and/or application of such models for estimation of stomatal traits from epidermal micrographs. In doing so, we hope to provide plant biologists with a foundational understanding of AI/ML and summarize the current capabilities and limitations of published tools. While most models show human-level performance for stomatal density (SD) quantification at superhuman speed, they are often likely to be limited in how broadly they can be applied across phenotypic diversity associated with genetic, environmental, or developmental variation. Other models can make predictions across greater phenotypic diversity and/or additional stomatal/epidermal traits, but require significantly greater time investment to generate ground-truth data. We discuss the challenges and opportunities presented by AI/ML-enabled computer vision analysis, and make recommendations for future work to advance accelerated stomatal phenotyping.

Plant Sciences↗

Global Futures: a multithreaded execution model for Global Arrays-based applications

We present Global Futures (GF), an execution model extension to Global Arrays, which is based on a PGAS-compatible Active Message-based paradigm. We describe the design and implementation of Global Futures and illustrate its use in a computational chemistry application benchmark (Hartree-Fock matrix construction using the Self-Consistent Field method). Our results show how we used GF to increase the scalability of the Hartree-Fock matrix build to up to 6,144 cores of an Infiniband cluster. We also show how GF's multithreaded execution has comparable performance to the traditional process-based SPMD model.

99 GENERAL AND MISCELLANEOUS↗

The COVID-19 High-Performance Computing Consortium

In March of 2020, recognizing the potential of High Performance Computing (HPC) to accelerate understanding and the pace of scientific discovery in the fight to stop COVID-19, the HPC community assembled the largest collection of worldwide HPC resources to enable COVID-19 researchers worldwide to advance their critical efforts. Amazingly, the COVID-19 HPC Consortium was formed within one week through the joint effort of the Office of Science and Technology Policy (OSTP), the U.S. Department of Energy (DOE), the National Science Foundation (NSF), and IBM to create a unique public–private partnership between government, industry, and academic leaders. This article is the Consortium's story–how the Consortium was created, its founding members, what it provides, how it works, and its accomplishments. We will reflect on the lessons learned from the creation and operation of the Consortium and describe how the features of the Consortium could be sustained as a National Strategic Computing Reserve to ensure the nation is prepared for future crises.

97 MATHEMATICS AND COMPUTING↗

A modern concept of Lagrangian hydrodynamics

Here, we offer a modern interpretation of Lagrangian hydrodynamics as employed in Lagrangian simulations of compressible fluid flow. Our main result is to show that artificial viscosity, traditionally viewed as a numerical artifice to control unphysical oscillations in flows with shocks, actually represents a physical process and is necessary to derive accurate simulations in any compressible flow. We begin by reviewing the origins of two numerical devices, artificial viscosity and finite-volume methods. We proceed to construct a mathematical (PDE) model that incorporates those numerics and in which a new length scale, the observer, arises representing the discretization. Associated with that length scale, there are new inviscid fluxes that are the artificial viscosity as first formulated by Richtmyer and an artificial heat flux postulated by Noh but typically not included in Lagrangian codes. We discuss the connection of our results to bivelocity hydrodynamics. We conclude with some speculation as to the direction of future developments in multidimensional Lagrangian codes as computers get faster and have larger memories.

97 MATHEMATICS AND COMPUTING↗

SU-F-T-122: 4Dand 5D Proton Dose Evaluation with Monte Carlo

Purpose: We evaluated uncertainties in therapeutic proton doses of a lung treatment, taking into account intra-fractional geometry changes, such as breathing, and inter-fractional changes, such as tumor shrinkage and weight loss. Methods: A Monte Carlo study was performed using four dimensional CT image sets (4DCTs) and weekly repeat imaging (5DCTs) to compute fixed RBE (1.1) and variable RBE weighted dose in an actual lung treatment geometry. The MC2 Monte Carlo system was employed to simulate proton energy deposition and LET distributions according to a thoracic cancer treatment plan developed with a 3D-CT in a commercial treatment planning system, as well as in each of the phases of 4DCT sets which were recorded weekly throughout the course of the treatment. A cumulative dose distribution in relevant structures was computed and compared to the predictions of the treatment planning system. Results: Using the Monte Carlo method, dose deposition estimates with the lowest possible uncertainties were produced. Comparison with treatment planning predictions indicates that significant uncertainties may be associated with therapeutic lung dose prediction from treatment planning systems, depending on the magnitude of inter- and intra-fractional geometry changes. Conclusion: As this is just a case study, a more systematic investigation accounting for a cohort of patients is warranted; however, this is less practical because Monte Carlo simulations of such cases require enormous computational resources. Hence our study and any future case studies may serve as validation/benchmarking data for faster dose prediction engines, such as the track repeating algorithm, FDC.

60 APPLIED LIFE SCIENCES↗

Python in the NERSC Exascale Science Applications Program for Data

We describe a new effort at the National Energy Re- search Scientific Computing Center (NERSC) in performance analysis and optimization of scientific Python applications targeting the Intel Xeon Phi (Knights Landing, KNL) many- core architecture. The Python-centered work outlined here is part of a larger effort called the NERSC Exascale Science Applications Program (NESAP) for Data. NESAP for Data focuses on applications that process and analyze high-volume, high-velocity data sets from experimental/observational science (EOS) facilities supported by the US Department of Energy Office of Science. We present three case study applications from NESAP for Data that use Python. These codes vary in terms of “Python purity” from applications developed in pure Python to ones that use Python mainly as a convenience layer for scientists without expertise in lower level programming lan- guages like C, C++ or Fortran. The science case, requirements, constraints, algorithms, and initial performance optimizations for each code are discussed. Our goal with this paper is to contribute to the larger conversation around the role of Python in high-performance computing today and tomorrow, highlighting areas for future work and emerging best practices

Ronaghi, Zahra↗