Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallelization”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Online multimedia retrieval on CPU–GPU platforms with adaptive work partition

Nearest neighbors search is a core operation found in several online multimedia services. These services have to handle very large databases, while, at the same time, they must minimize the query response times observed by users. This is specially complex because those services deal with fluctuating query workloads (rates). Consequently, they must adapt at run-time to minimize the response times as the load varies. In this paper, we address the aforementioned challenges with a distributed memory parallelization of the product quantization nearest neighbor search, also known as IVFADC, for hybrid CPU–GPU machines. Overall, our parallel IVFADC implements an out-of-GPU memory execution scheme to use the GPU for databases in which the index does not fit in its memory, which is crucial for searching in very large databases. The careful use of CPU and GPU with work stealing led to an average response time reduction of 2.4 as compared to using the GPU only. Also, our approach to adapt the system to fluctuating loads, called Dynamic Query Processing Policy (DQPP), attained a response time reduction of up to 5 vs. the best static (BS) policy for moderate loads. The system has attained high query processing rates and near-linear scalability in all experiments. We have evaluated our system on a machine with up to 256 NVIDIA V100 GPUs processing a database of 256 billion SIFT features vectors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Exploring thermal runaway propagation in Li-ion batteries through high-speed X-ray imaging and thermal analysis: Impact of cell chemistry and electrical connections

Battery safety design is important to consider from the individual Li-ion cell to the level of the macro-system. On the macro-level, failure in one single cell can lead to propagation of the thermal runaway and rapidly set a whole battery pack on fire. Factors that can impact the propagation outcome, such as cell model/chemistry and electrical connection are here investigated using a combination of measurements. Several abusive tests were conducted, combining two different cell models (Molicel P42A and LG M50, both 21700s) in series and parallel connections (16 tests per configuration). Overall, a propagation outcome of 56% was measured from the 32 conducted tests, a minimum temperature of 150 °C was required to initiate propagation, and the fastest propagation occurred in 123 s. Temperature measurements were higher in series connected cells, initiating the discussion of cell chemistry and internal resistance on this effect. The difference in current-flow during thermal runaway in series and parallel connections, and how this can affect the temperature evolution is further discussed. Spatio-temporal mapping of X-ray radiography allowed us to derive the speed of thermal runaway evolution inside the battery and has shown that series connected cells, in particular P42A, occur faster. It was further observed that deviant sidewall behaviors such as temperature-induced breaches and pressure-induced ruptures occurred in P42As only respective nail-penetrated cells only.

25 ENERGY STORAGE↗
DOE OSTI2025

Microstructure and CPO evolution of dynamically recrystallized olivine during complex deformation conditions: a full-field numerical modeling approach

The rheological properties of mantle rocks are strongly dependent on their crystallographic preferred orientation (CPO). Olivine CPO, defined by the orientation of seismically fast [100] axes parallel to flow direction, is also thought to be a dominant contributor to seismic anisotropy in the Earth's upper mantle. However, the amount of deformation needed to overprint a new CPO on a pre-existing fabric and the impact of the inherited CPOs on the transient microstructure evolution, remain unknown. This study employs a full-field numerical approach (VPFFT-ELLE) to explore the dynamic recrystallization and microstructural evolution of olivine polycrystalline aggregates under complex deformation conditions. We test four combinations of successive pure shear and simple shear boundary conditions. Findings indicate that inherited CPOs influence subsequent deformation in a manner dependent on the kinematic relationship between successive stages. In all cases, a minor strain increment (ε ∼0.3–0.6) is sufficient to erase the previous microstructure and CPO. However, when deformation conditions change dramatically (e.g., stretching direction changes orthogonally), the intensity of the new CPO developed is significantly lower and strain distribution are specially altered. During a transient strain stage, pre-existing microstructures undergo extensive reworking, especially when deformation conditions are changed dramatically, such as switching from simple shear to a pure shear condition with a parallel shortening direction relative to the stretching direction. We estimate the significance of these results in interpreting observations of seismic velocity anisotropy, concluding that P-wave seismic anisotropy is significantly and positively correlated with the evolution of olivine CPO with deformation history. This research underscores the transient nature of microstructural rearrangement in olivine aggregates and the necessity for caution in interpreting seismic anisotropy in regions with complex deformation histories, as inherited CPOs can influence current fabric development and induce deviation to the present deformation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Deformation behavior of nanoscale Al–Al 2 Cu eutectics studied by in situ micropillar compression

