Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel, scaling”

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,693 records · Page 94

DIMPLES: Distributed Influence Maximization for Pandemic pLanning on Exascale Systems

We study exascale parallel algorithms for the selection of intervention or monitoring strategies in massive realistic socio-technical networks through scalable Influence Maximization (InfMax) algorithms. We employ novel techniques to enable efficient scaling on up to 8k nodes of OLCF Frontier, with 65k AMD GPUs and 458k AMD CPU cores. Current state-of-the-art InfMax tools are limited to networks with only a few million actors (vertices) and a few hundred million interactions (edges). By overcoming these limitations, we show that our approach is capable of processing a realistic social contact network of the United States with 285 million nodes and about 8 billion edges. This two orders-of-magnitude improvement over the previous state-of-the-art is obtained by leveraging algorithmic advancements for the InfMax problem and designing several problem-specific approaches to overlap communication with computation, improve GPU efficiency, and lower the application’s memory requirements. We evaluate strong scaling for computing 10k most influential seeds using up to 8k nodes of an exascale system, and weak scaling from 128 to 8k system nodes for seed sets ranging from 625 to 40k seeds. We achieve the fastest-known runtime of 25 minutes while performing 48 million diffusion simulations totaling 2.31 petabytes to identify 40k influential seeds using 8k nodes, and take 5.75 minutes to identify 10k seeds while using 4k nodes.

