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 109 records · Page 6

Parallel structural optimization with different parallel analysis interfaces

The real benefit of structural optimization techniques is in the application of these techniques to large structures such as full vehicles or full aircraft. For these structures, however, the sequential computer's time and memory requirements prohibit the solutions. With the rapid development of parallel computers, parallel processing of large scale structural optimization problems is achievable. In this paper we discuss the parallel processing of structural optimization problems with parallel structural analysis. Two different types of interface between the optimization and analysis routines are developed and tested.

El-Sayed, Mohamed E. M.↗

Automated Euler and Navier-Stokes Database Generation for a Glide-Back Booster

The past two decades have seen a sustained increase in the use of high fidelity Computational Fluid Dynamics (CFD) in basic research, aircraft design, and the analysis of post-design issues. As the fidelity of a CFD method increases, the number of cases that can be readily and affordably computed greatly diminishes. However, computer speeds now exceed 2 GHz, hundreds of processors are currently available and more affordable, and advances in parallel CFD algorithms scale more readily with large numbers of processors. All of these factors make it feasible to compute thousands of high fidelity cases. However, there still remains the overwhelming task of monitoring the solution process. This paper presents an approach to automate the CFD solution process. A new software tool, AeroDB, is used to compute thousands of Euler and Navier-Stokes solutions for a 2nd generation glide-back booster in one week. The solution process exploits a common job-submission grid environment, the NASA Information Power Grid (IPG), using 13 computers located at 4 different geographical sites. Process automation and web-based access to a MySql database greatly reduces the user workload, removing much of the tedium and tendency for user input errors. The AeroDB framework is shown. The user submits/deletes jobs, monitors AeroDB's progress, and retrieves data and plots via a web portal. Once a job is in the database, a job launcher uses an IPG resource broker to decide which computers are best suited to run the job. Job/code requirements, the number of CPUs free on a remote system, and queue lengths are some of the parameters the broker takes into account. The Globus software provides secure services for user authentication, remote shell execution, and secure file transfers over an open network. AeroDB automatically decides when a job is completed. Currently, the Cart3D unstructured flow solver is used for the Euler equations, and the Overflow structured overset flow solver is used for the Navier-Stokes equations. Other codes can be readily included into the AeroDB framework.

Chaderjian, Neal M.↗

Development of a Robust, High Current, Low Power Field Emission Electron Gun for a Spaceflight Reflectron Time-of-Flight Mass Spectrometer

Carbon materials, including carbon nanotubes (CNTs) and nitrogen-incorporated ultrananocrystalline diamond (N-UNCD), have been of considerable interest for field emission applications for over a decade. In particular, robust field emission materials are compelling for space applications due to the low power consumption and potential for miniaturization. A reflectron time-of-flight mass spectrometer (TOF-MS) under development for in situ measurements on the Moon and other Solar System bodies uses a field emitter to generate ions from gaseous samples, using electron ionization. For these unusual environments, robustness, reliability, and long life are of paramount importance, and to this end, we have explored the field emission properties and lifetime of carbon nanotubes and nitrogen-incorporated ultrananocrystalline diamond (N-UNCD) thin films, the latter developed and patented by Argonne National Laboratory. We will present recent investigations of N-UNCD as a robust field emitter, revealing that this material offers stable performance in high vacuum for up to 1000 hours with threshold voltage for emission of about 3-4 V/lJm and current densities in the range of tens of microA. Optimizing the mass resolution and sensitivity of such a mass spectrometer has also been enabled by a parallel effort to scale up a CNT emitter to an array measuring 2 mm x 40 mm. Through simulation and experiment of the new extended format emitter, we have determined that focusing the electron beam is limited due to the angular spread of the emitted electrons. This dispersion effect can be reduced through modification of the electron gun geometry, but this reduces the current reaching the ionization region. By increasing the transmission efficiency of the electron beam to the anode, we have increased the anode current by two orders of magnitude to realize a corresponding enhancement in instrument sensitivity, at a moderate cost to mass resolution. We will report recent experimental and modeling results to describe the performance of a field emission electron gun as employed in the Volatile Analysis by Pyrolysis of Regolith (VAPoR) TOF-MS prototype.

