Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parallelism”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Experimental Study of an Easily Controlled Ultra-High-Resolution Pixel-Matched Parallel-Hole Collimator with a Small Cadmium Zinc Telluride Pixelated Gamma Camera System

The aim of this study was to develop an easily controlled, ultra-high-resolution, tungsten parallel-hole collimator based on a small pixelated gamma camera system. A small cadmium zinc telluride (CZT) pixelated semiconductor detector (eValuator-2500 detector [eV product, Saxonburg, PA]) was evaluated. This detector is composed of an array of 51.2 × 0.8 × 3-mm{sup 3} individual CZT crystal elements. The ultra-high-resolution, pixel-matched, parallel-hole collimators consisted of six layers, with the same between the hole and pixel size. The basic characteristics of the imaging system, such as sensitivity and spatial resolution, was measured using a {sup 57}Co point source. The measured averages of sensitivity and spatial resolution varied depending on the septal heights of the ultra-high-resolution parallel-hole collimator and source-to-collimator distances. When the 30-mm septal height was at 1-cm source-to-collimator distance, the spatial resolution was approximately 0.85 mm. Using 5-mm septal height, over 0.3 cps/kBq sensitivity was achieved. One advantage of our system is the use of stacked collimators that can select the best combination of system sensitivity and spatial resolution. Our results demonstrated that the developed CZT–pixelated gamma camera system using an ultra-high-resolution parallel-hole collimator of various collimator geometric designs has potential as an effective instrument.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Parallel Electrochemical Treatment System and Application for Identifying Acid-Stable Oxygen Evolution Electrocatalysts

There are many energy technologies require electrochemical stability or preactivation of functional materials. Due to the long experiment duration required for either electrochemical preactivation or evaluation of operational stability, parallel screening is required to enable high throughput experimentation. We found that imposing operational electrochemical conditions to a library of materials in parallel creates several opportunities for experimental artifacts. We discuss the electrochemical engineering principles and operational parameters that mitigate artifacts int he parallel electrochemical treatment system. We also demonstrate the effects of resistive losses within the planar working electrode through a combination of finite element modeling and illustrative experiments. Operation of the parallel-plate, membrane-separated electrochemical treatment system is demonstrated by exposing a composition library of mixed metal oxides to oxygen evolution conditions in 1M sulfuric acid for 2h. This application is particularly important because the electrolysis and photoelectrolysis of water are promising future energy technologies inhibited by the lack of highly active, acid-stable catalysts containing only earth abundant elements.

14 SOLAR ENERGY↗
DOE OSTI2012

Measurements of parallel electron velocity distributions using whistler wave absorption

We describe a diagnostic to measure the parallel electron velocity distribution in a magnetized plasma that is overdense ({omega}{sub pe} > {omega}{sub ce}). This technique utilizes resonant absorption of whistler waves by electrons with velocities parallel to a background magnetic field. The whistler waves were launched and received by a pair of dipole antennas immersed in a cylindrical discharge plasma at two positions along an axial background magnetic field. The whistler wave frequency was swept from somewhat below and up to the electron cyclotron frequency {omega}{sub ce}. As the frequency was swept, the wave was resonantly absorbed by the part of the electron phase space density which was Doppler shifted into resonance according to the relation {omega}-k{sub ||v||} = {omega}{sub ce}. The measured absorption is directly related to the reduced parallel electron distribution function integrated along the wave trajectory. The background theory and initial results from this diagnostic are presented here. Though this diagnostic is best suited to detect tail populations of the parallel electron distribution function, these first results show that this diagnostic is also rather successful in measuring the bulk plasma density and temperature both during the plasma discharge and into the afterglow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

Numerical investigation of edge plasma phenomena in an enhanced D-alpha discharge at Alcator C-Mod: Parallel heat flux and quasi-coherent edge oscillations

Reduced-model scrape-off layer turbulence (SOLT) simulations of an enhanced D-alpha (EDA) H-mode shot observed in the Alcator C-Mod tokamak were conducted to compare with observed variations in the scrape-off-layer (SOL) width of the parallel heat flux profile. In particular, the role of the competition between sheath- and conduction-limited parallel heat fluxes in determining that width was studied for the turbulent SOL plasma that emerged from the simulations. The SOL width decreases with increasing input power and with increasing separatrix temperature in both the experiment and the simulation, consistent with the strong temperature dependence of the parallel heat flux in balance with the perpendicular transport by turbulence and blobs. The particularly strong temperature dependence observed in the case analyzed is attributed to the fact that these simulations produce SOL plasmas which are in the conduction-limited regime for the parallel heat flux. A persistent quasi-coherent (QC) mode dominates the SOLT simulations and bears considerable resemblance to the QC mode observed in C-Mod EDA operation. The SOLT QC mode consists of nonlinearly saturated wave-fronts located just inside the separatrix that are convected poloidally by the mean flow, continuously transporting particles and energy and intermittently emitting blobs into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Effects of parallel electron dynamics on plasma blob transport

