Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heterogeneous 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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Union: An Automatic Workload Manager for Accelerating Network Simulation

With the rapid growth of the machine learning applications, the workloads of future HPC systems are anticipated to be a mix of scientific simulation, big data analytics, and machine learning applications. Simulation is a great research vehicle to understand the performance implications of co-running scientific applications with big data and machine learning workloads on large-scale systems. In this paper, we present Union, a workload manager that provides an automatic framework to facilitate hybrid workload simulation in CODES. Furthermore, we use Union, along with CODES, to investigate various hybrid workloads composed of traditional simulation applications and emerging learning applications on two dragonfly systems. The experiment results show that both message latency and communication time are important performance metrics to evaluate network interference. Network interference on HPC applications is more reflected by the message latency variation, whereas ML application performance depends more on the communication time.

heterogeneous workloads↗
DOE OSTI2022

Network Anomaly Detection Using Federated Learning

The internet is turning out to be an integral part of every-one's lives as more and more devices are being connected to serve societal needs. Our work is motivated by two ma-jor observations. Firstly, one drawback of connecting to the network is the threat of network attacks that can compromise users' private information, leading to data loss and adversely affecting productivity. There are several traditional security mechanisms to defend against these attacks, such as firewalls, virtual private networks (VPNs), demilitarized zones (DMZs), and vulnerability scanners. One way to prevent these attacks is early detection and prevention. However, these kinds of architecture do not scale very well because of their centralized nature. Secondly, we observe from heuristics and data set distributions that the majority of the requests made to a server are innocuous. Therefore, almost all server request data sets are highly imbalanced, weighted highly towards the harmless requests.

Marfo, William↗
DOE OSTI2022

Inverse Analysis with Variational Autoencoders: A Comparison of Shallow and Deep Networks

Inverse problems are applied to determine unknown properties by matching observational data with a physical model that often takes many parameters as input. To overcome the underconstrained nature of inverse problems and achieve good performance, an approach is presented involving regularization with a technique known as a variational autoencoder (VAE), which is trained to map a high-dimensional parameter space with a complex structure to a low-dimensional latent space with a simple structure. We apply this approach to unconditioned realizations of the parameters (heterogeneous hydraulic fields) for a hydrogeological inverse problem. Two types of hydraulic conductivity fields are used to evaluate the characterization for the different levels of heterogeneity complexity of the physical inputs. This approach keeps the computational cost of generating the training data low. The reason is unconditioned realizations neither rely on the observational data used to perform the inverse analysis nor require any groundwater flow model forward runs. In addition, this approach applies regularization on a low-dimensional latent space from the VAE and increases optimization efficiency through automatic differentiation. Furthermore, two different neural network (NN) structures are tested for their utility in using the VAE for inverse analysis. The performance of a deep, convolutional neural network strongly depends on the dimensionality of the latent space, which requires tuning. In contrast, a shallow, dense neural network provides consistently accurate characterization without tuning. Furthermore, our approach evaluates the advantages of the shallow, dense neural network over the deep, convolutional one and enables future application to a wide range of inverse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Probing catalytic activity in defect sites in transition metal oxides and sulfides using cluster models: A combined experimental and theoretical approach

Our research program has combined experimental and computational methods to study a range of cluster models for heterogeneous catalytic materials. The focus of our studies has been transition metal oxide and sulfide clusters in non-traditional oxidation states, and their chemical and physical interactions with water, to model steps in the electrocatalytic or photocatalytic hydrogen evolution reaction (HER). In the closing year of this project, we have also begun to explore addition systems in preparation for future grant proposals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Machine Learning Analysis of Hydrologic Exchange Flows and Transit Time Distributions in a Large Regulated River

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. High-resolution numerical models were often used to resolve the spatial and temporal variations of exchange flows, which are computationally expensive. In this study, we adopt Random Forest (RF) and Extreme Gradient Boosting (XGB) approaches for deriving reduced order models of hydrologic exchange flows and associated transit time distributions, with integrated field observations (e.g., bathymetry) and hydrodynamic simulation data (e.g., river velocity, depth). The setup allows an improved understanding of the influences of various physical, spatial, and temporal factors on the hydrologic exchange flows and transit times. The predictors also contain those derived using hybrid clustering, leveraging our previous work on river corridor system hydromorphic classification. The machine learning-based predictive models are developed and validated along the Columbia River Corridor, and the results show that the top parameters are the thickness of the top geological formation layer, the flow regime, river velocity, and river depth; the RF and XGB models can achieve 70% to 80% accuracy and therefore are effective alternatives to the computational demanding numerical models of exchange flows and transit time distributions. Each machine learning model with its favorable configuration and setup have been evaluated. The transferability of the models to other river reaches and larger scales, which mostly depends on data availability, is also discussed.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

