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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Toucan: A performance portable, scalable implementation of the DECA algorithm

In the field of additive manufacturing (AM), cellular automata (CA) is extensively used to simulate microstructural evolution during solidification. However, while traditional CA approaches are relatively fast, they still require a substantial number of time steps, are limited to moderate volumes, and are relatively difficult to improve through parallelism due to the highly localized nature of the solidification front. Here, to address these issues of time to solution and load balancing, we introduce Toucan, a parallel, performance-portable, and scalable code written in C++ with the Kokkos library that leverages the discrete event inspired cellular automata (DECA) algorithm to perform parallel-in-time (PinT) grain growth simulations. Toucan effectively mitigates load balancing issues by distributing the computational workload more evenly across processors, enhancing scalability and efficiency. We conduct both strong and weak scaling studies on up to 64 GPUs on the Frontier supercomputer, demonstrating that Toucan significantly outperforms the current state-of-the-art, time-stepped CA code, ExaCA, on both single and multi-GPU simulations. Even in AM-specific weak scaling scenarios, Toucan maintains near-ideal scaling, in contrast to the linear increase observed with ExaCA due to the moving laser raster pattern. This study highlights Toucan’s potential to transform microstructural simulations in AM by radically improving both efficiency and scalability over existing methods.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Graphene Oxide‐Template Controlled Cuboid‐Shaped High‐Capacity VS 4 Nanoparticles as Anode for Sodium‐Ion Batteries

Abstract Room‐temperature sodium‐ion batteries have attracted great attentions for large‐scale energy storage applications in renewable energy. However, exploring suitable anode materials with high reversible capacity and cyclic stability is still a challenge. The VS 4 , with parallel quasi‐1D chains structure of V 4+ (S 2 2− ) 2 , which provides large interchain distance of 5.83 Å and high capacity, has showed great potential for sodium storage. Here, the uniform cuboid‐shaped VS 4 nanoparticles are prepared as anode for sodium‐ion batteries by the controllable of graphene oxide (GO)‐template contents. It exhibits superb electrochemical performances of high‐specific charge capacity (≈580 mAh·g −1 at 0.1 A·g −1 ), long‐cycle‐life (≈98% retain at 0.5 A·g− 1 after 300 cycles), and high rates (up to 20 A·g −1 ). In addition, electrolytes are optimized to understand the sodium storage mechanism. It is thus demonstrated that the findings have great potentials for the applications in high‐performance sodium‐ion batteries.

Wang, Sizhe↗
DOE OSTI2023

Robust multi‐loop control of a floating wind turbine

Abstract A principal challenge facing the control of floating offshore wind turbines (FOWTs) is the problem of instability, or “negative damping,” when using blade pitch feedback to control generator speed. This closed‐loop instability can be attributed to non‐minimum phase zeros in the transfer function from blade pitch to generator speed. Standard approaches to improving stability and performance include robust tuning of control gains and introducing multiple feedback loops to respond to platform motion. Combining these approaches is nontrivial because multiple control loops complicate the impact of coupling in the system dynamics. The single‐loop approach to analyzing stability robustness neglects inter‐loop coupling, while a simplistic multi‐loop approach is highly sensitive to dimensional scaling and overestimates the robustness of the single‐loop controller. This work proposes a sensitivity representation that separates some of the natural FOWT dynamic coupling into a parallel feedback loop in the sensitivity function loop to address both of these concerns. The modified robustness measure is used with a simplified linear FOWT model to optimize scheduled multi‐loop control parameters in an automated tuning procedure. This controller is implemented for the 10‐MW Ultraflexible Smart FLoating Offshore Wind Turbine (USFLOWT) and compared against conventional single‐ and multi‐loop controllers tuned using frequency‐domain analysis and high‐fidelity OpenFAST simulations. The multi‐loop robust controller shows the highest overall performance in generator speed regulation and tower load reduction, though consideration of power quality, actuator usage, and other structural loading leads to additional trade‐offs.

Stockhouse, David↗
DOE OSTI2017

Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration

Turbulent Richtmyer–Meshkov instability (RMI) is investigated through a series of high resolution three-dimensional simulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regimes combined. Overall, the growth rate exponent is determined as θ=0.292±0.009, in good agreement with prior studies; however, the exponent is decaying slowly in time. Also, θ is shown to be relatively insensitive to the choice of mixing layer width measurements. The asymptotic integral molecular mixing measures Θ=0.792±0.014, Ξ=0.800±0.014, and Ψ=0.782±0.013 are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49% and 66% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single θ. In conclusion, quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Extended Kohler's rule of magnetoresistance in TaCo 2 Te 2

