Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed 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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Advanced Precipitation Radar Antenna to Measure Rainfall From Space

To support NASA s planned 20-year mission to provide sustained global precipitation measurement (EOS-9 Global Precipitation Measurement (GPM)), a deployable antenna has been explored with an inflatable thin-membrane structure. This design uses a 5.3 5.3-m inflatable parabolic reflector with the electronically scanned, dual-frequency phased array feeds to provide improved rainfall measurements at 2.0-km horizontal resolution over a cross-track scan range of up to 37 , necessary for resolving intense, isolated storm cells and for reducing the beam-filling and spatial sampling errors. The two matched radar beams at the two frequencies (Ku and Ka bands) will allow unambiguous retrieval of the parameters in raindrop size distribution. The antenna is inflatable, using rigidizable booms, deployable chain-link supports with prescribed curvatures, a smooth, thin-membrane reflecting surface, and an offset feed technique to achieve the precision surface tolerance (0.2 mm RMS) for meeting the low-sidelobe requirement. The cylindrical parabolic offset-feed reflector augmented with two linear phased array feeds achieves dual-frequency shared-aperture with wide-angle beam scanning and very low sidelobe level of -30 dB. Very long Ku and Ka band microstrip feed arrays incorporating a combination of parallel and series power divider lines with cosine-over-pedestal distribution also augment the sidelobe level and beam scan. This design reduces antenna mass and launch vehicle stowage volume. The Ku and Ka band feed arrays are needed to achieve the required cross-track beam scanning. To demonstrate the inflatable cylindrical reflector with two linear polarizations (V and H), and two beam directions (0deg and 30deg), each frequency band has four individual microstrip array designs. The Ku-band array has a total of 166x2 elements and the Ka-band has 166x4 elements with both bands having element spacing about 0.65 lambda(sub 0). The cylindrical reflector with offset linear array feeds reduces the complexity from "NxN" transmit/receive (T/R) modules of a conventional planar-phased array to just "N" T/R modules. The antenna uses T/R modules with electronic phase-shifters for beam steering. The offset reflector does not provide poor cross-polarization like a double- curved offset reflector would, and it allows the wide scan angle in one plane required by the mission. Also, the cylindrical reflector with two linear array feeds provides dual-frequency performance with a single, shared aperture. The aperture comprises a reflective surface with a focal length of 1.89 m and is made from aluminized Kapton film. The reflective surface is of uniform thickness in the range of a few thousandths of an inch and is attached to the chain-link support structure via an adjustable suspension system. The film aperture rolls up, together with the chain-link structure, for launch and can be deployed in space by the deployment of the chain-link structure.

Rahmat-Samii, Yahya↗
NASA NTRS1986

Spatial distributions of magnetic field fluctuations in the dayside magnetosheath

In a study that tests the hypothesis that magnetosheath magnetic fields are disturbed on plasma streamlines which are connected to the quasi-parallel bow shock, magnetometer observations from the ISEE 2 and IMP 8 spacecraft are used to investigate the dayside spatial distributions of fluctuating magnetosheath fields for different interplanetary field orientations. The results suggest that although other sources such as Kelvin-Helmholtz instabilities and flux transfer processes probably contribute to the fluctuations in the magnetosheath field, the quasi-parallel shock source is an important contributor in the dayside region.

Luhmann, J. G.↗
NASA NTRS1987

Stability of the polar wind

The classical polar wind is an ambipolar outflow of thermal plasma from the terrestrial ionosphere at high latitudes. At altitudes above about 3000 km, the H(+) flow becomes supersonic and collisionless, and the H(+) velocity distribution becomes non-Maxwellian. The non-Maxwellian features include a temperature anisotropy, with the parallel H(+) temperature greater than the perpendicular temperature, and an asymmetry, with an elongated tail in the upward direction. These distortions from a Maxwellian increase as the H(+) gas escapes in the diverging geomagnetic field, and at 10 earth radii, the parallel-to-perpendicular temperature ratio is about 50 and the elongated tail is sufficiently long to move the drift velocity point off the peak of the distribution function. The stability of these highly non-Maxwellian H(+) velocity distributions was studied with regard to the excitation of electrostatic waves, and the plasma was found to be remarkably stable for a wide range of electron temperatures. This indicates that the various macroscopic formulations of the classical polar wind are valid. The stability of a perpendicularly heated polar wind was also studied, assuming bulk perpendicular heating of H(+) in the cusp, followed by the subsequent convection of the heated plasma into the polar cap. Two regions of instability were found.