We report deformation behavior of nanoscale laser processed Al–Al 2 Cu eutectics at room temperature is characterized through in situ micro-pillar compression testing in a scanning transmission microscope. Interlamellar spacing of Al–Al 2 Cu eutectics varies from hundreds of nanometers to 20 nm. Three different sizes of micro-pillars are fabricated in order to study the deformation behaviors of single colony and multiple colonies, corresponding to the single crystal and polycrystal respectively. For single colonies, lamellar orientations parallel, normal or inclined to the loading direction were tested. The main findings are: 1) the plasticity mechanisms strongly depend on loading orientation: buckling and kinking in the parallel-loaded eutectics, planar sliding along Al–Al 2 Cu lamellar interfaces in the incline-loaded eutectics and localized shearing in the normal-loaded eutectics. 2) the incline-loaded eutectics exhibits the lowest compression flow strength, and the normal-loaded eutectics has the highest compression flow strength. 3) with decreasing inter-lamellar spacing, the strength increases and plasticity is uniformly distributed, as opposed to shear localization. Highest compressive plasticity is observed in polycrystalline eutectics with an ensemble of lamellar orientations with ~20 nm average spacing: 17.9% at flow stress of 1.63 GPa, and degenerate, bimodal morphology: 11.1% at flow stress of 1.36 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Self-equilibrated backstresses induce compensation between hardening and softening: Micromechanical and microstructural features

For crystal plasticity formulation neglecting the strain gradients at the dislocation scale, backstresses are necessary to satisfy the mesoscopic stress equilibrium and strain compatibility requirements. Backstresses in a material point and those between neighboring points are treated separately. Physically backstresses can dynamically induce either hardening or softening effects to adjust the dislocation hardening of slip systems, whereas the sum of the stresses resulting in backstresses in the representative volume element (RVE) is zero. Backstress hardening suppresses certain slip activities, while backstress softening triggers excess dislocations. A two-phase brass is studied. Backstresses make the hard phase deform smaller and the soft phase deform more than those without employing backstresses in a model. In the former, more homogeneous lattice strain and higher yield stress are obtained. Backstress has limited effects on the evolution of the critical resolved shear stress and crystallographic texture. Furthermore, it is a major cause of the banded structures featuring sharp crystallographic contrasts which are not parallel to the slip planes. Tests indicate that backstress induces a serrated behavior at the early stages of deformation and induces the geometrically necessary bands (GNBs), which are parallel to the macroscopic axes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Anisotropy evolution of elastomeric foams during uniaxial compression measured via in-situ X-ray computed tomography

We have characterized the three-dimensional evolution of microstructural anisotropy of a family of elastomeric foams during uniaxial compression via in-situ X-ray computed tomography. Flexible polyurethane foam specimens with densities of 136, 160 and 240 kg/m 3 were compressed in uniaxial stress tests both parallel and perpendicular to the foam rise direction, to engineering strains exceeding 70%. The uncompressed microstructures show slightly elongated ellipsoidal pores, with elongation aligned parallel to the foam rise direction. The evolution of this microstructural anisotropy during deformation is quantified based on the autocorrelation of the image intensity, and verified via the mean intercept length as well as the shape of individual pores. Trends are consistent across all three methods. In the rise direction, the material remains transversely anisotropic throughout compression. Anisotropy initially decreases with compression, reaches a minimum, then increases up to large strains, followed by a small decrease in anisotropy at the largest strains as pores collapse. Compression perpendicular to the foam rise direction induces secondary anisotropy with respect to the compression axis, in addition to primary anisotropy associated with the foam rise direction. In contrast to compression in the rise direction, primary anisotropy initially increases with compression, and shows a slight decrease at large strains. These surprising non-monotonic trends and qualitative differences in rise and transverse loading are explained based on the compression of initially ellipsoidal pores. Microstructural anisotropy trends reflect macroscopic stress-strain and lateral strain response. These findings provide novel quantitative connections between three-dimensional microstructure and anisotropy in moderate density polymer foams up to large deformation, with important implications for understanding complex three-dimensional states of deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Utilizing ink composition to tune bulk-electrode gas transport, performance, and operational robustness for a Fe-N-C catalyst in polymer electrolyte fuel cell

