Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallelization Strategy”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Simple adaptive control system design for a quadrotor with an internal PFC

The paper deals with an adaptive control system design problem for a four rotor helicopter or quadrotor. A simple adaptive control design scheme with a parallel feedforward compensator (PFC) in the internal loop of the considered quadrotor will be proposed based on the backstepping strategy. As is well known, the backstepping control strategy is one of the advanced control strategy for nonlinear systems. However, the control algorithm will become complex if the system has higher order relative degrees. We will show that one can skip some design steps of the backstepping method by introducing a PFC in the inner loop of the considered quadrotor, so that the structure of the obtained controller will be simplified and a high gain based adaptive feedback control system will be designed. The effectiveness of the proposed method will be confirmed through numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Efficient Delaunay Tessellation through K-D Tree Decomposition

Delaunay tessellations are fundamental data structures in computational geometry. They are important in data analysis, where they can represent the geometry of a point set or approximate its density. The algorithms for computing these tessellations at scale perform poorly when the input data is unbalanced. We investigate the use of k-d trees to evenly distribute points among processes and compare two strategies for picking split points between domain regions. Because resulting point distributions no longer satisfy the assumptions of existing parallel Delaunay algorithms, we develop a new parallel algorithm that adapts to its input and prove its correctness. We evaluate the new algorithm using two late-stage cosmology datasets. The new running times are up to 50 times faster using k-d tree compared with regular grid decomposition. Moreover, in the unbalanced data sets, decomposing the domain into a k-d tree is up to five times faster than decomposing it into a regular grid.

Morozov, Dmitriy↗
DOE OSTI2015

SU-E-T-33: A Feasibility-Seeking Algorithm Applied to Planning of Intensity Modulated Proton Therapy: A Proof of Principle Study

Purpose: To investigate the applicability of feasibility-seeking cyclic orthogonal projections to the field of intensity modulated proton therapy (IMPT) inverse planning. Feasibility of constraints only, as opposed to optimization of a merit function, is less demanding algorithmically and holds a promise of parallel computations capability with non-cyclic orthogonal projections algorithms such as string-averaging or block-iterative strategies. Methods: A virtual 2D geometry was designed containing a C-shaped planning target volume (PTV) surrounding an organ at risk (OAR). The geometry was pixelized into 1 mm pixels. Four beams containing a subset of proton pencil beams were simulated in Geant4 to provide the system matrix A whose elements a-ij correspond to the dose delivered to pixel i by a unit intensity pencil beam j. A cyclic orthogonal projections algorithm was applied with the goal of finding a pencil beam intensity distribution that would meet the following dose requirements: D-OAR < 54 Gy and 57 Gy < D-PTV < 64.2 Gy. The cyclic algorithm was based on the concept of orthogonal projections onto half-spaces according to the Agmon-Motzkin-Schoenberg algorithm, also known as ‘ART for inequalities’. Results: The cyclic orthogonal projections algorithm resulted in less than 5% of the PTV pixels and less than 1% of OAR pixels violating their dose constraints, respectively. Because of the abutting OAR-PTV geometry and the realistic modelling of the pencil beam penumbra, complete satisfaction of the dose objectives was not achieved, although this would be a clinically acceptable plan for a meningioma abutting the brainstem, for example. Conclusion: The cyclic orthogonal projections algorithm was demonstrated to be an effective tool for inverse IMPT planning in the 2D test geometry described. We plan to further develop this linear algorithm to be capable of incorporating dose-volume constraints into the feasibility-seeking algorithm.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Comparison of Some RANS Solvers

We will take a look at solving the Reynolds-averaged Navier-Stokes (RANS) equations that are encountered in the context of wind farm performance simulations and optimizations. We will compare some of the more popular ways to solve these equations with a focus on using iterative solvers for the linear solve. We will compare their performance and reliability to a direct solve as we scale the problem both by adding more and parallel resources and by increasing the size of the domain, both in two and three dimensions. There are many strategies that can be applied to solving the RANS equations, some are very efficient, while others are very insensitive to, for example, the Reynolds number. The first contender we will consider is the Pressure-Convection-Diffusion (PCD) preconditioner. Early results suggest that PCD is indeed a very efficient solver, in particular in two dimensions, as long as the Reynold's number remains small. Next we will try to reorder our degrees of freedom such that we can use GMRES with ILU for our linear solve. Another popular choice we will consider for solving the RANS equations is SIMPLE (and its derivatives). For all of our implementations we make use of either FEniCS or Firedrake, basing our work on both existing implementations of some of these solvers while also writing new extensions for others.