Southard, Adrian E.↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

Interactions Between Two Existential Threats: COVID-19 and Climate Change

The COVID-19 pandemic and climate change are complex existential threats, unpredictable in many ways and unprecedented in modern times. There are parallels between the scale and scope of their impacts and responses. Understanding shared drivers, coupled vulnerabilities, and criteria for effective responses will help societies worldwide prepare for the simultaneous threats of climate change and future pandemics. We summarize some shared characteristics of COVID-19 and climate change impacts and interventions and discuss key policy implications and recommendations.

Climate change adaptation↗

Automatic differentiation for design sensitivity analysis of structural systems using multiple processors

An automatic differentiation tool (ADIFOR) is incorporated into a finite element based structural analysis program for shape and non-shape design sensitivity analysis of structural systems. The entire analysis and sensitivity procedures are parallelized and vectorized for high performance computation. Small scale examples to verify the accuracy of the proposed program and a medium scale example to demonstrate the parallel vector performance on multiple CRAY C90 processors are included.

Nguyen, Duc T.↗

Anisotropy of Density Fluctuations in the Solar Wind at 1 au

A well-known property of solar wind plasma turbulence is the observed anisotropy of the autocorrelations, or equivalently the spectra, of velocity and magnetic field fluctuations. Here we explore the related but apparently not well-studied issue of the anisotropy of plasma density fluctuations in the energy-containing and inertial ranges of solar wind turbulence. Using 10 yr (1998–2008) of in situ data from the Advanced Composition Explorer mission, we find that for all but the fastest wind category, the density correlation scale is slightly larger in directions quasi-parallel to the large-scale mean magnetic field as compared to quasi-perpendicular directions. The correlation scale in fast wind is consistent with isotropic. The anisotropy as a function of the level of correlation is also explored. We find at small correlation levels, i.e., at energy-containing scales and larger, the density fluctuations are close to isotropy for fast wind, and slightly favor more rapid decorrelation in perpendicular directions for slow and medium winds. At relatively smaller (inertial range) scales where the correlation values are larger, the sense of anisotropy is reversed in all speed ranges, implying a more "slablike" structure, especially prominent in the fast wind samples. We contrast this finding with published results on velocity and magnetic field correlations.

Interplanetary turbulence↗

Characterizing parallel file-access patterns on a large-scale multiprocessor

Rapid increases in the computational speeds of multiprocessors have not been matched by corresponding performance enhancements in the I/O subsystem. To satisfy the large and growing I/O requirements of some parallel scientific applications, we need parallel file systems that can provide high-bandwidth and high-volume data transfer between the I/O subsystem and thousands of processors. Design of such high-performance parallel file systems depends on a thorough grasp of the expected workload. So far there have been no comprehensive usage studies of multiprocessor file systems. Our CHARISMA project intends to fill this void. The first results from our study involve an iPSC/860 at NASA Ames. This paper presents results from a different platform, the CM-5 at the National Center for Supercomputing Applications. The CHARISMA studies are unique because we collect information about every individual read and write request and about the entire mix of applications running on the machines. The results of our trace analysis lead to recommendations for parallel file system design. First the file system should support efficient concurrent access to many files, and I/O requests from many jobs under varying load conditions. Second, it must efficiently manage large files kept open for long periods. Third, it should expect to see small requests predominantly sequential access patterns, application-wide synchronous access, no concurrent file-sharing between jobs appreciable byte and block sharing between processes within jobs, and strong interprocess locality. Finally, the trace data suggest that node-level write caches and collective I/O request interfaces may be useful in certain environments.

Purakayastha, Apratim↗

High temperature oxidation of beta-NiAl

