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,279 records · Page 71

Mechanism Analysis of Wind Turbine Var Oscillations

Electromagnetic transient simulation of parallel connected 4-MW type-3 wind turbines based on original equipment manufacturer's real-code turbine model shows 1.2-Hz turbine-turbine oscillations in reactive power. This letter reveals why such oscillations occur in the individual var measurement, while being insignificant in the total var measurement, regardless of the varying grid impedance. We adopt two analysis approaches: open-loop single-input single-output analysis and network decomposition. The two approaches differ in their treatment of turbine-network interaction. The open-loop analysis shows that the turbine-turbine oscillation mode is due to an open-loop system pole being attracted to an open-loop system zero. Furthermore, we use network decomposition method to explain why this mode is observable in individual vars while not observable in the total var. The entire system of n -turbines can be viewed as n decoupled circuits. For the two-turbine case, the system has an aggregated mode and a turbine-turbine oscillation mode. Here, the aggregated mode is associated with a circuit associated with the total var, while the turbine-turbine oscillation mode is associated with the var difference and is insensitive to the grid parameters.

17 WIND ENERGY↗

True Load Balancing for Matricized Tensor Times Khatri-Rao Product

MTTKRP is the bottleneck operation in algorithms used to compute the CP tensor decomposition. For sparse tensors, utilizing the compressed sparse fibers (CSF) storage format and the CSF-oriented MTTKRP algorithms is important for both memory and computational efficiency on distributed-memory architectures. Existing intelligent tensor partitioning models assume the computational cost of MTTKRP to be proportional to the total number of nonzeros in the tensor. However, this is not the case for the CSF-oriented MTTKRP on distributed-memory architectures. We outline two deficiencies of nonzero-based intelligent partitioning models when CSF-oriented MTTKRP operations are performed locally: failure to encode processors' computational loads and increase in total computation due to fiber fragmentation. We focus on existing fine-grain hypergraph model and propose a novel vertex weighting scheme that enables this model encode correct computational loads of processors. We also propose to augment the fine-grain model by fiber nets for reducing the increase in total computational load via minimizing fiber fragmentation. In this way, the proposed model encodes minimizing the load of the bottleneck processor. In conclusion, parallel experiments with real-world sparse tensors on up to 1024 processors prove the validity of the outlined deficiencies and demonstrate the merit of our proposed improvements in terms of parallel runtimes.

97 MATHEMATICS AND COMPUTING↗

Design and Optimization of a 200-kW Medium-Frequency Transformer for Medium-Voltage SiC PV Inverters

This paper presents a design and optimization methodology for a 200-kW medium-frequency transformer (MFT) based on low-loss magnetic core (FINEMET FT-3TL). The proposed optimal design methodology consists of predesign, preliminary design, and optimal design. In the preliminary design, the parallel-concentric winding structure is selected to increase the current carrying capability and reduce the leakage inductance. Based on the parallel-concentric winding concept, a novel cooing and insulation structure with 3-D printed bobbins are proposed. The optimal designed MFT prototype achieves a power density higher than 19.23 kW/L. The electrical insulation system is tested at 12 kV ac peak voltage. In addition, the partial discharge (PD) test is conducted at 7.5 kV ac peak voltage to ensure the PD-free design. The high-frequency bipolar pulsewidth modulation voltage insulation (PD) test is the first time applied in MFT design and test process. Finally, the transformer is applied in a dual-active-bridge (DAB) converter with 200 kW rated power. The peak efficiency of the DAB converter is 99.53%, and the efficiency at 200 kW is 98.85%. The peak efficiency of MFT is 99.844%, and the efficiency at 200 kW is 99.842%.

42 ENGINEERING↗

Scalable Comparative Visualization of Ensembles of Call Graphs