Barakat, A. R.↗
NASA NTRS1947

The Monoplane as a Lifting Vortex Surface

In Prandtl's airfoil theory the monoplane was replaced by a single lifting vortex line and yielded fairly practical results. However, the theory remained restricted to the straight wing. Yawed wings and those curved in flight direction could not be computed with this first approximation; for these the chordwise lift distribution must be taken into consideration. For the two-dimensional problem the transition from the lifting line to the lifting surface has been explained by Birnbaum. In the present report the transition to the three-dimensional problem is undertaken. The first fundamental problem involves the prediction of flow, profile, and drag for prescribed circulation distribution on the straight rectangular wing, the yawed wing for lateral boundaries parallel to the direction of flight, the swept-back wing, and the rectangular wing in slipping, with the necessary series developments for carrying through the calculations, the practical range of convergence of which does not comprise the wing tips or the break point of the swept-back wing. The second problem concerns the calculation of the circulation distribution with given profile for a slipping rectangular monoplane with flat profile and aspect ratio 6, and a rectangular wing with cambered profile and variable aspect ratio-the latter serving as check of the so-called conversion formulas of the airfoil theory.

Blenk, Hermann↗
DOE OSTI2023

Distributed out-of-memory NMF on CPU/GPU architectures

We propose an efficient distributed out-of-memory implementation of the non-negative matrix factorization (NMF) algorithm for heterogeneous high-performance-computing systems. The proposed implementation is based on prior work on NMFk, which can perform automatic model selection and extract latent variables and patterns from data. In this work, we extend NMFk by adding support for dense and sparse matrix operation on multi-node, multi-GPU systems. The resulting algorithm is optimized for out-of-memory problems where the memory required to factorize a given matrix is greater than the available GPU memory. Memory complexity is reduced by batching/tiling strategies, and sparse and dense matrix operations are significantly accelerated with GPU cores (or tensor cores when available). Input/output latency associated with batch copies between host and device is hidden using CUDA streams to overlap data transfers and compute asynchronously, and latency associated with collective communications (both intra-node and inter-node) is reduced using optimized NVIDIA Collective Communication Library (NCCL) based communicators. Benchmark results show significant improvement, from 32X to 76x speedup, with the new implementation using GPUs over the CPU-based NMFk. Good weak scaling was demonstrated on up to 4096 multi-GPU cluster nodes with approximately 25,000 GPUs when decomposing a dense 340 Terabyte-size matrix and an 11 Exabyte-size sparse matrix of density 10 -6 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct numerical simulations for hybrid rocket boundary layers: Performance modeling and scaling

This paper presents a comprehensive performance and scaling analysis of direct numerical simulations for reacting boundary layers, focusing on slab burner configurations. Using a PETSc-based finite volume CFD framework, the study evaluates the scalability and computational cost of flow, chemistry, and radiation evaluations across 2D and 3D simulations. Polymethyl methacrylate (PMMA) is the fuel with pure O 2 as the oxidizer, modeled using a detailed chemical kinetics mechanism with 113 species and 660 reactions. A ray-tracing-based radiation solver, designed for distributed memory applications, is implemented to model radiation heat transfer. Parallel scalability is analyzed for the coupled flow, chemistry, and radiation heat transfer processes. Weak and strong scaling studies are conducted on up to 15,000 computational ranks, revealing robust performance when flow cells exceed 200 per rank. Chemistry evaluations dominate the computational cost in large 3D simulations, accounting for approximately 40% of the total runtime, while flow processes contribute around 35%, and radiation solver contributions remain below 10% due to reduced evaluation frequencies. GPU accelerated chemistry evaluation, implemented with Zero-RK, demonstrates significant promise, achieving up to a 4x speedup for workloads exceeding 30,000 cells per GPU. However, diminishing returns are observed for smaller workloads due to CPU-GPU communication overhead. This study identifies key challenges, including memory bottlenecks and the effects of domain partitioning on flow scalability, while highlighting the potential of GPU-accelerated chemistry to reduce computational costs. In conclusion, these findings provide realizable run configurations for 2D, 3D, and GPU-accelerated cases, offering insights for optimizing reactive flow solvers.