With lower site density and turnover frequency, platinum group metal (PGM)-free catalysts based electrodes are often greater than 50 mu m thick in order to increase performance across the fuel cell operating range. Consequently, PGM-free electrodes have an additional bulk electrode transport resistance beyond the local or aggregate level transport in thin platinum-based electrodes. In parallel to the development of more active and durable PGM-free catalysts, advancements in understanding the interplay between PGM-free electrode fabrication, bulk-electrode transport, proton conductivity and performance are needed. Here, the relationship between ionic and gas phase transport through the electrode thickness is modified by adjusting electrocatalyst and ionomer flocculation/interaction at the ink level. The influence of the ink composition (water/n-propanol content) is examined via various in-situ electrochemical and ex-situ characterization techniques and the resulting electrode structure/performance relationship contrasted with electrode performance robustness across a range of relative humidity (RH). For the electrocatalyst examined here, a water-rich (82 wt% H2O) ink formulation was favorable for operation at high RH due to improved molecular diffusion through larger electrode pores. In contrast, the improved interactions between ionomer and electrocatalyst enabled a more robust electrode and higher performance during low RH operation for the 50 wt% H2O content ink.

08 HYDROGEN↗
DOE OSTI2022

An overview of how parity-violating electron scattering experiments are performed at CEBAF

Parity-violating electron-scattering experiments represent an important focus of the nuclear physics experimental program at the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab. These experiments pose significant challenges because the scattering asymmetries can be very small, of the order parts-per-million and smaller. To succeed, the properties of the electron beam such as current, position, size and energy, must be very nearly identical in the two electron-polarization spin states (parallel and anti-parallel relative to the direction of beam motion at the scattering target). This paper describes the origins of unwanted helicity-correlated beam asymmetries present on the electron beam and methods to minimize them to acceptable levels.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Alignment of the CLAS12 central hybrid tracker with a Kalman Filter

Several factors can contribute to the difficulty of aligning the sensors of tracking detectors, including a large number of modules, multiple types of detector technologies, and non-linear strip patterns on the sensors. All three of these factors apply to the CLAS12 CVT, which is a hybrid detector consisting of planar silicon sensors with non-parallel strips, and cylindrical micromegas sensors with longitudinal and arc-shaped strips located within a 5 T superconducting solenoid. To align this detector, we used the Kalman Alignment Algorithm, which accounts for correlations between the alignment parameters without requiring the time-consuming inversion of large matrices. This is the first time that this algorithm has been adapted for use with hybrid technologies, non-parallel strips, and curved sensors. We present the results for the first alignment of the CLAS12 CVT using straight tracks from cosmic rays and from a target with the magnetic field turned off. After running this procedure, we achieved alignment at the level of 10 µm, and the widths of the residual spectra were greatly reduced. These results attest to the flexibility of this algorithm and its applicability to future use in the CLAS12 CVT and other hybrid or curved trackers, such as those proposed for the future Electron-Ion Collider.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Towards scaling community detection on distributed-memory heterogeneous systems