CFD↗
DOE OSTI2021

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

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

Work stealing for GPU-accelerated parallel programs in a global address space framework

Task parallelism is an attractive approach to automatically load balance the computation in a parallel system and adapt to dynamism exhibited by parallel systems. Exploiting task parallelism through work stealing has been extensively studied in shared and distributed-memory contexts. In this paper, we study the design of a system that uses work stealing for dynamic load balancing of task-parallel programs executed on hybrid distributed-memory CPU-graphics processing unit (GPU) systems in a global-address space framework. We take into account the unique nature of the accelerator model employed by GPUs, the significant performance difference between GPU and CPU execution as a function of problem size, and the distinct CPU and GPU memory domains. We consider various alternatives in designing a distributed work stealing algorithm for CPU-GPU systems, while taking into account the impact of task distribution and data movement overheads. These strategies are evaluated using microbenchmarks that capture various execution configurations as well as the state-of-the-art CCSD(T) application module from the computational chemistry domain

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Parallel Multi-cycle LES of an Optical Pent-roof DISI Engine Under Motored Operating Conditions

The use of Large-eddy Simulations (LES) has increased due to their ability to resolve the turbulent fluctuations of engine flows and capture the resulting cycle-to-cycle variability. One drawback of LES, however, is the requirement to run multiple engine cycles to obtain the necessary cycle statistics for full validation. The standard method to obtain the cycles by running a single simulation through many engine cycles sequentially can take a long time to complete. Recently, a new strategy has been proposed by our research group to reduce the amount of time necessary to simulate the many engine cycles by running individual engine cycle simulations in parallel. With modern large computing systems this has the potential to reduce the amount of time necessary for a full set of simulated engine cycles to finish by up to an order of magnitude. In this paper, the Parallel Perturbation Methodology (PPM) is used to simulate up to 35 engine cycles of an optically accessible, pent-roof Directinjection Spark-ignition (DISI) engine at two different motored engine operating conditions, one throttled and one un-throttled. Comparisons are made against corresponding sequential-cycle simulations to verify the similarity of results using either methodology. Mean results from the PPM approach are very similar to sequential-cycle results with less than 0.5% difference in pressure and a magnitude structure index (MSI) of 0.95. Differences in cycle-to-cycle variability (CCV) predictions are larger, but close to the statistical uncertainty in the measurement for the number of cycles simulated. PPM LES results were also compared against experimental data. Mean quantities such as pressure or mean velocities were typically matched to within 5- 10%. Pressure CCVs were under-predicted, mostly due to the lack of any perturbations in the pressure boundary conditions between cycles. Velocity CCVs for the simulations had the same average magnitude as experiments, but the experimental data showed greater spatial variation in the root-mean-square (RMS). Conversely, circular standard deviation results showed greater repeatability of the flow directionality and swirl vortex positioning than the simulations.

I.C. Engines↗
DOE OSTI2022

Strategies for conformal REBCO windings

A high-field winding can be fabricated from a cable of non-insulated REBCO tapes stacked face-to-face without twisting. If the cable is oriented within each turn of a winding so that the tape face is closely parallel to the magnetic field at its location, the supercurrent capacity of that cable is enhanced ~3x greater than in a transverse or twisting orientation. This concept for a conformal winding was presented in a previous paper pertinent to the body winding of a REBCO based high-field dipole. Strategies are presented and simulated for how the same orientation can be sustained in the flared ends of a high-field hybrid dipole.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Bioresorbable Primary Battery Anodes Built on Core–Double-Shell Zinc Microparticle Networks

Bioresorbable implantable electronics require power sources that are also bioresorbable with controllable electrical output and lifetime. In this paper, we report a bioresorbable zinc primary battery anode filament based on a zinc microparticle (MP) network coated with chitosan and Al 2 O 3 double shells. When discharged in 0.9% NaCl saline, a Zn MP filament with a 0.17 × 2 mm 2 cross-sectional area exhibited a stable voltage output of 0.55 V at a current of 0.01 mA. Covered by chitosan and Al 2 O 3 double shells, the zinc MP filament exhibited a directional dissolution behavior with a tunable lifetime approximately linear to its length. A stable 200 h discharging time was achieved with a 15 mm Zn MP filament. The maximum output power was found to be 12 μW at 0.03 mA for one filament. The linearity relationship between the current output and the filament cross-sectional area suggested a facile strategy to raise the power output at constant discharging voltage. The filaments could also be connected in series and in parallel to boost its overall voltage and current output, demonstrating their excellent integration capability. Furthermore, this work presents a promising pathway toward bioresorbable transient batteries with controllable lifetime and power output, demonstrating a great potential for powering transient implantable biomedical devices.