Minutoli, Marco [Pacific Northwest National Labora↗

Alstom’s Limestone Chemical Looping Gasification Process for High Hydrogen Syngas Generation

Alstom Power Inc., a wholly owned subsidiary of the General Electric Company (GE), is developing Limestone-based Chemical Looping (LCL™) Processes, utilizing limestone derived calcium sulfate as an oxygen carrier to transport oxygen from air to the fuel. The LCL™ process is a breakthrough technology with the potential to capture CO 2 and SO 2 from new and existing coal-fired power plants at reduced costs while maintaining high plant power generation efficiency. Chemical looping is a next generation advancement to oxy-combustion which would not require the significant investment in a cryogenic distillation type oxygen production plant. Instead of the energy-consuming oxygen plant, chemical looping uses a benign material to directly remove oxygen from the air and bring that oxygen into contact with the fuel. This "game changing" technology is the lowest cost-of-electricity technical approach that GE Power has identified to date for coal power with carbon capture. The technology has significantly lower capital costs compared to a conventional steam power plant with first generation carbon capture technologies. This new power plant concept is based on a hybrid combustion-gasification process utilizing high temperature chemical and thermal looping technology. In the air reactor, a solid oxygen carrier picks up oxygen from air through an oxidation reaction and leaves nitrogen behind. The hot oxygen-carrying solid carrier is then transported to a fuel reactor, where it releases oxygen. The released oxygen then converts the coal in the fuel reactor into combustion gases. The solid carrier is recycled back to the air reactor for regeneration completing the “chemical loop.” The regeneration step produces heat from which steam is produced for power generation while the process results in a concentrated carbon dioxide stream in which the carbon is captured for use or sequestration (CCUS). Limestone as a precursor for the oxygen carrier is abundant, benign, low cost, and is used in CFB and PC boilers for sulfur control. Challenges associated with this carrier include a complex reaction mechanism with possible side reactions leading to the release of SO 2 . The development of the Limestone Chemical Looping Gasification (LCL-G™) technology was co-funded by the Department of Energy under the DOE/NETL Cooperative Agreement No. DE- FE0023497, in parallel with the collaborative combustion (LCL-C™) program. In this report, the LCL-G™ experimental testing from the electrically-heated, once-through, 100 kW Pilot-Scale Test Facility (PSTF) is analysed and reviewed. The PSTF Oxidizer is a 50-foot tall, 100-mm diameter transport reactor, while the PSTF Reducer reactor is based upon a bubbling-bed design. The experimental data from four (4) test campaigns (denoted “Weeks”) are reviewed and discussed. The LCL-G™ (DE-FE0023497) project objectives for gasification are >90% carbon gasified in the Reducer and a H 2 :CO molar generation ratio of 2:1. These project performance metrics were not met by a single coal test, but they were individually approached in two separate tests.

01 COAL, LIGNITE, AND PEAT↗

Distributed Training for High Resolution Images: A Domain and Spatial Decomposition Approach

In this work we developed two Pytorch libraries using the PyTorch RPC interface for distributed deep learning approaches on high resolution images. The spatial decomposition library allows for distributedtraining on very large images, which otherwise won’t be possible on a single GPU. The domain parallelism library allows for distributed training across multiple domain unlabeled data, by leveraging the domain separation architecture. Both of those libraries were tested on the Summit supercomputer at a moderate scale, and we are releasing the code for both of them.

97 MATHEMATICS AND COMPUTING↗

BM3DORNL

BM3DORNL is a high-performance, open-source library for removing streak and ring artifacts from computed-tomography (CT) data, developed for neutron imaging at Oak Ridge National Laboratory's Spallation Neutron Source (VENUS beamline) and applicable to X-ray CT as well. Ring artifacts — concentric rings in reconstructed slices caused by detector pixel-to-pixel response non-uniformities — appear as vertical streaks in the sinogram and degrade both image quality and quantitative analysis. BM3DORNL operates in the sinogram domain using an adaptation of the BM3D (block-matching and 3D collaborative filtering) algorithm (Dabov et al., 2007). It provides a dedicated streak-removal mode, a true multi-scale BM3D variant (after Mäkinen et al., 2021) that suppresses wide streaks single-scale methods miss, and an alternative Fourier–SVD method (~2.6× faster) combining FFT-based energy detection with rank-1 SVD. The computationally intensive core is implemented in Rust with parallel (Rayon) block matching, integral-image pre-screening, and optimized transforms, and is exposed through a simple Python API (with an optional GUI) so it integrates directly into existing tomography reconstruction pipelines. It processes both 2D sinograms and 3D sinogram stacks, is pip-installable for Linux and macOS, and is documented at https://bm3dornl.readthedocs.io.

Zhang, Chen [Oak Ridge National Laboratory (ORNL),↗

Application of satellite time transfer in autonomous spacecraft clocks

The conceptual design of a spacecraft clock that will provide a standard time scale for experimenters in future spacecraft., and can be sychronized to a time scale without the need for additional calibration and validation is described. The time distribution to the users is handled through onboard computers, without human intervention for extended periods. A group parallel binary code, under consideration for onboard use, is discussed. Each group in the code can easily be truncated. The autonomously operated clock not only achieves simpler procedures and shorter lead times for data processing, but also contributes to spacecraft autonomy for onboard navigation and data packetization. The clock can be used to control the sensor in a spacecraft, compare another time signal such as that from the global positioning system, and, if the cost is not a consideration, can be used on the ground in remote sites for timekeeping and control.

Chi, A. R.↗

What happens to the large eddies when net drag reduction is achieved by outer flow manipulators?

The tandem-arranged parallel plate manipulator (TAPPM) presently used to generate a flow with net drag reduction is conditionally sampled and found to exhibit significant changes of the boundary layer's large-scale motions at 20 delta(0), where delta(0) is the boundary layer thickness at the upstream plate of the TAPPM. Flow fields in the 'valleys' are equally affected. Both turbulent and irrotational flows then reestablish themselves by 50 delta(0). Flow visualization indicates that the TAPPM wake is very coherent at 20 delta(0), and has not yet spread into the wall region, while large scale motions and the ambient response flow continue to exhibit dynamic changes. This indicates that the shielding effect of the TAPPM's wake prolongs suppression of large-scale motions, thereby reducing skin friction over a sufficiently long extent of the boundary layer to overcome device drag and yield net drag reductions.

Falco, R. E.↗

Dense gas and HII regions in the starburst galaxy NGC 253

The energetic activity in the nuclear barred region of NGC 253 is attributable to a burst of star formation. NGC 253 is in many ways a twin of the prototypical starburst galaxy M82; the strong non-thermal radio continuum, high far-infrared luminosity, and bright molecular emission of the central 1 Kpc parallel the morphology of the M82 starburst. Furthermore, the filamentary low ionization optical emission and extended x ray emission along the minor axis in NGC 253 is similar to a scaled down version of the well developed galactic bipolar wind in M82. The infrared luminosity of NGC 253, 3(exp 10) solar luminosity, is comparable to M82 but is emitted from a smaller region (Telesco and Harper 1980). This suggests that the NGC 253 starburst may be more intense and at an earlier evolutionary stage than M82. However, the presence of a non-stellar AGN in NGC 253 may complicate the comparison (Turner and Ho, 1985). Researchers used the Hat Creek millimeter interferometer to map emission from the J = 1 to 0 transitions of HCN and HCO(+) as well as 3 mm continuum emission, toward the nuclear region of NGC 253. The HCO(+) and continuum observations are sensitive to spatial scales from 6 to 45 seconds. The 2 minute field of view comfortably includes the entire starburst region (about 40 seconds; 650 pc). Because the longer baseline HCN observations are not yet complete, they are only sensitive to spatial scales from 15 to 45 seconds.

Carlstrom, J. E.↗

Application of computational aero-acoustics to real world problems

The application of computational aeroacoustics (CAA) to real problems is discussed in relation to the analysis performed with the aim of assessing the application of the various techniques. It is considered that the applications are limited by the inability of the computational resources to resolve the large range of scales involved in high Reynolds number flows. Possible simplifications are discussed. It is considered that problems remain to be solved in relation to the efficient use of the power of parallel computers and in the development of turbulent modeling schemes. The goal of CAA is stated as being the implementation of acoustic design studies on a computer terminal with reasonable run times.

Hardin, Jay C.↗

Emissions Prediction and Measurement for Liquid-Fueled TVC Combustor with and without Water Injection

An investigation is performed to evaluate the performance of a computational fluid dynamics (CFD) tool for the prediction of the reacting flow in a liquid-fueled combustor that uses water injection for control of pollutant emissions. The experiment consists of a multisector, liquid-fueled combustor rig operated at different inlet pressures and temperatures, and over a range of fuel/air and water/fuel ratios. Fuel can be injected directly into the main combustion airstream and into the cavities. Test rig performance is characterized by combustor exit quantities such as temperature and emissions measurements using rakes and overall pressure drop from upstream plenum to combustor exit. Visualization of the flame is performed using gray scale and color still photographs and high-frame-rate videos. CFD simulations are performed utilizing a methodology that includes computer-aided design (CAD) solid modeling of the geometry, parallel processing over networked computers, and graphical and quantitative post-processing. Physical models include liquid fuel droplet dynamics and evaporation, with combustion modeled using a hybrid finite-rate chemistry model developed for Jet-A fuel. CFD and experimental results are compared for cases with cavity-only fueling, while numerical studies of cavity and main fueling was also performed. Predicted and measured trends in combustor exit temperature, CO and NOx are in general agreement at the different water/fuel loading rates, although quantitative differences exist between the predictions and measurements.

Brankovic, A.↗

Using the GeoFEST Faulted Region Simulation System

GeoFEST (the Geophysical Finite Element Simulation Tool) simulates stress evolution, fault slip and plastic/elastic processes in realistic materials, and so is suitable for earthquake cycle studies in regions such as Southern California. Many new capabilities and means of access for GeoFEST are now supported. New abilities include MPI-based cluster parallel computing using automatic PYRAMID/Parmetis-based mesh partitioning, automatic mesh generation for layered media with rectangular faults, and results visualization that is integrated with remote sensing data. The parallel GeoFEST application has been successfully run on over a half-dozen computers, including Intel Xeon clusters, Itanium II and Altix machines, and the Apple G5 cluster. It is not separately optimized for different machines, but relies on good domain partitioning for load-balance and low communication, and careful writing of the parallel diagonally preconditioned conjugate gradient solver to keep communication overhead low. Demonstrated thousand-step solutions for over a million finite elements on 64 processors require under three hours, and scaling tests show high efficiency when using more than (order of) 4000 elements per processor. The source code and documentation for GeoFEST is available at no cost from Open Channel Foundation. In addition GeoFEST may be used through a browser-based portal environment available to approved users. That environment includes semi-automated geometry creation and mesh generation tools, GeoFEST, and RIVA-based visualization tools that include the ability to generate a flyover animation showing deformations and topography. Work is in progress to support simulation of a region with several faults using 16 million elements, using a strain energy metric to adapt the mesh to faithfully represent the solution in a region of widely varying strain.

Geophyical Finite Element Simulation Tool (GeoFEST↗

Cold Ion Demagnetization near the X-line of Magnetic Reconnection

Although the effects of magnetic reconnection in magnetospheres can be observed at planetary scales, reconnection is initiated at electron scales in a plasma. Surrounding the electron diffusion region, there is an Ion-Decoupling Region (IDR) of the size of the ion length scales (inertial length and gyroradius). Reconnection at the Earths magnetopause often includes cold magnetospheric (few tens of eV), hot magnetospheric (10 keV), and magnetosheath (1 keV) ions, with different gyroradius length scales. We report observations of a subregion inside the IDR of the size of the cold ion population gyroradius (approx. 15 km) where the cold ions are demagnetized and accelerated parallel to the Hall electric field. Outside the subregion, cold ions follow the E x B motion together with electrons, while hot ions are demagnetized. We observe a sharp cold ion density gradient separating the two regions, which we identify as the cold and hot IDRs.

Toledo-Redondo, Serio↗

Simulation-Based Parameter Optimization Framework for Large-Scale Hybrid Smart Grid Communications Systems Design

The design of reliable, dynamic, fault-tolerant hybrid smart grid communication networks is a challenge to achieve for autonomous power grids. Hybrid networks use different communications technologies for different area networks. A simulation-based parameter optimization framework is proposed to tune parameters of hybrid communication technologies to achieve the optimal network performance. It consists of three main components: a parallel executor used to speedup a list of simulations; a sampler running simulations using the parallel executor at each generation; and a hybrid stochastic optimization algorithm for tuning configurable parameters of hybrid designs and applications. The proposed hybrid metaheuristic optimization algorithm combines an evolutionary algorithm with a gradient method to quickly achieve an approximately global optimum solution. Three optimization test functions are employed to train the adjustable parameters of the hybrid algorithm. Results show the proposed parameter optimization framework can help the designer choose the right hybrid architecture with an optimal parameter set for a large-scale broadband PLC-WiMAX hybrid smart grid communication network.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Instabilities and subharmonic resonances of subsonic heated round jets, volume 2

When a jet is perturbed by a periodic excitation of suitable frequency, a large-scale coherent structure develops and grows in amplitude as it propagates downstream. The structure eventually rolls up into vortices at some downstream location. The wavy flow associated with the roll-up of a coherent structure is approximated by a parallel mean flow and a small, spatially periodic, axisymmetric wave whose phase velocity and mode shape are given by classical (primary) stability theory. The periodic wave acts as a parametric excitation in the differential equations governing the secondary instability of a subharmonic disturbance. The (resonant) conditions for which the periodic flow can strongly destabilize a subharmonic disturbance are derived. When the resonant conditions are met, the periodic wave plays a catalytic role to enhance the growth rate of the subharmonic. The stability characteristics of the subharmonic disturbance, as a function of jet Mach number, jet heating, mode number and the amplitude of the periodic wave, are studied via a secondary instability analysis using two independent but complementary methods: (1) method of multiple scales, and (2) normal mode analysis. It is found that the growth rates of the subharmonic waves with azimuthal numbers beta = 0 and beta = 1 are enhanced strongly, but comparably, when the amplitude of the periodic wave is increased. Furthermore, compressibility at subsonic Mach numbers has a moderate stabilizing influence on the subharmonic instability modes. Heating suppresses moderately the subharmonic growth rate of an axisymmetric mode, and it reduces more significantly the corresponding growth rate for the first spinning mode. Calculations also indicate that while the presence of a finite-amplitude periodic wave enhances the growth rates of subharmonic instability modes, it minimally distorts the mode shapes of the subharmonic waves.

Ng, Lian Lai↗

Influence of Turbulence Parameters, Reynolds Number, and Body Shape on Stagnation-Region Heat Transfer

The purpose of the present work was threefold: (1) to determine if a free-stream turbulence length scale existed that would cause the greatest augmentation in stagnation-region heat transfer over laminar levels; (2) to investigate the effect of velocity gradient on stagnation-region heat transfer augmentation by free-stream turbulence; and (3) to develop a prediction tool for stagnation heat transfer in the presence of free-stream turbulence. Heat transfer was measured in the stagnation region of four models with elliptical leading edges that had ratios of major to minor axes of 1:1, 1.5:1, 2.25:1, and 3:1. Five turbulence-generating grids were fabricated; four were square mesh, biplane grids made from square bars. The fifth grid was an array of fine parallel wires that were perpendicular to the model spanwise direction. Heat transfer data were taken at Reynolds numbers ranging from 37 000 to 228 000. Turbulence intensities were in the range of 1.1 to 15.9% while the ratio of integral length scale to leading-edge diameter ranged from 0.05 to 0.30. Stagnation-point velocity gradient was varied by nearly 50%. Stagnation-region heat transfer augmentation was found to increase with decreasing length scale but no optimum length scale was found. Heat transfer augmentation due to turbulence was found to be unaffected by the velocity gradient near the leading edge. A correlation was developed that fit heat transfer data for the square-bar grids to within +/- 4%.

Vanfossen, G. James↗

Scaling deep learning workloads: NVIDIA DGX-1/Pascal and Intel Knights Landing

Deep Learning (DL) algorithms have become ubiquitous in data analytics. As a result, major computing vendors --- including NVIDIA, Intel, AMD, and IBM --- have architectural road-maps influenced by DL workloads. Furthermore, several vendors have recently advertised new computing products as accelerating large DL workloads. Unfortunately, it is difficult for data scientists to quantify the potential of these different products. This paper provides a performance and power analysis of important DL workloads on two major parallel architectures: NVIDIA DGX-1 (eight Pascal P100 GPUs interconnected with NVLink) and Intel Knights Landing (KNL) CPUs interconnected with Intel Omni-Path or Cray Aries. Our evaluation consists of a cross section of convolutional neural net workloads: CifarNet, AlexNet, GoogLeNet, and ResNet50 topologies using the Cifar10 and ImageNet datasets. The workloads are vendor-optimized for each architecture. Our analysis indicates that although GPUs provide the highest overall performance, the gap can close for some convolutional networks; and the KNL can be competitive in performance/watt. We find that NVLink facilitates scaling efficiency on GPUs. However, its importance is heavily dependent on neural network architecture. Furthermore, for weak-scaling --- sometimes encouraged by restricted GPU memory --- NVLink is less important.

Gawande, Nitin A.↗

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗

Development of Slat-gap and Slat-cove Filler Treatments for Noise Reduction Assessment on the High Lift Common Research Model in the NASA LaRC 14x22

Airframe noise is an important issue during the approach and landing phases of flight because engines are throttled back, and airframe sources become an equal or greater contributor to environmental noise. The leading-edge slat component of the high-lift system for transport aircraft has been identified as a prominent source of airframe noise. Noise production mechanisms and potential treatments to mitigate airframe noise have been areas of active research for nearly three decades. Two promising technologies for slat noise reduction are the slat-gap filler (SGF) and slat-cove filler (SCF). Previous work was done using a 2D section of a transport class aircraft wing to demonstrate their feasibility via benchtop and computational models. A wind tunnel test campaign was recently undertaken to test the SGF and SCF technologies in a flight-like environment by installing dynamically scaled treatments on the 10%-scale High-Lift Common Research Model (CRM-HL), which is representative of a commercial twin-aisle transport aircraft. Advancement of the SGF and SCF technologies to overcome challenges of integration with 3D airframe structures, computational simulations of the CRM-HL predicting the aerodynamic and aeroacoustic performance and structural loads, and design of wind tunnel model hardware with features accommodating noise treatments were conducted in parallel. Details of innovations in the SGF and SCF technologies for 3D airframe feasibility will be described. The loads the treatments must sustain in flight and the influence of those loads on the design is also discussed. Finally, the design of wing leading-edge assemblies to treat the 10% CRM-HL with dynamically scaled SGF and SCF components is presented. The treatments were successfully tested in 2020/2021 providing a demonstration of leading-edge noise reduction using components that would allow slat articulation.

airframe noise reduction↗

Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes

The Advisory Group for Aerospace Research and Development (AGARD) has requested my participation in the lecture series entitled Parallel Computing in Computational Fluid Dynamics to be held at the von Karman Institute in Brussels, Belgium on May 15-19, 1995. In addition, a request has been made from the US Coordinator for AGARD at the Pentagon for NASA Ames to hold a repetition of the lecture series on October 16-20, 1995. I have been asked to be a local coordinator for the Ames event. All AGARD lecture series events have attendance limited to NATO allied countries. A brief of the lecture series is provided in the attached enclosure. Specifically, I have been asked to give two lectures of approximately 75 minutes each on the subject of parallel solution techniques for the fluid flow equations on unstructured meshes. The title of my lectures is "Parallel CFD Algorithms for Aerodynamical Flow Solvers on Unstructured Meshes" (Parts I-II). The contents of these lectures will be largely review in nature and will draw upon previously published work in this area. Topics of my lectures will include: (1) Mesh partitioning algorithms. Recursive techniques based on coordinate bisection, Cuthill-McKee level structures, and spectral bisection. (2) Newton's method for large scale CFD problems. Size and complexity estimates for Newton's method, modifications for insuring global convergence. (3) Techniques for constructing the Jacobian matrix. Analytic and numerical techniques for Jacobian matrix-vector products, constructing the transposed matrix, extensions to optimization and homotopy theories. (4) Iterative solution algorithms. Practical experience with GIVIRES and BICG-STAB matrix solvers. (5) Parallel matrix preconditioning. Incomplete Lower-Upper (ILU) factorization, domain-decomposed ILU, approximate Schur complement strategies.

Barth, Timothy J.↗