CFD Scalability↗
DOE OSTI2020

EQSIM—A multidisciplinary framework for fault-to-structure earthquake simulations on exascale computers, part II: Regional simulations of building response

The existing observational database of the regional-scale distribution of strong ground motions and measured building response for major earthquakes continues to be quite sparse. As a result, details of the regional variability and spatial distribution of ground motions, and the corresponding distribution of risk to buildings and other infrastructure, are not comprehensively understood. Utilizing high-performance computing platforms, emerging high-resolution, physics-based ground motion simulations can now resolve frequencies of engineering interest and provide detailed synthetic ground motions at high spatial density. This provides an opportunity for new insight into the distribution of infrastructure seismic demands and risk. In the work presented herein, the EQSIM fault-to-structure computational framework described in a companion paper, McCallen et al., is employed to investigate the regional-scale response of buildings to large earthquakes. A representative M = 7.0 strike-slip event is used to explore the distribution and amplitude of building demand, and comparisons are made between building response computed with fault-to-structure simulations and building response computed with existing measured near-fault earthquake records. New information on the distribution and variability of building response from high-performance parallel simulations is described and analyzed, and favorable first comparisons between building response predicted with both fault-to-structure simulations and real ground motions records are presented.

58 GEOSCIENCES↗
NASA NTRS1968

Thermal effects in face seals.

Incompressible fluid laminar and turbulent flow between face seal rotating parallel surfaces, considering conduction, convection and dissipation influence on thermal distribution

Sneck, H. J.↗
NASA NTRS1973

Topology of induced lunar magnetic fields

Using the asymmetric theory of lunar induction the total and induced magnetic field line structure within the Moon and the diamagnetic cavity were obtained. Total field distributions are shown for orientations of the oscillating interplanetary field parallel, perpendicular and at 45 deg to the cavity axis. Induced field lines are shown only for the orientations of the interplanetary field parallel and orthogonal to the cavity axis. When compared with the field lines derived using the long wavelength limit of spherically symmetric vacuum induction theory, the configurations obtained using the asymmetric theory exhibit significant distortion. For all orientations of the interplanetary field, the field lines are strongly compressed on the sunlit hemisphere because of the confining solar wind pressure at the lunar surface and the exclusion of the field by the lunar core.

Schwartz, K.↗
NASA NTRS1976

Molecular beam mass spectrometer development

An analytical model, based on the kinetics theory of a drifting Maxwellian gas is used to determine the nonequilibrium molecular density distribution within a hemispherical shell open aft with its axis parallel to its velocity. The concept of a molecular shield in terrestrial orbit above 200 km is also analyzed using the kinetic theory of a drifting Maxwellian gas. Data are presented for the components of the gas density within the shield due to the free stream atmosphere, outgassing from the shield and enclosed experiments, and atmospheric gas scattered off a shield orbiter system. A description is given of a FORTRAN program for computating the three dimensional transition flow regime past the space shuttle orbiter that employs the Monte Carlo simulation method to model real flow by some thousands of simulated molecules.

Brock, F. J.↗
NASA NTRS1975

Some correlation of rock exposure ages and regolith dynamics

Exposure age information on lunar rocks and regolith turnover rates are correlated. If plotted in a cumulative fashion, the distribution of spallogenic noble-gas exposure ages is remarkably parallel to the rate at which various fractions of the regolith surface are cratered and/or excavated. It appears that the rate at which lunar rocks are excavated from within the regolith is strongly controlled by the impact environment. Some suggestions for future refinement of regolith dynamics are presented.

Horz, F.↗
NASA NTRS1982

Electron heating within the earth's bow shock

