Search NASA⌕ Search

SEARCH · Search NASA

Results for “parallel methods”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Massively parallel phase-field simulations targeting exascale

The interface thickness in the phase-field (PF) method limits its simulation scales. Consequently, large-scale PF simulations become prohibitively expensive for resolving the extremely fine microstructures that typically form during rapid solidification processing. This challenge is significant in predicting microstructure evolution in metal additive manufacturing and has been identified by the United States Department of Energy’s Exascale Computing Project. Here, to address this, we develop a multi-GPU and MPI-based massively parallel simulation code, utilizing state-of-the-art algorithms, software, and libraries, for large-scale three-dimensional (3D) PF simulations. We report the first GPU-parallel PF simulations on Frontier (currently the second TOP500 exascale cluster) and Summit machines, taking dendritic growth as an example problem. We evaluate the parallel performance of our implementation using scaling studies with more than 24 000 GPUs (among the largest known computations to date) and the acceleration performance using large-scale simulations of dendritic growth in 3D. Finally, massively parallel GPUs in these supercomputers enabled the first coupled multiscale simulations of laser melting and subsequent dendritic solidification on the scale of a full melt-pool, demonstrating the feasibility of performing PF simulations with a point total over 2 billion grid points within an acceptable time.

Exascale↗
DOE OSTI2023

Multirate partitioned Runge–Kutta methods for coupled Navier–Stokes equations

Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. Further, this study explores multirate partitioned Runge-Kutta methods for the fluid-fluid interaction problem and demonstrates its parallel performance by using the PETSc library. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Impact of pellet preparation on reflectance measurements of ammonium sulfate by single-angle reflectance spectroscopy

The optical constants, n and k, are required to model the reflection, refraction and transmission of light at the first surface interface of a material, as well as its propagation through the material; here n corresponds to the real component and k the imaginary component of the refractive index. Several different spectroscopic methods have been used to determine the n and k values for different materials, including single angle reflectance spectroscopy and ellipsometry. The single-angle reflectance method quantitatively records the specular reflectance R(v) from a plane parallel face of the material and uses the Kramers-Kronig transform to extract the n and k values. For most compounds, however, it is difficult to obtain a single crystal or high-quality window of sufficient planarity and of the appropriate dimensions (several mm) to make the measurement. For this reason, we further investigate the use of pressed pellets of neat powdered substances to measure optical constants of these substances using the single-angle reflectance method. We have found that surface roughness can significantly influence the measured quantitative reflectance spectrum R(v) and, consequently, the derived n and k parameters. A collaborative study between Defence Research and Development Canada (DRDC) - Valcartier Research Center and Pacific Northwest National Laboratory (PNNL) has been carried out using different pellets of neat ammonium sulfate [(NH4)2SO4] to show how parameters such as the relative fraction of particle sizes for the component powder along with the grinding process of the pure compound all can affect the reflectance spectrum used to derive the optical constants. All pressed pellets were characterized by single-angle reflectance spectroscopy.

Diaz, Emmanuela↗
DOE OSTI2022

Sensor assemblies and methods for emulating interaction of entities within water systems

Sensor assemblies are provided for use in modeling water systems. These sensor assemblies can be used as sensor fish. These assemblies can include a circuit board supporting processing circuitry components on either or both opposing component support surfaces of the circuit board and a housing above the circuit board and the components, with the housing being circular about the circuit board in at least one cross section, and wherein the supporting surfaces of the circuit board are substantially parallel with the plane of the housing in the one cross section. Methods for emulating interaction of entities within water systems are provided. The methods can include introducing a sensor assembly into a water system. The sensor assembly can include: a circuit board supporting processing circuitry components on either or both of opposing component support surfaces of the circuit board; a housing about the circuit board and the components, the housing being circular about the circuit board in at least one cross section; and wherein the support surfaces of the circuit board are substantially parallel with the plane of the housing in the one cross section.

Deng, Z. Daniel↗
NASA NTRS2015

An Efficient Solution Method for Multibody Systems with Loops Using Multiple Processors