The oxidation of single crystal beta-NiAl has been studied primarily using electron microscopy. Oriented metastable Al2O3 phases form during transient oxidation at 800 C. Specific orientation relationships exist on all metal orientations studied and are a result of the small mismatch along aligned close-packed directions in the cation sublattices of the metal and oxide. Transformation of the metastable Al2O3 phases at 1100 C results in an oxide morphology described as the 'lacey' structure of alpha-Al2O3 scales. This structure results from impingement of oriented patches of alpha-Al2O3 as the transformation initiates and moves radially parallel to the surface. Scale growth occurs by diffusion along high angle grain boundaries. A drastic reduction in oxidation rate accompanies the change in oxide morphology.

Koychak, J. K.↗

A new parallel-vector finite element analysis software on distributed-memory computers

A new parallel-vector finite element analysis software package MPFEA (Massively Parallel-vector Finite Element Analysis) is developed for large-scale structural analysis on massively parallel computers with distributed-memory. MPFEA is designed for parallel generation and assembly of the global finite element stiffness matrices as well as parallel solution of the simultaneous linear equations, since these are often the major time-consuming parts of a finite element analysis. Block-skyline storage scheme along with vector-unrolling techniques are used to enhance the vector performance. Communications among processors are carried out concurrently with arithmetic operations to reduce the total execution time. Numerical results on the Intel iPSC/860 computers (such as the Intel Gamma with 128 processors and the Intel Touchstone Delta with 512 processors) are presented, including an aircraft structure and some very large truss structures, to demonstrate the efficiency and accuracy of MPFEA.

Qin, Jiangning↗

Spontaneous Hot Flow Anomalies at Quasi-Parallel Shocks: 2. Hybrid Simulations

Motivated by recent THEMIS observations, this paper uses 2.5-D electromagnetic hybrid simulations to investigate the formation of Spontaneous Hot Flow Anomalies (SHFA) upstream of quasi-parallel bow shocks during steady solar wind conditions and in the absence of discontinuities. The results show the formation of a large number of structures along and upstream of the quasi-parallel bow shock. Their outer edges exhibit density and magnetic field enhancements, while their cores exhibit drops in density, magnetic field, solar wind velocity and enhancements in ion temperature. Using virtual spacecraft in the simulation, we show that the signatures of these structures in the time series data are very similar to those of SHFAs seen in THEMIS data and conclude that they correspond to SHFAs. Examination of the simulation data shows that SHFAs form as the result of foreshock cavitons interacting with the bow shock. Foreshock cavitons in turn form due to the nonlinear evolution of ULF waves generated by the interaction of the solar wind with the backstreaming ions. Because foreshock cavitons are an inherent part of the shock dissipation process, the formation of SHFAs is also an inherent part of the dissipation process leading to a highly non-uniform plasma in the quasi-parallel magnetosheath including large scale density and magnetic field cavities.

Quasi-Parallel↗

Partitioning of unstructured problems for parallel processing

Many large-scale computational problems are based on unstructured computational domains. Primary examples are unstructured grid calculations based on finite volume methods in computational fluid dynamics, or structural analysis problems based on finite element approximations. The question of how to distribute such unstructured computational domains over a large number of processors in a MIMD machine with distributed memory is addressed. A graph theoretical framework for these problems is established. Based on this framework three decomposition algorithms are introduced. In particular a new decomposition algorithm is discussed, which is based on the computation of an eigenvector of the Laplacian matrix associated with the graph. Numerical comparisons on large-scale two- and three-dimensional problems demonstrate the superiority of the new spectral bisection algorithm.

Simon, H. D.↗

Dispersive and viscous scale lengths in the two-stage ion heating at quasi-parallel collisionless shocks

A simulation of quasi-parallel (QP) collisionless shock waves is used to deduce the dispersive and viscous scale lengths of QP shocks. The findings were found to be consistent with the observational results obtained by Scudder et al. (1984) and Kennel et al. (1984). It is shown that the upstream wavelength is the dispersive scale length which governs the magnetic-field structure of the shock. It is also shown that the two-stage ion heating obtained by Mandt and Kan (1988) is an intrinsic characteristic of the ion heating process in QP shocks.

Mandt, M. E.↗

Hot proton anisotropies and cool proton temperatures in the outer magnetosphere