Human capabilities: Long residence - An uncharted field

In the context of the planned NASA space station, it is anticipated that a number of problems will be posed by permanent residence in space. Not only will the intensive training which has characterized earlier manned programs diminish under the pressure of increased launch rates, but the new and larger crews will socially, educationally, psychologically and physiologically constitute a far more heterogeneous group than those of the past. In addition, advanced computer technology will result in a degree of reliance on onboard data processing which will alter the relations of ground and space operators, requiring careful integration of people and technologies. Attention is given to questions of food, robotics, extravehicular activity, hygiene and environmental noise.

Chambers, A.↗
NASA NTRS1990

Finite-element models of continental extension

Numerical models of the initial deformation of extending continental lithosphere, computed to investigate the control of preexisting thermal and mechanical heterogeneities on the style of deformation, are presented. The finite element method is used to calculate deformation with a viscoelastic-plastic model for the lithosphere. Comparisons of the results of analytic models and finite-element models using this method show that good results may be obtained by the numerical technique, even with elements containing both brittle and viscoelastic sampling points. It is shown that the gross style of initial extensional deformation is controlled by the depth and width of the initial heterogeneity which localizes deformation.

Lynch, H. David↗
DOE OSTI2013

Computational Mechanistic Studies of Acid-Catalyzed Lignin Model Dimers for Lignin Depolymerization

Lignin is a heterogeneous alkyl-aromatic polymer that constitutes up to 30% of plant cell walls, and is used for water transport, structure, and defense. The highly irregular and heterogeneous structure of lignin presents a major obstacle in the development of strategies for its deconstruction and upgrading. Here we present mechanistic studies of the acid-catalyzed cleavage of lignin aryl-ether linkages, combining both experimental studies and quantum chemical calculations. Quantum mechanical calculations provide a detailed interpretation of reaction mechanisms including possible intermediates and transition states. Solvent effects on the hydrolysis reactions were incorporated through the use of a conductor-like polarizable continuum model (CPCM) and with cluster models including explicit water molecules in the first solvation shell. Reaction pathways were computed for four lignin model dimers including 2-phenoxy-phenylethanol (PPE), 1-(para-hydroxyphenyl)-2-phenoxy-ethanol (HPPE), 2-phenoxy-phenyl-1,3-propanediol (PPPD), and 1-(para-hydroxyphenyl)-2-phenoxy-1,3-propanediol (HPPPD). Lignin model dimers with a para-hydroxyphenyl ether (HPPE and HPPPD) show substantial differences in reactivity relative to the phenyl ether compound (PPE and PPPD) which have been clarified theoretically and experimentally. The significance of these results for acid deconstruction of lignin in plant cell walls will be discussed.

09 BIOMASS FUELS↗
DOE OSTI2011

A New Approach to Computing Normal Tissue Complication Probability of an Intensity-Modulated Radiotherapy Treatment with Stereotactic Radiotherapy Boost of Nasopharyngeal Carcinoma: A Case Study

We attempted to develop a method to compute the normal tissue complication probability (NTCP) of various critical organs from combined intensity-modulated radiotherapy (IMRT) and stereotactic radiotherapy (SRT) boost treatment of nasopharyngeal carcinoma (NPC), with the aid of a nonlinear image registration method. The SRT's planning computed tomography (CT) of a NPC patient treated with IMRT was warped to the IMRT's planning CT using a nonlinear image registration. Because CT and dose were inherently in-register, the entire dose distribution could be deformed using the same deformation field derived from the two CT sets. Using the biologically effective dose concept and the linear-quadratic model, physical doses of IMRT and SRT were converted to a 2Gy-per-fraction equivalent dose to facilitate dose summation. The variation of organs' maximum doses of the combined treatments between traditional maximum dose sum and the proposed method was 1.5Gy {+-} 1.7Gy. After the correction of the effect of fractionation and dose heterogeneity within each organ, NTCP of each organ of interest was computed for the combined treatments. Based on the results of this case study, it is believed that dose registration could be a method for the NTCP computation of various critical organs when different fractionation schemes of radiation therapy treatment are instituted.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