High-temporal-resolution measurements of electron velocity distributions have been obtained for many transits through the earth's bow shock. Within all oblique shocks studied, the maximum of the electron velocity distribution is offset with respect to the ion rest frame parallel to the magnetic field vector and directed downstream. These observations indicate that electron thermalization within the bow shock consists first of a downstream acceleration parallel to the magnetic field vector by the macroscopic shock electric field, followed by beam-driven plasma instabilities.

Feldman, W. C.↗
NASA NTRS1985

Interplanetary magnetic field effects on high latitude ionospheric convection

Relations between the electric field and the electric current in the ionosphere can be established on the basis of a system of mathematical and physical equations provided by the equations of current continuity and Ohm's law. For this reason, much of the synthesis of electric field and plasma velocity data in the F-region is made with the aid of similar data sets derived from field-aligned current and horizontal current measurements. During the past decade, the development of a self-consistent picture of the distribution and behavior of these measurements has proceeded almost in parallel. The present paper is concerned with the picture as it applies to the electric field and plasma drift velocity and its dependence on the interplanetary magnetic field. Attention is given to the southward interplanetary magnetic field and the northward interplanetary magnetic field.

Heelis, R. A.↗
NASA NTRS1986

Ion specific differences in energetic field aligned upflowing ions at 1 earth radius

An ion-upflow event over the northern auroral zone during a small magnetic storm on August 24, 1976 is characterized on the basis of S3-3 satellite ion-mass-spectrometer observations at 0.5-16 keV. The data are presented in tables, graphs, and diagrams and analyzed in detail. The distributions of upward-flowing O(+) and H(+) associated with parallel electric fields below the satellite are shown to have transverse temperatures of about 1 keV and 0.1 keV, respectively, while the corresponding parallel temperatures were about 1 keV and 0.6 keV; the average energy of the O(+) ions is found to be 2-3 times greater than that of the H(+) ions in the acceleration regions. Possible interpretations of these findings are discussed.

Ghielmetti, A. G.↗
NASA NTRS1988

Cosmic-ray hydrodynamics at relativistic shocks

A test particle description of first-order Fermi acceleration of cosmic rays at plane-parallel relativistic shocks is presented in which the pitch angle distribution is obtained explicitly and used to determine the pressure anisotropy or variable Eddington factor as a function of position. The Eddington factor is in turn used to find the hydrodynamical cosmic-ray sound speed, taking the anisotropy into account. For upstream flows speeds above c/sq rt 3, it is shown that the cosmic-ray sound speed exceeds the shock speed in a region just upstream of the shock and decreases montonically with distance further upstream, reaching the value c/sq rt 3 far from the shock. Thus, the shock succeeds in maintaining an anisotropy in the cosmic-ray distribution which is sufficient to keep the flow subsonic with respect to the cosmic-ray sound speed both downstream and in a finite region upstream. The astrophysical implications of these findings are discussed.

Kirk, J. G.↗
NASA NTRS1992

Parallel arrays of Josephson junctions for submillimeter local oscillators

In this paper we discuss the influence of the DC biasing circuit on operation of parallel biased quasioptical Josephson junction oscillator arrays. Because of nonuniform distribution of the DC biasing current along the length of the bias lines, there is a nonuniform distribution of magnetic flux in superconducting loops connecting every two junctions of the array. These DC self-field effects determine the state of the array. We present analysis and time-domain numerical simulations of these states for four biasing configurations. We find conditions for the in-phase states with maximum power output. We compare arrays with small and large inductances and determine the low inductance limit for nearly-in-phase array operation. We show how arrays can be steered in H-plane using the externally applied DC magnetic field.

Pance, Aleksandar↗
NASA NTRS1991

Supercomputing '91; Proceedings of the 4th Annual Conference on High Performance Computing, Albuquerque, NM, Nov. 18-22, 1991

Various papers on supercomputing are presented. The general topics addressed include: program analysis/data dependence, memory access, distributed memory code generation, numerical algorithms, supercomputer benchmarks, latency tolerance, parallel programming, applications, processor design, networks, performance tools, mapping and scheduling, characterization affecting performance, parallelism packaging, computing climate change, combinatorial algorithms, hardware and software performance issues, system issues. (No individual items are abstracted in this volume)

Source record↗