The hot protons of the outer magnetosphere typically exhibit a temperature anistropy such that T(sub perp)/T(sub parallel) greater than 1, where perpendicular and parallel symbols denote directions relative to the background magnetic field. If this anisotrpy is sufficiently large, the electomagneitc proton cyclotron anistropy instability may be excited. This instability is studied using linear Vlasov theory and one-dimensional hybrid simulations for a homogeneous plasma model representative of conditions in the outer magnetosphere with a hot anisotropic proton component (denoted by subscript h) and a cool, initially isotropic proton component (subscript c). Linear theory yields an instability threshold condition on the hot proton temperature anistropy where as the simulations imply an upper bound on T(sub perp h)/T(sub parallel h); both the threshold and the upper bound have similar scaling with the maximum growth rate gamma (sub m), the parallel beta of the hot component, beta(sub parallel h), and the relative density of the hot component n(sub h)/n(sub e). An anlysis of plasma observations from the Los Alamos magnetospheric plasma analyzer (MPA) in geosynchronous orbits finds that the maximum value of the hot proton temperature anisotropy approximately satisfies the predicted scaling with beta(sub parallel h) and nu(sub h)/n(sub e) and yields the proportionality factor that quantifies this upper bound. The simulations are also used to examine the heating of the cool proton cyclotron instability. The simulations yield a scaling for the dimensionless late-time cool proton average temperature T(sub c)/T(sub parallel h) as (n(sub h)/n(sub e))/beta(sub parallel h exp 0.5). Analysis of MPA data shows that the observed values of T(sub c)/T(sub parallel h) have similar scaling and again yield the proportionality factor which quantifies this relationship.

Gary, S. Peter↗

Hierarchical Parallelism in Finite Difference Analysis of Heat Conduction

Based on the concept of hierarchical parallelism, this research effort resulted in highly efficient parallel solution strategies for very large scale heat conduction problems. Overall, the method of hierarchical parallelism involves the partitioning of thermal models into several substructured levels wherein an optimal balance into various associated bandwidths is achieved. The details are described in this report. Overall, the report is organized into two parts. Part 1 describes the parallel modelling methodology and associated multilevel direct, iterative and mixed solution schemes. Part 2 establishes both the formal and computational properties of the scheme.

Padovan, Joseph↗

A Superposition Technique for Deriving Photon Scattering Statistics in Plane-Parallel Cloudy Atmospheres

Photon transport in a multiple scattering medium is critically dependent on scattering statistics, in particular the average number of scatterings. A superposition technique is derived to accurately determine the average number of scatterings encountered by reflected and transmitted photons within arbitrary layers in plane-parallel, vertically inhomogeneous clouds. As expected, the resulting scattering number profiles are highly dependent on cloud particle absorption and solar/viewing geometry. The technique uses efficient adding and doubling radiative transfer procedures, avoiding traditional time-intensive Monte Carlo methods. Derived superposition formulae are applied to a variety of geometries and cloud models, and selected results are compared with Monte Carlo calculations. Cloud remote sensing techniques that use solar reflectance or transmittance measurements generally assume a homogeneous plane-parallel cloud structure. The scales over which this assumption is relevant, in both the vertical and horizontal, can be obtained from the superposition calculations. Though the emphasis is on photon transport in clouds, the derived technique is applicable to any scattering plane-parallel radiative transfer problem, including arbitrary combinations of cloud, aerosol, and gas layers in the atmosphere.

Platnick, S.↗

Parallel and Distributed Computational Fluid Dynamics: Experimental Results and Challenges

This paper describes several results of parallel and distributed computing using a large scale production flow solver program. A coarse grained parallelization based on clustering of discretization grids combined with partitioning of large grids for load balancing is presented. An assessment is given of its performance on distributed and distributed-shared memory platforms using large scale scientific problems. An experiment with this solver, adapted to a Wide Area Network execution environment is presented. We also give a comparative performance assessment of computation and communication times on both the tightly and loosely-coupled machines.

Djomehri, Mohammad Jahed↗