Optimizing the performance of large-scale parallel codes is critical for efficient utilization of computing resources. Code developers often explore various execution parameters, such as hardware configurations, system software choices, and application parameters, and are interested in detecting and understanding bottlenecks in different executions. They often collect hierarchical performance profiles represented as call graphs, which combine performance metrics with their execution contexts. The crucial task of exploring multiple call graphs together is tedious and challenging because of the many structural differences in the execution contexts and significant variability in the collected performance metrics (e.g., execution runtime). In this paper, we present Ensemble CallFlow to support the exploration of ensembles of call graphs using new types of visualizations, analysis, graph operations, and features. We introduce ensemble-Sankey , a new visual design that combines the strengths of resource-flow (Sankey) and box-plot visualization techniques. Whereas the resource-flow visualization can easily and intuitively describe the graphical nature of the call graph, the box plots overlaid on the nodes of Sankey convey the performance variability within the ensemble. Our interactive visual interface provides linked views to help explore ensembles of call graphs, e.g., by facilitating the analysis of structural differences, and identifying similar or distinct call graphs. Finally, we demonstrate the effectiveness and usefulness of our design through case studies on large-scale parallel codes.

97 MATHEMATICS AND COMPUTING↗

Efficient Probabilistic Visualization of Local Divergence of 2D Vector Fields with Independent Gaussian Uncertainty

This work focuses on visualizing uncertainty of local divergence of two-dimensional vector fields. Divergence is one of the fundamental attributes of fluid flows, as it can help domain scientists analyze potential positions of sources (positive divergence) and sinks (negative divergence) in the flow. However, uncertainty inherent in vector field data can lead to erroneous divergence computations, adversely impacting downstream analysis. While Monte Carlo (MC) sampling is a classical approach for estimating divergence uncertainty, it suffers from slow convergence and poor scalability with increasing data size and sample counts. Thus, we present a two-fold contribution that tackles the challenges of slow convergence and limited scalability of the MC approach. (1) We derive a closed-form approach for highly efficient and accurate uncertainty visualization of local divergence, assuming independently Gaussian-distributed vector uncertainties. (2) We further integrate our approach into Viskores, a platform-portable parallel library, to accelerate uncertainty visualization. In our results, we demonstrate significantly enhanced efficiency and accuracy of our serial analytical (speed-up up to 1946×) and parallel Viskores (speed-up up to 19698×) algorithms over the classical serial MC approach. We also demonstrate qualitative improvements of our probabilistic divergence visualizations over traditional mean-field visualization, which disregards uncertainty. We validate the accuracy and efficiency of our methods on wind forecast and ocean simulation datasets.

Ouermi, Timbwaoga [University of Utah]↗

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long‐standing interaction among land plants and fungi

SUMMARY The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant–fungal interactions today, it is hypothesized these kingdoms have a long‐standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free‐living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia , interact with P. patens in coculture. We also assess how Mollicute ‐related or Burkholderia ‐related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high‐throughput phenomics, microscopy, RNA‐sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high‐throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant–fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens .

59 BASIC BIOLOGICAL SCIENCES↗

DESIGN AND TESTING OF DIMES CARBON ABLATION RODS IN THE DIII-D TOKAMAK

We present the design of ATJ graphite rods developed for ablation experiments under high heat flux (up to 50 MW/m2) in the lower divertor of the DIII-D tokamak, a magnetic plasma confinement device. This work is motivated by the need to test ablation models relevant to carbon-based thermal shields used in high-speed spacecraft atmospheric entries, where the heat fluxes encountered can be comparable to those achieved in the DIII-D divertor plasma. Several different designs for the flowfacing side of the rod are analyzed, including "sharp nose," "blunt," and "concave". The last shape is studied for its potential to lower heat fluxes at the rod surface by increased radiation from trapped neutrals and reduced parallel plasma pressure. We also analyze the possibility of applying a thin (approximately 30 microns) layer of silicon carbide (SiC) to the exposed part of several carbon ablation rods to benchmark its erosion calculations and lifetime predictions. Such calculations are of interest as SiC represents a promising material for both thermal protection systems (TPS) and a fusion plasma-facing material (PFM). Preliminary results from the DIII-D rod ablation experiments are also discussed.

Orlov, Dmitri M.↗

High temperature isotropic and anisotropic etching of silicon carbide using forming gas

Plasma-etched micropillars in 4H-SiC were etched in forming gas (4% H2, 96% N2) at 1500, 1550, and 1600 °C at 2.4 and 9.4 standard liters per min (slm). At 2.4 slm, oxygen from the aluminum oxide components of the tube furnace oxidized the SiC surface, and the pillars etched isotropically. Additionally, at 9.4 slm, the pillars etched crystallographically at 1500 and 1550 °C, and sharp 4H-SiC needless with tips as narrow as 15 nm were produced. The lateral etch rates of both 4H- and 6H-SiC were measured with the a-plane etching faster than the m-plane at 1500 and 1550 °C. At 1600 °C, the m-plane and a-plane etch at comparable rates. Due to the difference in etch rates between the m-plane and a-plane, convex surfaces tend to produce facets parallel to the a-plane, while concave surfaces produce facets parallel to the m-plane.