TaCo 2 Te 2 is recently reported to be an air-stable, high mobility van der Waals material with probable magnetic order. Here we investigate the scaling behavior of its magnetoresistance. Here we measured both the longitudinal (p xx ) and Hall (p xy ) magnetoresistivities of TaCo 2 Te 2 crystals in magnetic fields parallel to the c axis and found that the magnetoresistance violates the Kohler's rule MR similar to f [H/p 0 ] while obeying the extended Kohler's rule MR ~ f [H/(n T p 0 )], where MR ~ [p xx (H) - p 0 ]/p 0 , H is the magnetic field, n T is a thermal factor, and p xx (H) and p 0 are the resistivities at H and zero field, respectively. Further, while deviating from those of the densities of electrons (n e ) and holes (n h ) obtained from the two-band model analysis of the magnetoconductivities, the temperature dependence of n T is close to that of the Hall carrier densities n H calculated from the slopes of p xy (H) curves at low magnetic fields, providing a different way to obtain the thermal factor in the extended Kohler's rule.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HPDR: High-Performance Portable Scientific Data Reduction Framework

The rapid growth in scientific data generation is outpacing advancements in computing systems necessary for efficient storage, transfer, and analysis, particularly in the context of exascale computing. With the deployment of first-generation exascale computing systems and next-generation experimental facilities, this gap is widening and necessitates effective data reduction techniques to manage enormous data volumes. Over the past decade, various data reduction methods, including lossless compression, error-controlled lossy compression, and data refactoring, have been developed to accelerate I/O in scientific workflows. Despite significant reductions in data volume, these methods introduce considerable computational overhead, which can become the new bottleneck in data processing. To mitigate this, GPU-accelerated data reduction algorithms have been introduced. However, challenges remain in their integration into exascale workflows, including limited portability across different GPU architectures, substantial memory transfer overhead, and reduced scalability on dense multi-GPU systems. To address these challenges, we propose HPDR, a high-performance and portable data reduction framework. HPDR is designed to enable the execution of state-of-the-art reduction algorithms across diverse processor architectures while reducing memory transfer overhead to 2.3 % of the original, resulting in up to 3.5× faster throughput compared to existing solutions. It also achieves up to 96% of the theoretical speedup in multi-GPU settings. In addition, evaluations on accelerating I/O operations at scale up to 1,024 nodes of the Frontier supercomputer demonstrate that HPDR can achieve up to 103 TB/s reduction throughput, providing up to 4× acceleration in parallel I/O performance compared to existing data reduction routines. This work highlights the potential of HPDR to significantly enhance data reduction efficiency in exascale computing environments.

Chen, Jieyang [University of Oregon]↗
DOE OSTI2022

Attribution of North American Subseasonal Precipitation Prediction Skill

The skill of NOAA’s official monthly U.S. precipitation forecasts (issued in the middle of the prior month) has historically been low, having shown modest skill over the southern United States, but little or no skill over large portions of the central United States. The goal of this study is to explain the seasonal and regional variations of the North American subseasonal (weeks 3–6) precipitation skill, specifically the reasons for its successes and its limitations. The performances of multiple recent-generation model reforecasts over 1999–2015 in predicting precipitation are compared to uninitialized simulation skill using the atmospheric component of the forecast systems. This parallel analysis permits attribution of precipitation skill to two distinct sources: one due to slowly evolving ocean surface boundary states and the other to faster time-scale initial atmospheric weather states. A strong regionality and seasonality in precipitation forecast performance is shown to be analogous to skill patterns dictated by boundary forcing constraints alone. The correspondence is found to be especially high for the North American pattern of the maximum monthly skill that is achieved in the reforecast. The boundary forcing of most importance originates from tropical Pacific SST influences, especially those related to El Niño–Southern Oscillation. Furthermore, we discuss physical constraints that may limit monthly precipitation skill and interpret the performance of existing models in the context of plausible upper limits.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

AthenaK: A Performance-portable Version of the Athena++ Adaptive Mesh Refinement Framework