This paper describes a multibody dynamics algorithm formulated for parallel implementation on multiprocessor computing platforms using the divide-and-conquer approach. The system of interest is a general topology of rigid and elastic articulated bodies with or without loops. The algorithm divides the multibody system into a number of smaller sets of bodies in chain or tree structures, called "branches" at convenient joints called "connection points", and uses an Order-N (O (N)) approach to formulate the dynamics of each branch in terms of the unknown spatial connection forces. The equations of motion for the branches, leaving the connection forces as unknowns, are implemented in separate processors in parallel for computational efficiency, and the equations for all the unknown connection forces are synthesized and solved in one or several processors. The performances of two implementations of this divide-and-conquer algorithm in multiple processors are compared with an existing method implemented on a single processor.

Multibody dynamics↗
DOE OSTI2015

Nitrogen-modified nano-titania: True phase composition, microstructure and visible-light induced photocatalytic NO{sub x} abatement

Titanium dioxide (TiO{sub 2}) is a popular photocatalyst used for many environmental and anti-pollution applications, but it normally operates under UV light, exploiting ∼5% of the solar spectrum. Nitrification of titania to form N-doped TiO{sub 2} has been explored as a way to increase its photocatalytic activity under visible light, and anionic doping is a promising method to enable TiO{sub 2} to harvest visible-light by changing its photo-absorption properties. In this paper, we explore the insertion of nitrogen into the TiO{sub 2} lattice using our green sol–gel nanosynthesis method, used to create 10 nm TiO{sub 2} NPs. Two parallel routes were studied to produce nitrogen-modified TiO{sub 2} nanoparticles (NPs), using HNO{sub 3}+NH{sub 3} (acid-precipitated base-peptised) and NH{sub 4}OH (totally base catalysed) as nitrogen sources. These NPs were thermally treated between 450 and 800 °C. Their true phase composition (crystalline and amorphous phases), as well as their micro-/nanostructure (crystalline domain shape, size and size distribution, edge and screw dislocation density) was fully characterised through advanced X-ray methods (Rietveld-reference intensity ratio, RIR, and whole powder pattern modelling, WPPM). As pollutants, nitrogen oxides (NO{sub x}) are of particular concern for human health, so the photocatalytic activity of the NPs was assessed by monitoring NO{sub x} abatement, using both solar and white-light (indoor artificial lighting), simulating outdoor and indoor environments, respectively. Results showed that the onset of the anatase-to-rutile phase transformation (ART) occurred at temperatures above 450 °C, and NPs heated to 450 °C possessed excellent photocatalytic activity (PCA) under visible white-light (indoor artificial lighting), with a PCA double than that of the standard P25 TiO{sub 2} NPs. However, higher thermal treatment temperatures were found to be detrimental for visible-light photocatalytic activity, due to the effects of four simultaneous occurrences: (i) loss of OH groups and water adsorbed on the photocatalyst surface; (ii) growth of crystalline domain sizes with decrease in specific surface area; (iii) onset and progress of the ART; (iv) the increasing instability of the nitrogen in the titania lattice. - Graphical abstract: Nitrogen modified TiO{sub 2} synthesised via a green aqueous sol–gel method showed to degrade nitrogen oxides (NO{sub x}) under visible white-light (indoor artificial lighting), with a photocatalytic activity double than that of the standard P25 TiO{sub 2} NPs. - Highlights: • N–TiO{sub 2} synthesised via a green aqueous sol–gel method. • Advanced X-ray methods used to detect both crystalline and amorphous contents. • Microstructure fully addressed via XRPD and whole powder pattern modelling. • Photocatalytic NO{sub x} removal assessed using both solar and visible-light lamps.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Boundary layer receptivity and control

Receptivity processes initiate natural instabilities in a boundary layer. The instabilities grow and eventually break down to turbulence. Consequently, receptivity questions are a critical element of the analysis of the transition process. Success in modeling the physics of receptivity processes thus has a direct bearing on technological issues of drag reduction. The means by which transitional flows can be controlled is also a major concern: questions of control are tied inevitably to those of receptivity. Adjoint systems provide a highly effective mathematical method for approaching many of the questions associated with both receptivity and control. The long term objective is to develop adjoint methods to handle increasingly complex receptivity questions, and to find systematic procedures for deducing effective control strategies. The most elementary receptivity problem is that in which a parallel boundary layer is forced by time-harmonic sources of various types. The characteristics of the response to such forcing form the building blocks for more complex receptivity mechanisms. The first objective of this year's research effort was to investigate how a parallel Blasius boundary layer responds to general direct forcing. Acoustic disturbances in the freestream can be scattered by flow non-uniformities to produce Tollmien-Schlichting waves. For example, scattering by surface roughness is known to provide an efficient receptivity path. The present effort is directed towards finding a solution by a simple adjoint analysis, because adjoint methods can be extended to more complex problems. In practice, flows are non-parallel and often three-dimensional. Compressibility may also be significant in some cases. Recent developments in the use of Parabolized Stability Equations (PSE) offer a promising possibility. By formulating and solving a set of adjoint parabolized equations, a method for mapping the efficiency with which external forcing excites the three-dimensional motions of a non-parallel boundary layer was developed. The method makes use of the same computationally efficient formulation that makes the PSE currently so appealing. In the area of flow control, adjoint systems offer a powerful insight into the effect of control forces. One of the simplest control strategies for boundary layers involves the application of localized mean wall suction.