The 3D effects on sheath connected plasma blobs that result from parallel electron dynamics are studied by allowing for the variation of blob density and potential along the magnetic field line and using collisional Ohm's law to model the parallel current density. The parallel current density from linear sheath theory, typically used in the 2D model, is implemented as parallel boundary conditions. This model includes electrostatic 3D effects, such as resistive drift waves and blob spinning, while retaining all of the fundamental 2D physics of sheath connected plasma blobs. If the growth time of unstable drift waves is comparable to the 2D advection time scale of the blob, then the blob's density gradient will be depleted resulting in a much more diffusive blob with little radial motion. Furthermore, blob profiles that are initially varying along the field line drive the potential to a Boltzmann relation that spins the blob and thereby acts as an addition sink of the 2D potential. Basic dimensionless parameters are presented to estimate the relative importance of these two 3D effects. The deviation of blob dynamics from that predicted by 2D theory in the appropriate limits of these parameters is demonstrated by a direct comparison of 2D and 3D seeded blob simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Electron cyclotron current drive modelling with parallel momentum correction for tokamaks and stellarators

The electron cyclotron current drive (ECCD) is described by the solution of a drift-kinetic equation (DKE) with the quasi-linear heating term which is highly localised in phase space. Within the adjoint approach, an equivalent DKE must be solved without this highly localised source, and the ECCD is described by a straightforward convolution in phase space. The parallel momentum correction technique [Maassberg et al., Phys. Plasmas 16, 072504 (2009)] reduces the problem to the solution of a simplified mono-energetic DKE. The mono-energetic distribution functions are calculated by the DKES code [W. I. van Rij and S. P. Hirshman, Phys. Fluids B 1, 563 (1989)]. Parallel momentum correction, requiring only an energy-dependent weighting factor, is applied to these distribution functions allowing for the direct estimation of the EC-driven current. For small yet finite collisionalities, a rather simple model for ECCD is presented. In this approach, the symmetric portion of the electron distribution function with respect to v{sub Parallel-To} (the Ware pinch contribution) is neglected and an 'off-set' only in the passing particle domain adds to the collisionless (anti-symmetric) distribution function. For this approximation, only the mono-energetic parallel conductivity coefficient as function of the collisionality is needed. The impact of small yet finite collisionalities on ECCD is shown for X2- and O2-scenarios at W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Parallel computation of multigroup reactivity coefficient using iterative method

One of the research activities to support the commercial radioisotope production program is a safety research target irradiation FPM (Fission Product Molybdenum). FPM targets form a tube made of stainless steel in which the nuclear degrees of superimposed high-enriched uranium. FPM irradiation tube is intended to obtain fission. The fission material widely used in the form of kits in the world of nuclear medicine. Irradiation FPM tube reactor core would interfere with performance. One of the disorders comes from changes in flux or reactivity. It is necessary to study a method for calculating safety terrace ongoing configuration changes during the life of the reactor, making the code faster became an absolute necessity. Neutron safety margin for the research reactor can be reused without modification to the calculation of the reactivity of the reactor, so that is an advantage of using perturbation method. The criticality and flux in multigroup diffusion model was calculate at various irradiation positions in some uranium content. This model has a complex computation. Several parallel algorithms with iterative method have been developed for the sparse and big matrix solution. The Black-Red Gauss Seidel Iteration and the power iteration parallel method can be used to solve multigroup diffusion equation system and calculated the criticality and reactivity coeficient. This research was developed code for reactivity calculation which used one of safety analysis with parallel processing. It can be done more quickly and efficiently by utilizing the parallel processing in the multicore computer. This code was applied for the safety limits calculation of irradiated targets FPM with increment Uranium.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Massively parallel implementations of coupled-cluster methods for electron spin resonance spectra. I. Isotropic hyperfine coupling tensors in large radicals

