SEARCH · Search NASA

Results for “Time parallelization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Nonlinear damping of a finite amplitude whistler wave due to modified two stream instability

A two-dimensional, fully kinetic, particle-in-cell simulation is used to investigate the nonlinear development of a parallel propagating finite amplitude whistler wave (parent wave) with a wavelength longer than an ion inertial length. The cross field current of the parent wave generates short-scale whistler waves propagating highly oblique directions to the ambient magnetic field through the modified two-stream instability (MTSI) which scatters electrons and ions parallel and perpendicular to the magnetic field, respectively. The parent wave is largely damped during a time comparable to the wave period. The MTSI-driven damping process is proposed as a cause of nonlinear dissipation of kinetic turbulence in the solar wind.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗

DOE OSTI2023

A simple, flexible technique for RF cavity wake-field calculations

It is typical in the accelerator field to model machine components, especially RF cavities, as parallel RLC resonators. In the interest of simulating and diagnosing beam instabilities, knowledge of the time-domain voltage waveform over an equivalent resonator by a bunch current often proves useful. This waveform may be found by convolving the bunch current with the RLC resonator impulse response. While analytical and quasi-analytical expressions are available in this regime for common distributions such as the Gaussian, analogous results for less standard distributions are difficult to obtain using direct methods, which opens the door for the development of a more generalized technique. In this paper, a formulation is created that allows for the simple computation of the time-domain voltage waveform of and RLC resonator. The formulation uses the Cauchy Residue Theorem to extract the convolution result from the Fourier Domain, and it only requires that the current distribution Fourier Transform be holomorphic and known at one specific evaluation point. This greatly simplifies the computation of the time domain voltage for a large amount of bunch distributions both common and uncommon.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Flow of the normal component of He-II about bluff objects observed with He$^*_2$ excimers

Flow of the normal fluid component of He II about a cylinder and flat plate in a channel with a square cross-section was observed by tracking clouds of He$^*_2$ excimers. Flow was produced using a lithographically patterned heater. The direction and speed of flow projected onto a plane parallel to the axis of the channel and perpendicular to the heater exhibited significant change with time even after the heater was turned off. Flow was recorded in movies. Excimers moved at least as fast as could be recorded by the camera. Velocity vector field maps of flow suggest formation of structures downstream of the cylinder object that are consistent with parts of eddies. Here, in this paper, we establish a foundation to observe normal component flow of a quantum fluid over centimeters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Floquet-driven indirect exchange interaction mediated by topological insulator surface states

Light drives offer a potential tool for the dynamical control of magnetic interactions in matter. We theoretically investigate the indirect exchange coupling between two parallel chains of magnetic impurities on the surface of a topological insulator, driven by a time-periodic circularly polarized light field in the high-frequency, off-resonant regime. We derive a closed-form analytic expression for the spin susceptibility of the photon-dressed topological insulator surface states and obtain the irradiation dependence of the Ising, Heisenberg, and Dzyaloshinskii-Moriya exchange couplings between the impurity chains. Our results show a two-pronged modification of these exchange couplings by periodic drives. First, the Ruderman-Kittel-Kasuya-Yosida (RKKY) oscillation period of the exchange couplings can be extended by enhancing the driving strength. Second, increasing driving strength enhances the envelope of RKKY oscillations of the Ising type while suppressing those of the Heisenberg type and Dzyaloshinskii-Moriya type. Furthermore, our work provides useful insights for realizing Floquet engineering of collinear and noncollinear indirect exchange interactions in topological insulating systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2017

Concurrent conversation modeling and parallel simulation of the naming game in social networks

The Naming Game is an effective self-organization model to understand the emergence of linguistic consensus and to investigate the system dynamics in a variety of phenomena over social networks of autonomous agents. The Naming Game is an effective description for the evolution of consensus despite the absence of any central coordination or specialized initialization even in large-scale networks. While the classical game is effective in description, it was defined with inherently sequential evaluation semantics over the entire network. Here, we develop a new concurrent model as a relaxation of the classical formulation and express it in a discrete event style of evaluation. Further, with the uncovered concurrency that was absent in the classical algorithm, we map the concurrent model to parallel discrete event simulation. Using a prototype implementation, we present an initial parallel performance study on networks containing hundreds of thousands of individuals, with a decrease in simulation time in the best-observed case from 4800 seconds down to 1400 seconds.

Perumalla, Kalyan↗

DOE OSTI2019

Avoiding Communication in Primal and Dual Block Coordinate Descent Methods