We describe AthenaK: a new implementation of the Athena++ block-based adaptive mesh refinement framework using the Kokkos programming model. Finite volume methods for Newtonian, special relativistic, and general relativistic (GR) hydrodynamics and magnetohydrodynamics (MHD), and GR-radiation hydrodynamics and MHD, as well as a module for evolving Lagrangian tracer or charged test particles (e.g., cosmic rays) are implemented using the framework. In two companion papers, we describe (1) a new solver for the Einstein equations based on the Z4c formalism, and (2) a GRMHD solver in dynamical spacetimes also implemented using the framework, enabling new applications in numerical relativity. By adopting Kokkos, the code can be run on virtually any hardware, including CPUs, GPUs from multiple vendors, and emerging Advanced RISC Machine processors. AthenaK shows excellent performance and weak scaling, achieving over 1 billion cell updates per second for hydrodynamics in three dimensions on a single NVIDIA Grace Hopper processor. It does this with a typical parallel efficiency of 80% on 65,536 AMD GPUs on the OLCF Frontier system. Such performance portability enables AthenaK to leverage modern exascale computing systems for challenging applications in astrophysical fluid dynamics, numerical relativity, and multimessenger astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1974

Assessment of ERTS-1 imagery as a tool for regional geological analysis in New York State

The author has identified the following significant results. Linear anomalies dominate the new geological information derived from ERTS-1 imagery, total lengths now exceeding 26,500 km. Maxima on rose diagrams for ERTS-1 anomalies correspond well with those for mapped faults and topographic lineaments. Multi-scale analysis of linears shows that single topographic linears at 1:2,500,000 may become dashed linears at 1:1,000,000 aligned zones of shorter parallel, en echelon, or conjugate linears at 1:5000,000, and shorter linears lacking any conspicuous zonal alignment at 1:250,000. Field work in the Catskills suggests that the prominent new NNE lineaments may be surface manifestations of dip slip faulting in the basement, and that it may become possible to map major joint sets over extensive plateau regions directly on the imagery. Most circular features found were explained away by U-2 airfoto analysis but several remain as anomalies. Visible glacial features include individual drumlins, drumlinoids, eskers, ice-marginal drainage channels, glacial lake shorelines, sand plains, and end moraines.

Isachsen, Y. W.↗
NASA NTRS1982

Opportunities for composites in commercial transport structures

Manufacturers are developing composite versions of structural components on existing aircraft. Development involves testing of various material options before selecting one and then extensive testing to develop an adequate data base of material strength and stiffness properties. Design options are narrowed through analysis and a varied spectrum of development tests on small and large subcomponents. In parallel with this, a suitable production process including economical ply preparation and cure at high temperature and pressure is evolved, tools are designed and fabricated, and full scale components are then manufactured for ground qualification tests, flight tests, and airline service. The various tests include many that are required by the FAA for flight certification, which must precede airline service. Inspection and repair methods to insure adequate maintenance in service are also developed.

Herman L. Bohon↗
NASA NTRS1982

Major disruptions, inverse cascades, and the Strauss equations

Current-carrying plasmas in a strong d.c. magnetic field are subject to violent disruptions above certain thresholds. At present difficult to verify, explanations are typically sought in terms of 'tearing modes'. An alternative explanation is in terms of inverse magnetic helicity cascades, generated from a variety of possible sources of small-scale MHD turbulence. Strongly anisotropic MHD plasmas may be described by the Strauss equations. Indications of turbulent inverse cascade behavior for the Strauss equations are sought, in parallel with earlier examples from MHD and fluid mechanics.

Montgomery, D.↗
NASA NTRS1983

Numerical modeling of turbulent flow in a channel

Two-dimensional incompressible turbulent flow in a channel with a backward-facing step was studied numerically by Chorin's Random Vortex Method (RVM), an algorithm capable of tracing the action of elementary turbulent eddies and their cumulative effects without imposing any restrictions upon their motions. The step occurs in one side of a channel with otherwise flat, parallel walls; its height equals 1/3, 1/4 or 1/5 the width of the channel downstream. The main objective was to investigate the behavior of the large-scale turbulent eddies in a flow and the flow characteristics in the separated shear layer, the reattached zone, and the rebuilding boundary layer after reattachment. The unsteady vorticity field and the distribution of time-averaged turbulent statistics were obtained. The effects of expansion step height and initial boundary layer state were also studied. Comparisons were made with the available experimental results. The agreement is satisfactory in the velocity profiles and in the reattachment length, and fairly good in the turbulence profiles. Also a mechanism of the development of the reattaching turbulent flow was suggested by the numerical results.

Dai, Y. W.↗
NASA NTRS1986

Optimal design of structures with multiple design variables per group and multiple loading conditions on the personal computer