25 ENERGY STORAGE↗
DOE OSTI2012

THE PANCHROMATIC HUBBLE ANDROMEDA TREASURY

The Panchromatic Hubble Andromeda Treasury is an ongoing Hubble Space Telescope Multi-Cycle Treasury program to image {approx}1/3 of M31's star-forming disk in six filters, spanning from the ultraviolet (UV) to the near-infrared (NIR). We use the Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) to resolve the galaxy into millions of individual stars with projected radii from 0 to 20 kpc. The full survey will cover a contiguous 0.5 deg{sup 2}area in 828 orbits. Imaging is being obtained in the F275W and F336W filters on the WFC3/UVIS camera, F475W and F814W on ACS/WFC, and F110W and F160W on WFC3/IR. The resulting wavelength coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for most spectral types. The data produce photometry with a signal-to-noise ratio of 4 at m{sub F275W} = 25.1, m{sub F336W} = 24.9, m{sub F475W} = 27.9, m{sub F814W} = 27.1, m{sub F110W} = 25.5, and m{sub F160W} = 24.6 for single pointings in the uncrowded outer disk; in the inner disk, however, the optical and NIR data are crowding limited, and the deepest reliable magnitudes are up to 5 mag brighter. Observations are carried out in two orbits per pointing, split between WFC3/UVIS and WFC3/IR cameras in primary mode, with ACS/WFC run in parallel. All pointings are dithered to produce Nyquist-sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, astrometry, and data products available for the survey, along with extensive testing of photometric stability, crowding errors, spatially dependent photometric biases, and telescope pointing control. We also report on initial fits to the structure of M31's disk, derived from the density of red giant branch stars, in a way that is independent of assumed mass-to-light ratios and is robust to variations in dust extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star formation, but is instead a dynamical structure containing a significant overdensity of stars with ages >1 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

AdaParse

SF-25-127AdaParse (Adaptive Parallel PDF Parsing and Resource Scaling Engine) enables scalable, high-accuracy PDF parsing. AdaParse is a data-driven strategy that assigns an appropriate parser to each document, offering high accuracy for any computational budget. Moreover, it offers a workflow of various PDF parsing software that includes extraction tools: PyMuPDF, pypdf traditional OCR: Tesseract, modern OCR (e.g., Vision Transformers): Nougat and Marker

Siebenschuh, Carlo↗
DOE OSTI2023

Nanoparticle Surface Passivation Strategies to Improve the Cycle and Calendar Lifetime of Silicon Anodes in Standard Carbonate Electrolytes

Silicon Lithium alloys (SiLix) as the anode active material in a Li-Ion battery configuration offer possible energy densities paralleled only by pure lithium metal. However, the extreme mechanical deformation of alloying and dealloying SiLix through charge/discharge cycles paired with the highly reactive interface of SiLix which generates an interfacial layer known as the solid-electrolyte-interphase are large barriers to industrial adoption of high-silicon-content negative electrodes. Moreover, these challenges are magnified when the thickness of the electrode is brought to relevant levels (>3 mg/cm2). Here, I describe our efforts to address these challenges by utilizing single-nanometer-scale silicon nanoparticles to reduce capacity fade related to mechanical failure. I also detail our efforts to modify the silicon surface through interfacial chemical engineering strategies to passivate the silicon surface. This development three different majority silicon electrodes (50-74 wt%) to achieve cycle capacity retention of greater than 73% through 1000 charge/discharge cycles against capacity-matched NMC-based cathodes and calendar lifetimes greater than one year. This research is a part of the multi-national lab Silicon Consortium Project.

anodes↗
DOE OSTI2022

A massively parallel time-domain coupled electrodynamics–micromagnetics solver

We present a high-performance coupled electrodynamics–micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain approach for Maxwell’s equations to a magnetization model described by the Landau–Lifshitz–Gilbert equation. The algorithm is implemented in the Exascale Computing Project software framework, AMReX, which provides effective scalability on manycore and GPU-based supercomputing architectures. Furthermore, the code leverages ongoing developments of the Exascale Application Code, WarpX, which is primarily being developed for plasma wakefield accelerator modeling. Our temporal coupling scheme provides second-order accuracy in space and time by combining the integration steps for the magnetic field and magnetization into an iterative sub-step that includes a trapezoidal temporal discretization for the magnetization. The performance of the algorithm is demonstrated by the excellent scaling results on NERSC multicore and GPU systems, with a significant (59×) speedup on the GPU using a node-by-node comparison. We demonstrate the utility of our code by performing simulations of an electromagnetic waveguide and a magnetically tunable filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Parallel deterministic transport sweeps of structured and unstructured meshes with overloaded mesh decompositions