42 ENGINEERING↗

X-ray spectroscopic identification of strain and structure-based resonances in a series of saturated carbon-cage molecules: Adamantane, twistane, octahedrane, and cubane

Novel nanocarbons such as fullerenes, nanotubes, graphene, and nanodiamond reside at the cutting edge of nanoscience and technology. Along with chemical functionalization, geometric constraints (such as extreme curvature in nanotubes or defects within or at the surfaces of diamond nanoparticles) significantly alter the electronic states of the nanocarbon material. Understanding the effects of steric strain on the electronic structure is critical to developing nanoelectronic applications based on these materials. This paper presents a fundamental study of how strain affects the electronic structure in a benchmark series of some fundamental saturated carbon cage compounds. Adamantane, C10H16, the smallest diamondoid and arguably the smallest nanodiamond crystallite, has carbon atoms essentially commensurate with diamond lattice positions and possesses by far the least molecular strain of this series. Twistane also is a C10H16 isomer but the fixed cyclohexane twist conformation of the central ring introduces additional strain into the cage. Octahedrane [(CH)12] and cubane [(CH)8] are considerably more strained, culminating in cubane where carbon–carbon bonds lie either parallel or orthogonal to one another. Using gas-phase near-edge x-ray absorption fine structure spectroscopy to probe the unoccupied electronic states, we observe two major progressions across this series. First, a broad C–C σ* resonance in the absorption splits into two more narrow and intense resonances with increasing strain. Second, the first manifold of states previously associated with tertiary C–H σ* in the diamondoid series appears to broaden and shift to lower energy. This feature is more than twice as intense in cubane than in octahedrane, even though these two molecules have only tertiary carbons, with the chemical formula (CH)x. The spectral differences are entirely due to the shape of the molecules; in particular, in cubane, the features arise from a high degree of p-p interaction between parallel C–C bonds. In contrast to the conventional wisdom that near-edge x-ray absorption is primarily an atomically localized spectroscopy, molecular shape and associated strain lead to the dominant features in spectra acquired from this fundamental series of carbon cage structures.

36 MATERIALS SCIENCE↗

Atomic-resolution lithography with an on-chip scanning tunneling microscope

In this work, atomic-resolution lithography with a Microelectromechanical-System (MEMS) based Scanning Tunneling Microscope (STM) is demonstrated for the first time. The microscope consists of a commercial UltraHigh-Vacuum (UHV) STM whose regular tip is replaced with a 1-Degree-of-Freedom (1-DOF) MEMS nanopositioner. This results in a hybrid STM system where XY-plane motions are provided by the piezotube of the original system and Z-axis motion by the MEMS with a higher bandwidth. Sharp tips made of Pt or W are added to the MEMS devices with postfabrication techniques. With this hybrid system, STM-based lithography is demonstrated on an H-passivated Si (100)-2×1 sample under UHV condition. Results prove the capability of the hybrid STM system for atomic-scale lithography. This capability, paired with the small footprint of the MEMS device, makes this approach a candidate for building a high-throughput parallel STM lithography platform by incorporating an array of 1-DOF MEMS devices that perform lithography in parallel.

Alipour, Afshin (ORCID:0000000263252067)↗

Bias-Modulated ALD of ZnO: Insights into Precursor-Surface Interactions for ZnO Films

Atomic layer deposition (ALD) is widely used to deposit conformal thin films but is often limited in the tunability of the resulting material’s properties. Substrate bias and electric fields alter precursor-surface interactions and provide means to tune material properties. To explore this, we performed zinc oxide (ZnO) ALD using diethylzinc (DEZ) and water on silicon native oxide substrates at 150 °C in a sample holder designed to create a static electrical field by biasing one plate of a parallel plate capacitor-style sample holder during deposition. ZnO films prepared in an electric field/on a biased sample holder were thinner, changed relative crystalline composition, and contained more carbon compared to samples grown in identical sample holders without bias. The thickness was independent of the magnitude of the eletric field between plates, indicating that the primary driver for the change was substrate biasing not the electric field between plates of the parallel plate capacitor-style sample holder. Density functional theory calculations showed enhanced electron migration between dissociatively adsorbed DEZ molecules and the ZnO (002) facet with increasing force from an electric field at the substrate surface, which strengthens the electronic interactions between the surface and the adsorbate. These models offer a compelling explanation for the inhibited growth, changes in the crystallinity, and increase in carbon content of films grown in an electric field/on biased plates.