Hill, D. C.↗
NASA NTRS1982

A parallel VLSI architecture for a digital filter of arbitrary length using Fermat number transforms

A parallel architecture for computation of the linear convolution of two sequences of arbitrary lengths using the Fermat number transform (FNT) is described. In particular a pipeline structure is designed to compute a 128-point FNT. In this FNT, only additions and bit rotations are required. A standard barrel shifter circuit is modified so that it performs the required bit rotation operation. The overlap-save method is generalized for the FNT to compute a linear convolution of arbitrary length. A parallel architecture is developed to realize this type of overlap-save method using one FNT and several inverse FNTs of 128 points. The generalized overlap save method alleviates the usual dynamic range limitation in FNTs of long transform lengths. Its architecture is regular, simple, and expandable, and therefore naturally suitable for VLSI implementation.

Truong, T. K.↗
DOE OSTI2023

Virtual Time III, Part 1: Unified Virtual Time Synchronization for Parallel Discrete Event Simulation

Algorithms for synchronization of parallel discrete event simulation have historically been divided between conservative methods that require lookahead but not rollback, and optimistic methods that require rollback but not lookahead. In this paper we present a new approach in the form of a framework called Unified Virtual Time (UVT) that unifies the two approaches, combining the advantages of both within a single synchronization theory. Whenever timely lookahead information is available, a logical process (LP) executes conservatively using an irreversible event handler. When lookahead information is not available the LP does not block, as it would in a classical conservative execution, but instead executes optimistically using a reversible event handler. The switch from conservative to optimistic synchronization and back is decided on an event-by-event basis by the simulator, transparently to the model code. UVT treats conservative synchronization algorithms as optional accelerators for an underlying optimistic synchronization algorithm, enabling the speed of conservative execution whenever it is applicable, but otherwise falling back on the generality of optimistic execution. We describe UVT in a novel way, based on fundamental invariants, monotonicity requirements, and synchronization rules. UVT permits zero-delay messages and pays careful attention to tie-handling using superposition. We prove that under fairly general conditions a UVT simulation always makes progress in virtual time. This is Part 1 of a trio of papers describing the UVT framework for PDES, mixing conservative and optimistic synchronization and integrating throttling control.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2003

FLAPS (Fatigue Life Analysis Programs): Computer Programs to Predict Cyclic Life Using the Total Strain Version of Strainrange Partitioning and Other Life Prediction Methods

This manual presents computer programs FLAPS for characterizing and predicting fatigue and creep-fatigue resistance of metallic materials in the high-temperature, long-life regime for isothermal and nonisothermal fatigue. The programs use the Total Strain version of Strainrange Partitioning (TS-SRP), and several other life prediction methods described in this manual. The user should be thoroughly familiar with the TS-SRP and these life prediction methods before attempting to use any of these programs. Improper understanding can lead to incorrect use of the method and erroneous life predictions. An extensive database has also been developed in a parallel effort. The database is probably the largest source of high-temperature, creep-fatigue test data available in the public domain and can be used with other life-prediction methods as well. This users' manual, software, and database are all in the public domain and can be obtained by contacting the author. The Compact Disk (CD) accompanying this manual contains an executable file for the FLAPS program, two datasets required for the example problems in the manual, and the creep-fatigue data in a format compatible with these programs.

Arya, Vinod K.↗
DOE OSTI2019

Recovering noise-free quantum observables