A finite element based programming system for minimum weight design of a truss-type structure subjected to displacement, stress, and lower and upper bounds on design variables is presented. The programming system consists of a number of independent processors, each performing a specific task. These processors, however, are interfaced through a well-organized data base, thus making the tasks of modifying, updating, or expanding the programming system much easier in a friendly environment provided by many inexpensive personal computers. The proposed software can be viewed as an important step in achieving a 'dummy' finite element for optimization. The programming system has been implemented on both large and small computers (such as VAX, CYBER, IBM-PC, and APPLE) although the focus is on the latter. Examples are presented to demonstrate the capabilities of the code. The present programming system can be used stand-alone or as part of the multilevel decomposition procedure to obtain optimum design for very large scale structural systems. Furthermore, other related research areas such as developing optimization algorithms (or in the larger level: a structural synthesis program) for future trends in using parallel computers may also benefit from this study.

Nguyen, D. T.↗
NASA NTRS1987

Anomalous transport in discrete arcs and simulation of double layers in a model auroral circuit

The evolution and long-time stability of a double layer in a discrete auroral arc requires that the parallel current in the arc, which may be considered uniform at the source, be diverted within the arc to charge the flanks of the U-shaped double-layer potential structure. A simple model is presented in which this current re-distribution is effected by anomalous transport based on electrostatic lower hybrid waves driven by the flank structure itself. This process provides the limiting constraint on the double-layer potential. The flank charging may be represented as that of a nonlinear transmission. A simplified model circuit, in which the transmission line is represented by a nonlinear impedance in parallel with a variable resistor, is incorporated in a 1-d simulation model to give the current density at the DL boundaries. Results are presented for the scaling of the DL potential as a function of the width of the arc and the saturation efficiency of the lower hybrid instability mechanism.

Smith, Robert A.↗
NASA NTRS1988

Vectorization of a particle simulation method for hypersonic rarefied flow

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.

Mcdonald, Jeffrey D.↗
NASA NTRS1992

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗
NASA NTRS1992

Characteristics of obliquely propagating, nonlinear Alfven waves

Consideration is given to properties of obliquely propagating (with respect to the large-scale magnetic field) nonlinear Alfven and fast magnetosonic waves as described by the derivative nonlinear Schroedinger equation. The focus is on those waves which in the parallel-propagating limit are circularly polarized, exact solutions to the MHD equations. Results are presented of a study of the changes in the structure of these waves as the propagation angle is increased. Attention is also given to the evolution of initial conditions which do not correspond to stationary wave propagations. It is found that, as the propagation angle with respect to the mean field increases, the functional forms of the wave magnetic field components in the two directions orthogonal to the propagation direction become increasingly dissimilar. The angle about the mean field within which stationary waves resemble parallel-propagating, circularly polarized waves becomes smaller as the wavelength increases. If the 'initial condition' wave train is close to parallel-propagating, then there is little subsequent evolution in its form.

Spangler, Steven R.↗
NASA NTRS2003

The FUSE Survey of 0 VI in the Galactic Halo

This paper summarizes the results of the Far-Ultraviolet Spectroscopic Explorer (FUSE) program to study 0 VI in the Milky Way halo. Spectra of 100 extragalactic objects and two distant halo stars are analyzed to obtain measures of O VI absorption along paths through the Milky Way thick disk/halo. Strong O VI absorption over the velocity range from -100 to 100 km/s reveals a widespread but highly irregular distribution of O VI, implying the existence of substantial amounts of hot gas with T approx. 3 x 10(exp 5) K in the Milky Way thick disk/halo. The overall distribution of O VI is not well described by a symmetrical plane-parallel layer of patchy O VI absorption. The simplest departure from such a model that provides a reasonable fit to the observations is a plane-parallel patchy absorbing layer with an average O VI mid-plane density of n(sub 0)(O VI) = 1.7 x 10(exp -2)/cu cm, a scale height of approx. 2.3 kpc, and a approx. 0.25 dex excess of O VI in the northern Galactic polar region. The distribution of O VI over the sky is poorly correlated with other tracers of gas in the halo, including low and intermediate velocity H I, Ha emission from the warm ionized gas at approx. l0(exp 4) K, and hot X-ray emitting gas at approx. l0(exp 6) K . The O VI has an average velocity dispersion, b approx. 60 km/s and standard deviation of 15 km/s. Thermal broadening alone cannot explain the large observed profile widths. A combination of models involving the radiative cooling of hot fountain gas, the cooling of supernova bubbles in the halo, and the turbulent mixing of warm and hot halo gases is required to explain the presence of O VI and other highly ionized atoms found in the halo. The preferential venting of hot gas from local bubbles and superbubbles into the northern Galactic polar region may explain the enhancement of O VI in the North.

Sonneborn, George↗