Jones, Jessica C. (ORCID:0000000174754620)↗

Femtojoule optical nonlinearity for deep learning with incoherent illumination

Optical neural networks (ONNs) are a promising computational alternative for deep learning due to their inherent massive parallelism for linear operations. However, the development of energy-efficient and highly parallel optical nonlinearities, a critical component in ONNs, remains an outstanding challenge. Here, we introduce a nonlinear optical microdevice array (NOMA) compatible with incoherent illumination by integrating the liquid crystal cell with silicon photodiodes at the single-pixel level. We fabricate NOMA with more than half a million pixels, each functioning as an optical analog of the rectified linear unit at ultralow switching energy down to 100 femtojoules per pixel. With NOMA, we demonstrate an optical multilayer neural network. Our work holds promise for large-scale and low-power deep ONNs, computer vision, and real-time optical image processing.

36 MATERIALS SCIENCE↗

Stacking-engineered ferroelectricity in bilayer boron nitride

Stacking a ferroelectric Properties of layered van der Waals structures can depend sensitively on the stacking arrangement of constituent layers. This phenomenon has been exploited to engineer superconducting, correlated insulator, and magnetic states. Two groups now show that ferroelectricity can also be engineered through stacking: Parallel-stacked bilayers of hexagonal boron nitride exhibit ferroelectric switching even though the bulk material is not ferroelectric (see the Perspective by Tsymbal). To explore these phenomena, Yasuda et al. used transport measurements, whereas Vizner Stern et al. focused on atomic force microscopy. Science , abd3230 and abe8177, this issue p. 1458 and p. 1462 ; see also abi7296, p. 1389

Science & Technology - Other Topics↗

Hybridization capture sequencing for Vibrio spp. and associated virulence factors

ABSTRACT Proliferation ofVibriospp. in aquatic ecosystems is associated with climate change and, concomitantly, increased incidence of vibriosis. They are autochthonous to aquatic environments globally, but traditional metagenomic methods for detecting and typing pathogenicVibriospp. are challenged by their presence in relatively low abundance and ability to persist in a viable but nonculturable state. In the study reported here, hybridization capture sequencing (HCS) was employed to profile low-abundanceVibriospp. in environmental samples. The HCS panel targeted a family of molecular chaperones (CPN60) specific to 69Vibriospp. and 162Vibrio-specific virulence factors. This approach was evaluated in parallel with traditional whole-community shotgun sequencing in a metagenomic analysis of water and oyster samples collected from the Chesapeake Bay. In addition,Vibrio parahaemolyticusandVibrio vulnificusstrains isolated from the samples were subjected to whole-genome sequencing to determine the genetic characteristics of pathogenicVibriospp. circulating in an aquatic environment. HCS, employed to determine the incidence and characterization of specificVibriospp., yielded significantly greater metagenomic insight, notably a variety of otherVibriospp., including detection ofVibrio cholerae,Vibrio fluvialis, andVibrio aestuarianus, in addition toVibrio parahaemolyticusandVibrio vulnificus, and also important virulence factors not detectable using traditional molecular methods. Thus, pathogenicVibriospp. in aquatic ecosystems may be far more common than currently understood. It is concluded that environmental surveillance should include HCS, a valuable tool for the detection and characterization of pathogenic agents in aquatic ecosystems, notably vibrios. IMPORTANCE The increasing prevalence of pathogenicVibriospp. in aquatic ecosystems, driven by climate change, is closely linked to a rise in cholera and vibriosis cases, emphasizing the need for improved environmental surveillance. Vibrios are naturally occurring in aquatic environments globally, but traditional metagenomic methods for detecting and typing pathogenicVibriospp. are challenged by their presence in relatively low abundance and ability to persist in a viable but nonculturable state. In the study reported here, hybridization capture sequencing was employed to profile low-abundanceVibriospp. in metagenomic samples, namely water and oysters collected from the Chesapeake Bay. This approach was evaluated in parallel with traditional whole-community shotgun sequencing and whole-genome sequencing ofVibrio parahaemolyticusandVibrio vulnificusstrains isolated from the samples. Results suggest pathogenicVibriospp. in aquatic ecosystems may be far more common than currently understood, when multiple methods are considered for environmental surveillance.