Distributed multi-GPU systems pose significant challenges and opportunities for efficient execution of parallel applications. Graph algorithms are generally characterized by irregular memory accesses, low computation to communication ratios, and load balancing problems that are especially hard to address on multi-GPU systems. Graph community detection is an important problem in the emerging domain of graph analytics with numerous applications. In this paper, we present our ongoing work on distributed-memory multi-GPU implementation for graph community detection. Our work parallelizes the widely used (albeit serial) Louvain method on distributed multi-GPU platforms. Supported by an extensive set of experiments on a multi-GPU enabled supercomputer (OLCF Summit) and a single compute node (Nvidia DGX-2®), we demonstrate competitive performance to existing distributed-memory CPU-based implementation, and up to 6.5 better results than Nvidia RAPIDS® CUGRAPH. To the best of our knowledge, this work represents the first effort for community detection on distributed multi-GPU systems. Our approach and related findings can be extended to numerous other iterative graph algorithms on multi-GPU systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Perpendicular in-plane negative magnetoresistance in ZrTe 5

The unique band structure in topological materials frequently results in unusual magneto-transport phenomena, one of which is in-plane longitudinal negative magnetoresistance (NMR) with the magnetic field aligned parallel to the electrical current direction. This NMR is widely considered as a hallmark of chiral anomaly in topological materials. Here we report the observation of in-plane NMR in the topological material ZrTe 5 when the in-plane magnetic field is both parallel and perpendicular to the current direction, revealing an unusual case of quantum transport beyond the chiral anomaly. We find that a general theoretical model, which considers the combined effect of Berry curvature and orbital moment, can quantitatively explain this in-plane NMR. In conclusion, our results provide new insights into the understanding of in-plane NMR in topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

The effects of thickness, polarization, and strain on vibrational modes of 2D Fe 3 GeTe 2

In this study, we investigated the effects of thickness, light polarization, and strain on the Raman spectra of two-dimensional (2D) Fe 3 GeTe 2 (FGT) crystals synthesized via chemical vapor transport. The crystals are thoroughly characterized using a combination of microscopic, diffraction, and spectroscopic techniques. Particularly, a systematic angle-resolved polarized Raman spectroscopy study reveals a clear polarization dependence of the Raman intensity in both parallel and crossed polarization directions, with the $A$$^{1}_{g}$ mode completely disappearing in the crossed polarization direction. The angle-dependent intensity of both the $A$$^{1}_{g}$ and $E$$^{2}_{2g}$ modes in parallel polarization and the intensity of the $E$$^{2}_{2g}$ mode in the crossed polarization remain constant at all angles. These findings align with predictions from Raman tensor analysis, providing compelling evidence for the unambiguous assignment of the $A$$^{1}_{g}$ and $E$$^{2}_{2g}$ modes to specific peaks observed in the Raman spectrum of FGT, resolving existing confusion in the literature regarding their assignment. Furthermore, we examine the effect of strain on the Raman spectrum of 2D FGT in-situ using a bending device. Our study, conducted on a monolayer to few-layer 2D FGT deposited onto polyethylene terephthalate and subjected to outward (inward), i.e., tensile (compressive) bending, demonstrates appreciable downshifting (upshifting) of the Raman peak position of both $A$$^{1}_{g}$ and $E$$^{2}_{2g}$ modes. Furthermore, these findings are particularly significant given that strain engineering represents an effective approach to modulate the magnetic properties of FGT and other 2D magnetic materials.

2D materials↗
DOE OSTI2024

Utilizing waste heat in wastewater treatment plants for water desalination: Modeling and Multi-Objective optimization of a Multi-Effect desalination system using Decision Tree Regression and Pelican optimization algorithm

