Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel performance”

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 379 records · Page 21

Explosive Soot Challenge (Final Report)

This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Separation of Lanthanide Isotopes from Mixed Fission Product Samples

The measurement of radioactive fission products from nuclear events has important implications for nuclear data production, environmental monitoring, and nuclear forensics. In a previous paper, the authors reported the optimization of an intra-group lanthanide separation using LN extraction resin from Eichrom Technologies®, Inc. and a nitric acid gradient. In this work, the method was demonstrated for the separation and quantification of multiple short-lived fission product lanthanide isotopes from a fission product sample produced from the thermal irradiation of highly enriched uranium. The separations were performed in parallel in quadruplicate with reproducible results and high decontamination factors for 153 Sm, 156 Eu, and 161 Tb. Based on the results obtained here, the fission yields for 144 Ce, 153 Sm, 156 Eu, and 161 Tb are consistent with published fission yields. This work demonstrates the effectiveness of the separations for the intended application of short-lived lanthanide fission product analysis requiring high decontamination factors.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Scintillating Bubble Chambers for Rare Event Searches

The Scintillating Bubble Chamber (SBC) collaboration is developing liquid-noble bubble chambers for the detection of sub-keV nuclear recoils. These detectors benefit from the electron recoil rejection inherent in moderately-superheated bubble chambers with the addition of energy reconstruction provided from the scintillation signal. The ability to measure low-energy nuclear recoils allows the search for GeV-scale dark matter and the measurement of coherent elastic neutrino-nucleus scattering on argon from MeV-scale reactor antineutrinos. The first physics-scale detector, SBC-LAr10, is in the commissioning phase at Fermilab, where extensive engineering and calibration studies will be performed. In parallel, a functionally identical low-background version, SBC-SNOLAB, is being built for a dark matter search underground at SNOLAB. SBC-SNOLAB, with a 10 kg-yr exposure, will have sensitivity to a dark matter–nucleon cross section of 2×10 –42 cm 2 at 1 GeV/c 2 dark matter mass, and future detectors could reach the boundary of the argon neutrino fog with a tonne-yr exposure. In addition, the deployment of an SBC detector at a nuclear reactor could enable neutrino physics investigations including measurements of the weak mixing angle and searches for sterile neutrinos, the neutrino magnetic moment, and the light Z’ gauge boson.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗

Unpowered to powered failure rate ratio - A key reliability parameter

It is shown that the initial assumption of the ratio of unpowered to powered failure rates can have a strong influence on the design of a modular system intended for space missions. The analysis is performed for parallel systems and for triple modular redundant systems. Parallel systems are shown to be much more sensitive to the unpowered to powered failure rate ratio than the TMR/Spares systems, however, regardless of which standby redundancy technique is considered, the dependence of the system reliability on this ratio increases as the number of standby spares increases.

Taylor, D. S.↗

Miniature triaxial metastable ionization detector for gas chromatographic trace analysis of extraterrestrial volatiles

The present investigation is concerned with a miniature metastable ionization detector featuring an unconventional electrode configuration, whose performance characteristics parallel those of traditional design. The ionization detector is to be incorporated in a flight gas chromatograph (GC) for use in the Space Shuttle. The design of the detector is discussed, taking into account studies which verified the sensitivity of the detector. The triaxial design of the detector is compared with a flat-plate style. The obtained results show that the principal goal of developing a miniature, highly sensitive ionization detector for flight applications was achieved. Improved fabrication techniques will utilize glass-to-metal seals and brazing procedures.

Woeller, F. H.↗

A 406 MHz geostationary satellite SARSAT experiment

A 406 MHz SARSAT experiment using the GOES system to demonstrate the utility of integrating a near-instantaneous distress alerting capability achievable by the geostationary satellites with the location capability of the low altitude, near-polar orbiting satellite system is described. Tests include the evaluation of technical performance of a variety of maritime, aviation, and land-based SARSAT beacons under operational environment conditions as well as a demonstration of system effectiveness in pseudo real distress incident situations. The technical performance tests parallel the CCIR coordinated trials program conducted at L-band frequencies using the INMARSAT satellites. Tests, including schedules, are outlined.

Friedman, M. L.↗

Anomalous transport by Kelvin-Helmholtz instabilities

A simulation of magnetohydrodynamic Kelvin-Helmholtz instabilities has been performed for parallel and transverse configurations, modeling high latitude (or downstream flanks) and dayside low latitude magnetospheric boundaries. In the parallel configuration, a super-Alfvenic and transsonic shear flow develops into small eddies, which strongly compresses, twists, and hence amplifies the magnetic field by the dynamo action with an amplification factor of M(A)/2. In the nonlinear stage, however large the initial Alfven mach number M(A) may be, the magnetic field amplified and twisted by the hydromagnetic flow vortices reacts back upon the flow evolution, and the flow vortices cascade into smaller structures. In the transverse configuration, the instability leads to the formation of a fast shock discontinuity from an initially subfast shear flow. Anomalous tangential stress by the instability in the transverse configuration reaches 1 percent of the magnetosheath momentum flux, but for the parallel configuration, the anomalous transport is 2-3 times larger than the anomalous transport in the transverse configuration. The anomalous transport for both configurations satisfies the requirement of the viscous-like interaction at the magnetospheric boundary.