Primal and dual block coordinate descent methods are iterative methods for solving regularized and unregularized optimization problems. Distributed-memory parallel implementations of these methods have become popular in analyzing large machine learning datasets. However, existing implementations communicate at every iteration, which, on modern data center and supercomputing architectures, often dominates the cost of floating-point computation. Recent results on communication-avoiding Krylov subspace methods suggest that large speedups are possible by re-organizing iterative algorithms to avoid communication. We show how applying similar algorithmic transformations can lead to primal and dual block coordinate descent methods that only communicate every $s$ iterations---where $s$ is a tuning parameter---instead of every iteration for the regularized least-squares problem. We show that the communication-avoiding variants reduce the number of synchronizations by a factor of $s$ on distributed-memory parallel machines without altering the convergence rate and attain strong scaling speedups of up to $6.1\times$ over the “standard algorithm" on a Cray XC30 supercomputer.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Dynamical regimes of vortex flow in type-II superconductors with parallel twin boundaries

We explore the dynamics of driven magnetic flux lines in disordered type-II superconductors in the presence of twin boundaries oriented parallel to the direction of the applied magnetic field, using a three-dimensional elastic line model simulated with Langevin molecular dynamics. The lines are driven perpendicular to the planes to model the effect of an electric current applied parallel to the planes and perpendicular to the magnetic field. A study of the long-time non-equilibrium steady states for varying sample orientation and thickness reveals a rich collection of dynamical regimes spanning the depinning crossover region that separates the pinned and moving-lattice states of vortex matter. We observe the emergence of a preferred direction for the ordering of the Abrikosov lattice in the free- owing vortex regime due to asymmetric pinning by the planar defects. Furthermore, we have performed novel direct measurements of flux line excitations such as half-loops and double kinks to aid the characterization of the topologically rich flux flow profile.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Electron Coherence and Interactions in Nanostructures (Final Report)

The goal of this program was to develop the tools of quantum many-body theory suitable for studying electronic properties of nanostructures, to apply them to recent or contemporaneous experiments, and to chart the directions of future development in the physics of low-dimensional quantum materials and devices. Over the duration of this Program, we introduced and developed the notion of Nonlinear Luttinger liquid, which is viewed by peers as one of a few major developments in the theory of one-dimensional quantum liquids since the notion of Tomonaga-Luttinger liquids was introduced by Tomonaga in 1950. Later, motivated by the raise of the ideas of topological quantum computing, we contributed to the theory of structure and dynamics of Majorana states in superconducting devices. In parallel, we worked to understand the sources of decoherence of experimentally-implemented superconducting qubits. Over time, we developed a comprehensive theory of the qubit relaxation by quasiparticles. This theory is widely used in the interpretation of experiments and in developing practical ways to mitigate the qubit decoherence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

VARIANCE ANISOTROPY IN KINETIC PLASMAS

Solar wind fluctuations admit well-documented anisotropies of the variance matrix, or polarization, related to the mean magnetic field direction. Typically, one finds a ratio of perpendicular variance to parallel variance of the order of 9:1 for the magnetic field. Here we study the question of whether a kinetic plasma spontaneously generates and sustains parallel variances when initiated with only perpendicular variance. We find that parallel variance grows and saturates at about 5% of the perpendicular variance in a few nonlinear times irrespective of the Reynolds number. For sufficiently large systems (Reynolds numbers) the variance approaches values consistent with the solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Material fabrication using acoustic radiation forces

Apparatus and methods for using acoustic radiation forces to order particles suspended in a host liquid are described. The particles may range in size from nanometers to millimeters, and may have any shape. The suspension is placed in an acoustic resonator cavity, and acoustical energy is supplied thereto using acoustic transducers. The resulting pattern may be fixed by using a solidifiable host liquid, forming thereby a solid material. Patterns may be quickly generated; typical times ranging from a few seconds to a few minutes. In a one-dimensional arrangement, parallel layers of particles are formed. With two and three dimensional transducer arrangements, more complex particle configurations are possible since different standing-wave patterns may be generated in the resonator. Fabrication of periodic structures, such as metamaterials, having periods tunable by varying the frequency of the acoustic waves, on surfaces or in bulk volume using acoustic radiation forces, provides great flexibility in the creation of new materials. Periodicities may range from millimeters to sub-micron distances, covering a large portion of the range for optical and acoustical metamaterials.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2018

Material fabrication using acoustic radiation forces

36 MATERIALS SCIENCE↗

DOE OSTI2022

Novel Solver Algorithms for Nearly Singular Linear Systems Arising in Combustion Modelling