Here, the efficiency of discrete ordinates transport sweeps depends on the scheduling algorithm, the domain decomposition, the problem to be solved, and the computational platform. Sweep scheduling algorithms may be categorized by their approach to several issues. In this paper we examine the strategy of domain overloading for mesh partitioning as one of the components of such algorithms. In particular, we extend the domain overloading strategy, previously defined and analyzed for structured meshes, to the general case of unstructured meshes. We also present computational results for both the structured and unstructured domain overloading cases. We find that an appropriate amount of domain overloading can greatly improve the efficiency of parallel sweeps for both structured and unstructured partitionings of the test problems examined on up to 10 5 processor cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Partitioning and Communication Strategies for Sparse Non-negative Matrix Factorization

Non-negative matrix factorization (NMF), the problem of finding two non-negative low-rank factors whose product approximates an input matrix, is a useful tool for many data mining and scientific applications such as topic modeling in text mining and unmixing in microscopy. In this paper, we focus on scaling algorithms for NMF to very large sparse datasets and massively parallel machines by employing effective algorithms, communication patterns, and partitioning schemes that leverage the sparsity of the input matrix. We consider two previous works developed for related problems, one that uses a fine-grained partitioning strategy using a point-to-point communication pattern and one that uses a Cartesian, or checkerboard, partitioning strategy using a collective-based communication pattern. We show that a combination of the previous approaches balances the demands of the various computations within NMF algorithms and achieves high efficiency and scalability. From the experiments, we see that our proposed strategy runs up to 10x faster than the state of the art on real-world datasets.

Kaya, Oguz↗
DOE OSTI2020

Efficient solvers for hybridized three-field mixed finite element coupled poromechanics

We consider a mixed hybrid finite element formulation for coupled poromechanics. A stabilization strategy based on a macro-element approach is advanced to eliminate the spurious pressure modes appearing in undrained/incompressible conditions. The efficient solution of the stabilized mixed hybrid block system is addressed by developing a class of block triangular preconditioners based on a Schur-complement approximation strategy. Robustness, computational efficiency and scalability of the proposed approach are theoretically discussed and tested using challenging benchmark problems on massively parallel architectures.

42 ENGINEERING↗
DOE OSTI2017

Comprehensive End-to-End Design of Novel High Energy Density Materials: I. Synthesis and Characterization of Oxadiazole Based Heterocycles

A methodology to design novel energetic materials by means of a holistic approach that links synthesis, experimental characterization, quantum-chemical modeling, and statistical empirical evaluation is proposed. An analysis of the revealed structure–property–function correlations in the LLM compound series (oxadiazole-based heterocyclic energetics), BNFF, BNFF-1, LLM-172, LLM-191, and LLM-192, led us to predict, obtain, and characterize a new member in the materials family, LLM-200, which exhibits attractive energetic characteristics compared to known conventional high energy density materials. Finally, while the applied strategy convincingly demonstrated feasibility of the end-to-end design of high energy density materials, there are certain limitations in parallel improvements of sensitivity and performance within a single compound.

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

Supramolecular Assemblies on Surfaces: Nanopatterning, Functionality, and Reactivity

Understanding how molecules interact to form large-scale hierarchical structures on surfaces holds promise for building designer nanoscale constructs with defined chemical and physical properties. Here, we describe early advances in this field and highlight upcoming opportunities and challenges. Both direct intermolecular interactions and those that are mediated by coordinated metal centers or substrates are discussed. These interactions can be additive, but they can also interfere with each other, leading to new assemblies in which electrical potentials vary at distances much larger than those of typical chemical interactions. Earlier spectroscopic and surface measurements have provided partial information on such interfacial effects. In the interim, scanning probe microscopies have assumed defining roles in the field of molecular organization on surfaces, delivering deeper understanding of interactions, structures, and local potentials. Furthermore, self-assembly is a key strategy to form extended structures on surfaces, advancing nanolithography into the chemical dimension and providing simultaneous control at multiple scales. In parallel, the emergence of graphene and the resulting impetus to explore 2D materials have broadened the field, as surface-confined reactions of molecular building blocks provide access to such materials as 2D polymers and graphene nanoribbons. In this Review, we describe recent advances and point out promising directions that will lead to even greater and more robust capabilities to exploit designer surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