ARES v1.x - Performance Portable Tool to Simulate Supernovae based on Parthenon Framework

Historically, codes for simulating supernovae (such as Arepo, FLASH or LEAFS) have been at the forefront of scientific high-performance computing to the immense computational resources required for full 3D simulations. However, given the shift towards heterogenous HPC architectures, many current-generation codes are at the risk of losing their competitiveness as they are only designed to run on homogeneous CPU-only systems. There exist several efforts to enable these codes for GPU’s, however, these efforts only consider specific architectures or vendors (e.g., implement only CUDA or HIP), limiting themselves to a small range of exascale computing systems. Frameworks such as Kokkos aim to provide a framework which is agnostic of the targeted architecture, enabling the development of performant and portable code. In the Ares code, we develop a performance portable tool to simulate supernovae based on the Parthenon Framework, which in turn uses Kokkos in the background. Here, the Parthenon Framework provides an interface to the underlying mesh-refinement routines, which form the backbone of our code. In addition, we incorporate the already existing Singularity-EOS toolkit to provide us with various equations of state, primarily the Helmholtz equation of state. We also include the JINA Reaclib as a basis for our nuclear network solver. Finally, we implement a gravity solver to complete the required physics. This setup will provide us with a minimal code base to simulate supernova in a similar style to the tried-and-tested Arepo code, but in a futureproof performance portable framework.

Lim, Hyun↗
DOE OSTI2016

A moment-convergence method for stochastic analysis of biochemical reaction networks

Traditional moment-closure methods need to assume that high-order cumulants of a probability distribution approximate to zero. However, this strong assumption is not satisfied for many biochemical reaction networks. Here, we introduce convergent moments (defined in mathematics as the coefficients in the Taylor expansion of the probability-generating function at some point) to overcome this drawback of the moment-closure methods. As such, we develop a new analysis method for stochastic chemical kinetics. This method provides an accurate approximation for the master probability equation (MPE). In particular, the connection between low-order convergent moments and rate constants can be more easily derived in terms of explicit and analytical forms, allowing insights that would be difficult to obtain through direct simulation or manipulation of the MPE. In addition, it provides an accurate and efficient way to compute steady-state or transient probability distribution, avoiding the algorithmic difficulty associated with stiffness of the MPE due to large differences in sizes of rate constants. Applications of the method to several systems reveal nontrivial stochastic mechanisms of gene expression dynamics, e.g., intrinsic fluctuations can induce transient bimodality and amplify transient signals, and slow switching between promoter states can increase fluctuations in spatially heterogeneous signals. The overall approach has broad applications in modeling, analysis, and computation of complex biochemical networks with intrinsic noise.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Mapping the mechanical properties of rocks using automated microindentation tests

A microindentation scanner is constructed that measures the spatial fluctuation in the elastic properties of natural rocks. This novel instrument performs automated indentation tests on the surface of a rock slab and outputs 2-D maps of the indentation modulus at submillimeter resolution. Maps obtained for clean, well-consolidated, sandstone are presented and demonstrate the capabilities of the instrument. We postulate that the elastic structure of sandstones correlates well with their visual appearance. Moreover, we observe that the probability distribution of the indentation modulus fluctuations across the slab surfaces can be modeled using a lognormal probability density function. To illustrate possible use of the data obtained with the microindentation scanner, we use roughly 10 cm × 10 cm scans with millimeter resolution over four sandstone planar slabs to numerically compute the overall drained elastic moduli for each sandstone sample. We demonstrate that such numerically computed moduli are well modeled using the multicomponent form of the Hashin-Shtrikman lower bound that employs the observed lognormal probability distribution for the mesoscopic-scale moduli (the geometric mean works almost the same). We also compute the seismic attenuation versus frequency associated with wave-induced fluid flow between the heterogeneities in the scanned sandstones and report relatively small values for the inverse quality factor (Q -1 <10 -2 ) in the seismic frequency band 10 2 Hz 4 Hz. The numerically computed frequency dependence in the attenuation varies from one type of sandstone to another, and we determine vital anisotropy in the attenuation associated with waves propagating in different directions.