Microbiology↗

Butterfly Factorization Via Randomized Matrix-Vector Multiplications

This paper presents an adaptive randomized algorithm for computing the butterfly factorization of an m × n matrix with m ≈ n provided that both the matrix and its transpose can be rapidly applied to arbitrary vectors. The resulting factorization is composed of O(log n) sparse factors, each containing O(n) nonzero entries. The factorization can be attained using O(n 3/2 log n) computation and O(n log n) memory resources. Furthermore, the proposed algorithm can be implemented in parallel and can apply to matrices with strong or weak admissibility conditions arising from surface integral equation solvers as well as multi-frontal-based finite-difference, finite-element, or finite-volume solvers. A distributed-memory parallel implementation of the algorithm demonstrates excellent scaling behavior.

97 MATHEMATICS AND COMPUTING↗

A New Semistructured Algebraic Multigrid Method

Multigrid methods are well suited to large massively parallel computer architectures because they are mathematically optimal and display good parallelization properties. Since current architecture trends are favoring regular compute patterns to achieve high performance, the ability to express structure has become much more important. The hypre software library provides high-performance multigrid preconditioners and solvers through conceptual interfaces, including a semistructured interface that describes matrices primarily in terms of stencils and logically structured grids. This paper presents a new semistructured algebraic multigrid (SSAMG) method built on this interface. The numerical convergence and performance of a CPU implementation of this method are evaluated for a set of semistructured problems. In conclusion, SSAMG achieves significantly better setup times than hypre’s unstructured AMG solvers and comparable convergence. In addition, the new method is capable of solving more complex problems than hypre’s structured solvers.

97 MATHEMATICS AND COMPUTING↗

Additive Polynomial Time Integrators, Part I: Framework and Fully Implicit-Explicit Collocation Methods

In this paper we generalize the polynomial time integration framework to additively partitioned initial value problems. The framework we present is general and enables the construction of many new families of additive integrators with arbitrary order-of-accuracy and varying degree of implicitness. In this first work, we focus on a new class of implicit-explicit polynomial block methods that are based on fully implicit Runge–Kutta methods with Radau nodes and that possess high stage order. Here, we show that the new fully implicit-explicit (FIMEX) integrators have improved stability compared to existing IMEX Runge–Kutta methods, while also being more computationally efficient due to recent developments in preconditioning techniques for solving the associated systems of nonlinear equations. For PDEs on periodic domains where the implicit component is trivial to invert, we will show how parallelization of the right-hand side evaluations can be exploited to obtain significant speedup compared to existing serial IMEX Runge–Kutta methods. For parallel (in space) finite element discretizations, the new methods can achieve orders of magnitude better accuracy than existing IMEX Runge–Kutta methods and/or achieve a given accuracy several times times faster in terms of computational runtime.

97 MATHEMATICS AND COMPUTING↗

A Performance Portable, Fully Implicit Landau Collision Operator with Batched Linear Solvers

Modern accelerators use hierarchical parallel programming models that enable massive multithreading within a processing element (PE), with multiple PEs per device driven by traditional processes. Batching is a technique for exposing PE-level parallelism in algorithms that have traditionally run on MPI processes or multiple threads within a single process. Opportunities for batching arise in, for example, kinetic discretizations of magnetized plasmas where collisions are advanced in velocity space at each spatial point independently. This paper builds on previous work on a high-performance, fully nonlinear, Landau collision operator by batching the linear solver, as well as batching the spatial point problems and adding new support for multiple grids for multiscale, multispecies problems. An anisotropic relaxation verification test that agrees well with previously published results and analytical models is presented. The performance results from NVIDIA A100 and AMD MI250X nodes are presented with hardware utilization analysis for each architecture. Finally, the entire implicit Landau operator time advance is implemented in Kokkos for performance portability, running entirely on the device and is available in the PETSc numerical library.

97 MATHEMATICS AND COMPUTING↗