Miura, A.↗

Negative base encoding in optical linear algebra processors

In the digital multiplication by analog convolution algorithm, the bits of two encoded numbers are convolved to form the product of the two numbers in mixed binary representation; this output can be easily converted to binary. Attention is presently given to negative base encoding, treating base -2 initially, and then showing that the negative base system can be readily extended to any radix. In general, negative base encoding in optical linear algebra processors represents a more efficient technique than either sign magnitude or 2's complement encoding, when the additions of digitally encoded products are performed in parallel.

Perlee, C.↗

Galerkin/Runge-Kutta discretizations for semilinear parabolic equations

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for semilinear parabolic initial boundary value problems. Unlike any classical counterpart, this class offers arbitrarily high, optimal order convergence. In support of this claim, error estimates are proved, and computational results are presented. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗

Galerkin/Runge-Kutta discretizations for parabolic equations with time dependent coefficients

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗

Euler flutter analysis of airfoils using unstructured dynamic meshes

Modifications to a two-dimensional unsteady Euler code for the aeroelastic analysis of airfoils are described. The modifications involve including the structural equations of motion and their simultaneous time-integration with the governing flow equations. A novel aspect of the capability is that the solutions are obtained using unstructured grids made up of triangles. Comparisons are made with parallel calculations performed using linear theory and a structured grid Euler code to assess the accuracy of the unstructured grid Euler results. Results are presented for a flat plate airfoil and the NACA 0012 airfoil to demonstrate applications of the Euler code for generalized force computations and aeroelastic analysis. In these comparisons, two different finite-volume discretizations of the Euler equations on unstructured meshes were employed. Sensitivity of the Euler results to changes in numerical parameters were also investigated.

Rausch, Russ D.↗

Galerkin/Runge-Kutta discretizations for parabolic equations with time-dependent coefficients

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗

FFTs in external or hierarchical memory

A description is given of advanced techniques for computing an ordered FFT on a computer with external or hierarchical memory. These algorithms (1) require as few as two passes through the external data set, (2) use strictly unit stride, long vector transfers between main memory and external storage, (3) require only a modest amount of scratch space in main memory, and (4) are well suited for vector and parallel computation. Performance figures are included for implementations of some of these algorithms on Cray supercomputers. Of interest is the fact that a main memory version outperforms the current Cray library FFT routines on the Cray-2, the Cray X-MP, and the Cray Y-MP systems. Using all eight processors on the Cray Y-MP, this main memory routine runs at nearly 2 Gflops.

Bailey, David H.↗

Theoretical and experimental investigation of the aeroelastic stability of an advanced bearingless rotor in hover and forward flight

The aeroelastic stability of a shaft-fixed, 1/8th Froude scaled bearingless rotor was investigated in a series of wind tunnel experiments simulating a wide range of operating conditions. A finite element formulation was used to perform a parallel theoretical analysis, with the goal of determining whether a bearingless rotor system could be made aeroelastically stable without the incorporation of auxilliary dampers. A quick estimate of lag mode damping was provided by a refined moving-block analysis implemented in real time which predicted similar damping values. Model rotor and blade properties were also determined, and these properties were used as inputs for a newly refined bearingless rotor analysis. Predicted results were compared with experimental results in hover and forward flight. Results indicated that soft pitch link stiffness increases pitch-lag coupling and stabilizes lag mode stability in hover and at low advance ratios, but destabilizes at higher advance ratios.

Wang, James M.↗

Ferroelectric/Optoelectronic Memory/Processor

Proposed hybrid optoelectronic nonvolatile analog memory and data processor comprises planar array of microscopic photosensitive ferroelectric capacitors performing massively parallel analog computations. Processors overcome electronic crosstalk and limitations on number of input/output contacts inherent in electronic implementations of large interconnection arrays. Used in general optical computing, recognition of patterns, and artificial neural networks.

Thakoor, Sarita↗

Self-checking self-repairing computer nodes using the mirror processor

Circuitry added to fault-tolerant systems for concurrent error deduction usually reduces performance. Using a technique called micro rollback, it is possible to eliminate most of the performance penalty of concurrent error detection. Error detection is performed in parallel with intermodule communication, and erroneous state changes are later undone. The author reports on the design and implementation of a VLSI RISC microprocessor, called the Mirror Processor (MP), which is capable of micro rollback. In order to achieve concurrent error detection, two MP chips operate in lockstep, comparing external signals and a signature of internal signals every clock cycle. If a mismatch is detected, both processors roll back to the beginning of the cycle when the error occurred. In some cases the erroneous state is corrected by copying a value from the fault-free processor to the faulty processor. The architecture, microarchitecture, and VLSI implementation of the MP, emphasizing its error-detection, error-recovery, and self-diagnosis capabilities, are described.

Tamir, Yuval↗