Direct Numerical Simulations of realistic combustion devices are extremely challenging due to the wide separation of scales in the simulation, for example an internal combustion (IC) engine chamber, and the flame thickness of a high-pressure flame. The PeleLMeX solver uses adaptive mesh refinement (AMR) to evolve multi-species reacting flows in the low Mach number limit at the Exascale and relies on an embedded boundary (EB) approach to represent complex geometries. In that framework, the EB geometries often give rise to very small cut-cells along the boundary, which translate into extreme ill-conditioning of the pressure-projection, with eigenvalues that span 15-16 orders of magnitude. In this talk, we focus on the case of a typical IC piston bowl geometry for which we present on a novel approach towards solving these nearly singular linear systems with ILU-based, C-AMG smoothers on massively parallel architectures. In particular, we use scaling and equilibration algorithms to handle the non-normality of the upper triangular factors. This enables us to approximate the highly sequential triangular solve algorithm, embedded in the AMG smoothing-solve phase, with Jacobi iterations. This approximation can be written as a convergent Neumann series whose terms are composed of highly parallel sparse matrix vector multiplications. The result is an algorithm that substantially decreases setup and solve time, compared to state-of-the-art, for these challenging linear systems.

combustion modelling↗

DOE OSTI2018

Short-Term Load Forecasting Based Automatic Distribution Network Reconfiguration

In a traditional dynamic network reconfiguration study, the optimal topology is determined at every scheduled time point by using the real load data measured at that time. The development of the load forecasting technique can provide an accurate prediction of the load power that will happen in a future time and provide more information about load changes. With the inclusion of load forecasting, the optimal topology can be determined based on the predicted load conditions during a longer time period instead of using a snapshot of the load at the time when the reconfiguration happens; thus, the distribution system operator can use this information to better operate the system reconfiguration and achieve optimal solutions. This paper proposes a short-term load forecasting approach to automatically reconfigure distribution systems in a dynamic and pre-event manner. Specifically, a short-term and high-resolution distribution system load forecasting approach is proposed with a forecaster based on support vector regression and parallel parameters optimization. The network reconfiguration problem is solved by using the forecasted load continuously to determine the optimal network topology with the minimum amount of loss at the future time. The simulation results validate and evaluate the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Packaging a 650V/400A GaN Half-bridge Power Module with Ultra-low Parasitics for Electric Vehicle Drive Applications

This paper proposes a compact and efficient half-bridge power module with three 650 V / 150 A GaN dies in parallel. The power module incorporates a main power printed circuit board (PCB), an interface PCB, and a flex PCB to achieve low parasitics in both power loop and gate-side connection, resolving the issue of high parasitics typically encountered with wire bonding in high-current applications. Additionally, the interface PCB decouples the design constraints between the power loop and the gate loops. The proposed design is optimized with a vertical loop configuration to reduce power loop inductance through magnetic flux cancellation. Finite element analysis indicates that the power loop inductance is 0.58 nH at 100 MHz, while the maximum die junction temperature reaches 131 °C under an ambient temperature of 65 °C and a load current of 385 A. The proposed multi-piece PCB structure reduces the inductance of the drive circuit to minimize EMI and to mitigate false triggering. At the same time, it reduces impedance mismatches across different driver circuits, thereby achieving dynamic current sharing in multi-chip parallel configurations. Under simulation conditions of 400 V / 385 A, the current imbalance among chips was limited to 5 A. A 400 V / 385 A double-pulse test was conducted to experimentally validate the performance of the proposed power module.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2020

Parallel algorithms for hyperdynamics and local hyperdynamics

Hyperdynamics (HD) is a method for accelerating the timescale of standard molecular dynamics (MD). It can be used for simulations of systems with an energy potential landscape that is a collection of basins, separated by barriers, where transitions between basins are infrequent. HD enables the system to escape from a basin more quickly while enabling a statistically accurate renormalization of the simulation time, thus effectively boosting the timescale of the simulation. In [Kim, Perez, Voter, J Chem Phys, 139:144110, 2013)1, a local version of HD was formulated, which exploits the intrinsic locality characteristic typical of most systems to mitigate the poor scaling properties of standard HD as the system size is increased. In this paper, we discuss how both HD and local HD can be formulated to run efficiently in parallel. We have implemented these ideas in the LAMMPS MD code, which means HD can be used with any interatomic potential LAMMPS supports. Together, these parallel methods allow simulations of any size to achieve the time acceleration offered by HD (which can be orders of magnitude), at a cost 3-5x that of standard MD. As examples, we performed two simulations of a million-atom system to model the diffusion and clustering of Pt adatoms on a large patch of Pt(100) surface for 80 and 160 μs.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