Coupled cluster (CC) methods provide highly accurate predictions of molecular properties, but their high computational cost has precluded their routine application to large systems. Fortunately, recent computational developments in the ACES III program by the Bartlett group [the OED/ERD atomic integral package, the super instruction processor, and the super instruction architecture language] permit overcoming that limitation by providing a framework for massively parallel CC implementations. In that scheme, we are further extending those parallel CC efforts to systematically predict the three main electron spin resonance (ESR) tensors (A-, g-, and D-tensors) to be reported in a series of papers. In this paper inaugurating that series, we report our new ACES III parallel capabilities that calculate isotropic hyperfine coupling constants in 38 neutral, cationic, and anionic radicals that include the {sup 11}B, {sup 17}O, {sup 9}Be, {sup 19}F, {sup 1}H, {sup 13}C, {sup 35}Cl, {sup 33}S,{sup 14}N, {sup 31}P, and {sup 67}Zn nuclei. Present parallel calculations are conducted at the Hartree-Fock (HF), second-order many-body perturbation theory [MBPT(2)], CC singles and doubles (CCSD), and CCSD with perturbative triples [CCSD(T)] levels using Roos augmented double- and triple-zeta atomic natural orbitals basis sets. HF results consistently overestimate isotropic hyperfine coupling constants. However, inclusion of electron correlation effects in the simplest way via MBPT(2) provides significant improvements in the predictions, but not without occasional failures. In contrast, CCSD results are consistently in very good agreement with experimental results. Inclusion of perturbative triples to CCSD via CCSD(T) leads to small improvements in the predictions, which might not compensate for the extra computational effort at a non-iterative N{sup 7}-scaling in CCSD(T). The importance of these accurate computations of isotropic hyperfine coupling constants to elucidate experimental ESR spectra, to interpret spin-density distributions, and to characterize and identify radical species is illustrated with our results from large organic radicals. Those include species relevant for organic chemistry, petroleum industry, and biochemistry, such as the cyclo-hexyl, 1-adamatyl, and Zn-porphycene anion radicals, inter alia.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Quantum transport through the system of parallel quantum dots with Majorana bound states

We study the tunneling transport properties through a system of parallel quantum dots which are coupled to Majorana bound states (MBSs). The conductance and spectral function are computed using the retarded Green's function method based on the equation of motion. The conductance of the system is 2e{sup 2}/h at zero Fermi energy and is robust against the coupling between the MBSs and the quantum dots. The dependence of the Fermi energy on the spectral function is different for the first dot (dot1) than for the second dot (dot2) with fixed dot2-MBSs coupling. The influence of the Majorana bound states on the spectral function was studied for the series and parallel configurations of the system. It was found that when the configuration is in series, the Majorana bound states play an important role, resulting in a spectral function with three peaks. However, the spectral function shows two peaks when the system is in a parallel configuration. The zero Fermi energy spectral function is always 1/2 not only in series but also in the parallel configuration and robust against the coupling between the MBSs and the quantum dots. The phase diagram of the Fermi energy versus the quantum dot energy levels was also investigated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Evolution of Kelvin-Helmholtz instability at Venus in the presence of the parallel magnetic field

Two-dimensional MHD simulations were performed to study the evolution of the Kelvin-Helmholtz (KH) instability at the Venusian ionopause in response to the strong flow shear in presence of the in-plane magnetic field parallel to the flow direction. The physical behavior of the KH instability as well as the triggering and occurrence conditions for highly rolled-up vortices are characterized through several physical parameters, including Alfvén Mach number on the upper side of the layer, the density ratio, and the ratio of parallel magnetic fields between two sides of the layer. Using these parameters, the simulations show that both the high density ratio and the parallel magnetic field component across the boundary layer play a role of stabilizing the instability. In the high density ratio case, the amount of total magnetic energy in the final quasi-steady status is much more than that in the initial status, which is clearly different from the case with low density ratio. We particularly investigate the nonlinear development of the case that has a high density ratio and uniform magnetic field. Before the instability saturation, a single magnetic island is formed and evolves into two quasi-steady islands in the non-linear phase. A quasi-steady pattern eventually forms and is embedded within a uniform magnetic field and a broadened boundary layer. The estimation of loss rates of ions from Venus indicates that the stabilizing effect of the parallel magnetic field component on the KH instability becomes strong in the case of high density ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

DGDFT: A massively parallel method for large scale density functional theory calculations