58 GEOSCIENCES↗
DOE OSTI2017

To address surface reaction network complexity using scaling relations machine learning and DFT calculations

Surface reaction networks involving hydrocarbons exhibit enormous complexity with thousands of species and reactions for all but the very simplest of chemistries. We present a framework for optimization under uncertainty for heterogeneous catalysis reaction networks using surrogate models that are trained on the fly. The surrogate model is constructed by teaching a Gaussian process adsorption energies based on group additivity fingerprints, combined with transition-state scaling relations and a simple classifier for determining the rate-limiting step. The surrogate model is iteratively used to predict the most important reaction step to be calculated explicitly with computationally demanding electronic structure theory. Applying these methods to the reaction of syngas on rhodium(111), we identify the most likely reaction mechanism. Lastly, propagating uncertainty throughout this process yields the likelihood that the final mechanism is complete given measurements on only a subset of the entire network and uncertainty in the underlying density functional theory calculations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Modeling Frequency–Independent Q Viscoacoustic Wave Propagation in Heterogeneous Media

Quantifying the attenuation of seismic waves propagating in the Earth interior is critical to study the subsurface structure. Previous studies have proposed fractional anelastic wave equations to model the frequency–independent Q seismic wave propagation. Such wave equations involve fractional derivatives that pose computational challenges for the numerical schemes in terms of accuracy and efficiency when dealing with heterogeneous Earth media. To tackle these challenges, here we derive a new viscoacoustic wave equation, where the power terms of the fractional Laplacian operators are spatially independent, thus accurate and efficient methods (e.g., the Fourier pseudospectral method) can be adopted. Our derivation enables the resultant equation to capture both amplitude and phase signatures of the anelastic wave propagation by matching the complex wave numbers for all the frequencies of interest. We verify the derivation by comparing the dispersion curves of both the attenuation factor and the phase velocity produced by the new wave equation with their theoretical values as well as the Pierre Shale in situ measurements. Following that, we use a synthetic attenuating gas chimney model to demonstrate the attenuation effects on seismic waveforms and then construct the Q–compensated reverse time migration to undo these effects for seismic image enhancement. Finally, we find that our forward modeling results can characterize the spatiotemporal attenuation effects revealed in the Frio–II CO2 injection time–lapse seismic monitoring data. Here, we expect this proposed equation to be useful to quantify the attenuation in seismic data to push the resolution limits of seismic imaging and inversion.

25 ENERGY STORAGE↗
DOE OSTI2016

A fast marching algorithm for the factored eikonal equation

The eikonal equation is instrumental in many applications in several fields ranging from computer vision to geoscience. This equation can be efficiently solved using the iterative Fast Sweeping (FS) methods and the direct Fast Marching (FM) methods. However, when used for a point source, the original eikonal equation is known to yield inaccurate numerical solutions, because of a singularity at the source. In this case, the factored eikonal equation is often preferred, and is known to yield a more accurate numerical solution. One application that requires the solution of the eikonal equation for point sources is travel time tomography. This inverse problem may be formulated using the eikonal equation as a forward problem. While this problem has been solved using FS in the past, the more recent choice for applying it involves FM methods because of the efficiency in which sensitivities can be obtained using them. However, while several FS methods are available for solving the factored equation, the FM method is available only for the original eikonal equation. In this paper we develop a Fast Marching algorithm for the factored eikonal equation, using both first and second order finite-difference schemes. Our algorithm follows the same lines as the original FM algorithm and requires the same computational effort. In addition, we show how to obtain sensitivities using this FM method and apply travel time tomography, formulated as an inverse factored eikonal equation. Numerical results in two and three dimensions show that our algorithm solves the factored eikonal equation efficiently, and demonstrate the achieved accuracy for computing the travel time. We also demonstrate a recovery of a 2D and 3D heterogeneous medium by travel time tomography using the eikonal equation for forward modeling and inversion by Gauss–Newton.

97 MATHEMATICS AND COMPUTING↗