This paper examines the feasibility of using waste heat from wastewater treatment plants (WWTPs) for water desalination. A model was developed to utilize waste heat from the gensets at As Samra WWTP in Jordan, using real data and TRNSYS® software to calculate available waste heat. The desalination process was then modeled with ASPEN PLUS® software, focusing on multi-effect desalination (MED). Both series and parallel configurations for the MED system were compared. The study investigated the effects of system feeding flow rate, feeding pressure, and heat input on productivity, performance ratio, and recovery ratio. The study also introduces a novel optimization technique combining machine learning and modern optimization algorithms to maximize system productivity and performance. Initially, a decision tree regression (DTR) model is developed to establish relationships between key independent variables (flow rate, feed pressure, and heat input) and dependent variables (productivity, performance ratio, and recovery ratio). The Pelican Optimization Algorithm (POA) is then used to identify the optimal values of the independent variables for maximum productivity and performance. The results show that using a series configuration yields a system productivity of 3984.2 kg/hr, a performance ratio of 3.78, and a recovery ratio of 0.991 at a feed flow rate of 4000 kg/hr, feed pressure of 3 bars, and heat input of 719 kW. Optimal productivity (4421 kg/hr), performance ratio (3.81), and recovery ratio (0.851) are achieved at a feed flow rate of 5166 kg/hr, feed pressure of 3.2 bars, and heat input of 794 kW. In conclusion, the techno-economic assessment indicates a levelized cost of water of 1.63 USD/m 3 for parallel configurations and 1.65 USD/m 3 for series configurations, with a payback period of less than two years.

42 ENGINEERING↗
DOE OSTI2026

Experimental and simulation study of target biasing effects on plasma transport in linear plasma device MPS-LD

Linear plasma devices (LPDs) are important experimental platforms for investigating plasma–material interactions (PMI). In PMI experiments, it has been found that applying a target bias not only effectively modifies the incident ion energy, but also induces significant changes in the electron density and electron temperature, whereby the evolution of these plasma parameters is primarily governed by plasma transport processes. However, at present, the physical process and mechanism underlying such bias-induced variations remain unclear. In this work, biasing experiments under argon plasma discharge conditions were first carried out on the MPS-LD device. For the corresponding experiments, an electric potential model was newly developed based on the BOUT++ LPD module, enabling self-consistent simulations of plasma transport under biased conditions. Numerical simulations were then performed to reproduce the experimental results and to validate the accuracy of the proposed model. Finally, by combining experimental measurements with numerical simulations, a bias-voltage scan was performed to investigate how the electron density and electron temperature vary with the bias voltage (U bias ). The results show that applying negative bias decreases the target electron density (n e,T ) while increasing the target electron temperature (T e,T ). In contrast, positive bias increases both n e,T and T e,T ; however, at high positive bias, n e,T first reaches a maximum and subsequently decreases with further increases in U bias . The underlying physical mechanisms are analyzed using particle flux, momentum, and energy conservation. It indicates that the applied bias regulates the parallel electric field, thereby changing ion and electron velocities, and consequently affecting the electron density. At high positive bias, the ion velocity is further influenced by ion viscosity, leading to the reversal in n e,T . Meanwhile, the enhanced parallel electric field drives stronger currents, significantly increasing ion–electron frictional work and converting the input bias power into electron energy, which raises the electron temperature. In conclusion, these results contribute to a deeper understanding of the effects and mechanisms of biasing on plasma transport in the MPS-LD device.

BOUT++ simulation↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2024

Bifurcation of equilibrium positions for ellipsoidal particles in inertial shear flows between two walls

We conducted a systematic numerical investigation of spherical, prolate and oblate particles in an inertial shear flow between two parallel walls, using smoothed particle hydrodynamics (SPH). It was previously shown that above a critical Reynolds number, spherical particles experience a supercritical pitchfork bifurcation of the equilibrium position in shear flow between two parallel walls, namely that the central equilibrium position becomes unstable, leading to the emergence of two new off-centre stable positions (Foxet al.,J. Fluid Mech., vol. 915, 2021). This phenomenon was unexpected given the symmetry of the system. In addition to confirming this finding, we found, surprisingly, that ellipsoidal particles can also return to the centre position from the off-centre positions when the particle Reynolds number is further increased, while spherical particles become unstable under this increased Reynolds number. By utilizing both SPH and the finite element method for flow visualization, we explained the underlining mechanism of this reverse of bifurcation by altered streamwise vorticity and symmetry breaking of pressure. Furthermore, we expanded our investigation to include asymmetric particles, a novel aspect that had not been previously modelled, and we observed similar trends in particle dynamics for both symmetric and asymmetric ellipsoidal particles. While further validation through laboratory experiments is necessary, our research paves the road for development of new focusing and separation methods for shaped particles.

Mechanics↗