We describe a massively parallel implementation of the recently developed discontinuous Galerkin density functional theory (DGDFT) method, for efficient large-scale Kohn-Sham DFT based electronic structure calculations. The DGDFT method uses adaptive local basis (ALB) functions generated on-the-fly during the self-consistent field iteration to represent the solution to the Kohn-Sham equations. The use of the ALB set provides a systematic way to improve the accuracy of the approximation. By using the pole expansion and selected inversion technique to compute electron density, energy, and atomic forces, we can make the computational complexity of DGDFT scale at most quadratically with respect to the number of electrons for both insulating and metallic systems. We show that for the two-dimensional (2D) phosphorene systems studied here, using 37 basis functions per atom allows us to reach an accuracy level of 1.3 × 10{sup −4} Hartree/atom in terms of the error of energy and 6.2 × 10{sup −4} Hartree/bohr in terms of the error of atomic force, respectively. DGDFT can achieve 80% parallel efficiency on 128,000 high performance computing cores when it is used to study the electronic structure of 2D phosphorene systems with 3500-14 000 atoms. This high parallel efficiency results from a two-level parallelization scheme that we will describe in detail.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Influence of equilibrium shear flow in the parallel magnetic direction on edge localized mode crash

The influence of the parallel shear flow on the evolution of peeling-ballooning (P-B) modes is studied with the BOUT++ four-field code in this paper. The parallel shear flow has different effects in linear simulation and nonlinear simulation. In the linear simulations, the growth rate of edge localized mode (ELM) can be increased by Kelvin-Helmholtz term, which can be caused by the parallel shear flow. In the nonlinear simulations, the results accord with the linear simulations in the linear phase. However, the ELM size is reduced by the parallel shear flow in the beginning of the turbulence phase, which is recognized as the P-B filaments' structure. Then during the turbulence phase, the ELM size is decreased by the shear flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Reduction of blob-filament radial propagation by parallel variation of flows: Analysis of a gyrokinetic simulation

Data from the XGC1 gyrokinetic simulation is analyzed to understand the three-dimensional spatial structure and the radial propagation of blob-filaments generated by quasisteady turbulence in the tokamak edge pedestal and scrape-off layer plasma. Spontaneous toroidal flows vary in the poloidal direction and shear the filaments within a flux surface resulting in a structure that varies in the parallel direction. Here, this parallel structure allows the curvature and grad-B induced polarization charge density to be shorted out via parallel electron motion. As a result, it is found that the blob-filament radial velocity is significantly reduced from estimates which neglect parallel electron kinetics, broadly consistent with experimental observations. Conditions for when this charge shorting effect tends to dominate blob dynamics are derived and compared with the simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗
DOE OSTI2024

Kinetic study of shock formation and particle acceleration in laser-driven quasi-parallel magnetized collisionless shocks

Quasi-parallel magnetized collisionless shocks are believed to be one of the most efficient accelerators in the universe. Compared to quasi-perpendicular shocks, quasi-parallel shocks are more difficult to form in the laboratory and to simulate because of their large spatial scales and long formation times. Our two-dimensional particle-in-cell simulations show that the early stages of quasi-parallel shock formation are achievable in experiments planned for the National Ignition Facility and that particles accelerated by diffusive shock acceleration (DSA) are expected to be observable in the experiment. Repetitive ion acceleration by crossings of the shock front, a key feature of DSA, is seen in the simulations. Other characteristic features of quasi-parallel shocks such as upstream wave excitation by energetic ions are also observed, and energy partition between the ions and the electrons in the downstream of the shock is briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Parallel interior-point solver for block-structured nonlinear programs on SIMD/GPU architectures

Here, we investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses SIMD/GPU accelerators locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations. We demonstrate our method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the reduction method is 50x faster than the sparse linear solver HSL MA57 running in serial on the CPU, and 6x faster than Pardiso running in parallel on CPU on the same number of processes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A PARALLEL MONTE CARLO CODE FOR SIMULATING COLLISIONAL N-BODY SYSTEMS

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N {approx} 10{sup 7} particles. Our code is based on the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures and the introduction of a parallel random number generation scheme as well as a parallel sorting algorithm required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. Our implementation uses the Message Passing Interface library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude from 10{sup 5} to 10{sup 7}. We find that our results are in good agreement with self-similar core-collapse solutions, and the core-collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within {approx}< 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N = 10{sup 5}, 128 for N = 10{sup 6} and 256 for N = 10{sup 7}. The runtime reaches saturation with the addition of processors beyond these limits, which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60 Multiplication-Sign , 100 Multiplication-Sign , and 220 Multiplication-Sign , respectively.

79 ASTRONOMY AND ASTROPHYSICS↗