We introduce a technique for recovering noise-free observables in noisy quantum systems by combining the results of many slightly different experiments. Our approach is applicable to a variety of quantum systems, but we illustrate it with applications to quantum information and quantum sensing. The approach corresponds to repeating the same quantum evolution many times with known variations on the underlying systems' error properties, e.g., the spontaneous emission and dephasing times $T_1$ and $T_2$. As opposed to standard quantum error correction methods, which have an overhead in the number of qubits (many physical qubits must be added for each logical qubit), our method has only an overhead in the number of evaluations, allowing the overhead to, in principle, be hidden via parallelization. We show that the effective spontaneous emission $T_1$ and dephasing $T_2$ times can be increased using this method in both simulation and experiments on an actual quantum computer. We also show how to correct more complicated entangled states and how Ramsey fringes relevant to quantum sensing can be significantly extended in time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Decentralized Interleaving of Paralleled Dc-Dc Buck Converters: Preprint

We present a decentralized control strategy that yields switch interleaving among parallel connected dc-dc buck converters without communication. The proposed method is based on the digital implementation of the dynamics of a nonlinear oscillator circuit as the controller. Each controller is fully decentralized, i.e., it only requires the locally measured output current to synthesize the pulse width modulation (PWM) carrier waveform. By virtue of the intrinsic electrical coupling between converters, the nonlinear oscillator-based controllers converge to an interleaved state with uniform phase-spacing across PWM carriers. To the knowledge of the authors, this work represents the first fully decentralized strategy for switch interleaving of paralleled dc-dc buck converters.

14 SOLAR ENERGY↗
NASA NTRS1988

A Strassen-Newton algorithm for high-speed parallelizable matrix inversion

Techniques are described for computing matrix inverses by algorithms that are highly suited to massively parallel computation. The techniques are based on an algorithm suggested by Strassen (1969). Variations of this scheme use matrix Newton iterations and other methods to improve the numerical stability while at the same time preserving a very high level of parallelism. One-processor Cray-2 implementations of these schemes range from one that is up to 55 percent faster than a conventional library routine to one that is slower than a library routine but achieves excellent numerical stability. The problem of computing the solution to a single set of linear equations is discussed, and it is shown that this problem can also be solved efficiently using these techniques.

Bailey, David H.↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Understanding performance variability in standard and pipelined parallel Krylov solvers

In this work, we collect data from runs of Krylov subspace methods and pipelined Krylov algorithms in an effort to understand and model the impact of machine noise and other sources of variability on performance. We find large variability of Krylov iterations between compute nodes for standard methods that is reduced in pipelined algorithms, directly supporting conjecture, as well as large variation between statistical distributions of runtimes across iterations. Based on these results, we improve upon a previously introduced nondeterministic performance model by allowing iterations to fluctuate over time. We present our data from runs of various Krylov algorithms across multiple platforms as well as our updated non-stationary model that provides good agreement with observations. We also suggest how it can be used as a predictive tool.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

Mixing of supersonic jets including the effects of supersonic jets including the effects of transverse pressure gradient using difference methods

The usual boundary-layer equations describing the steady-state mixing of parallel jets are supplemented by the momentum equation in the direction normal to the flow. This allows detailed computation of the flow field in the mixing region and simultaneous computation of the outer inviscid flow. An explicit and an implicit finite difference scheme have been developed and applied in several illustrative examples. The examples include mixing of planar and axisymmetric supersonic jets of different composition with both matched and unmatched static pressures. Numerical results were compared with available experimental data obtained for the unmatched pressure case.

Kurkov, A. P.↗
NASA NTRS1988

On multigrid methods for the Navier-Stokes Computer

The overall architecture of the multipurpose parallel-processing Navier-Stokes Computer (NSC) being developed by Princeton and NASA Langley (Nosenchuck et al., 1986) is described and illustrated with extensive diagrams, and the NSC implementation of an elementary multigrid algorithm for simulating isotropic turbulence (based on solution of the incompressible time-dependent Navier-Stokes equations with constant viscosity) is characterized in detail. The present NSC design concept calls for 64 nodes, each with the performance of a class VI supercomputer, linked together by a fiber-optic hypercube network and joined to a front-end computer by a global bus. In this configuration, the NSC would have a storage capacity of over 32 Gword and a peak speed of over 40 Gflops. The multigrid Navier-Stokes code discussed would give sustained operation rates of about 25 Gflops.

Nosenchuck, D. M.↗