Search NASA⌕ Search

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.

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

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Carburization and decarburization behavior of Grade 91 ferritic-martensitic steel in liquid sodium environments

Here, this paper presents a study of carbon transfer and its effect on microstructure and tensile properties of Grade 91 (G91) ferritic-martensitic steel exposed to sodium at 550–650 °C. Sodium exposure tests were conducted in Argonne's forced convection sodium loops up to exposure times of ~40,000 h. Thermal aging study of G91 steel was conducted in parallel to isolate the thermal aging effect from the sodium effect. It was found that sodium exposures at 650 °C dissolved M 23 C 6 carbides, eliminated the martensite subgrain structure resulting in excessive grain growth and reduced the tensile strength by >50%, while sodium exposures at 550 and 600 °C had an insignificant effect on its microstructure and tensile properties. These effects were attributed to the carburization/decarburization process of G91 steel in sodium environments. Carbon concentrations in sodium were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8–1.2 ppm in the SMT-1 loop and 0.3–0.7 ppm in the SMT-2 loop. Thermodynamic analysis of the carburization – decarburization process was conducted for G91 steels exposed in sodium environments. The carbon activity-concentration relationship for G91 was evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC and VC carbides dominated the process at low carbon activities. The calculated carburization-decarburization boundary showed that G91 would undergo decarburization at 650 °C and carburization at 550 °C in the sodium loop environments, which was consistent with our experimental observations. This experimental and theoretical analysis provided a basis for predicting the effect of carbon transfer on the integrity of reactor components in sodium environments and for the design of new alloys used in sodium-cooled fast reactors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer

A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green’s functions in real space and imaginary time. Furthermore, these integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol –1 for the C 60 dimer at the center distance of 9.748 Å.

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

Whistler wave generation by anisotropic tail electrons during asymmetric magnetic reconnection in space and laboratory

Whistler wave generation near the magnetospheric separatrix during reconnection at the dayside magnetopause is studied with data from the Magnetospheric Multiscale (MMS) mission. The dispersion relation of the whistler mode is measured for the first time near the reconnection region in space, which shows that whistler waves propagate nearly parallel to the magnetic field line. A linear analysis indicates that the whistler waves are generated by temperature anisotropy in the electron tail population. This is caused by loss of electrons with a high velocity parallel to the magnetic field to the exhaust region. Here, there is a positive correlation between activities of whistler waves and the lower–hybrid drift instability (LHDI) both in laboratory and space, indicating the enhanced transport by LHDI may be responsible for the loss of electrons with a high parallel velocity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.

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

The flush-mounted rail Langmuir probe array designed for the Alcator C-Mod vertical target plate divertor

An array of flush-mounted and toroidally elongated Langmuir probes (henceforth called rail probes) have been specifically designed for the Alcator C-Mod’s vertical target plate divertor and operated over multiple campaigns. The “flush” geometry enables the tungsten electrodes to survive high heat flux conditions in which traditional “proud” tungsten electrodes suffer damage from melting. The toroidally elongated rail-like geometry reduces the influence of sheath expansion, which is an important effect to consider in the design and interpretation of flush-mounted Langmuir probes. The new rail probes successfully operated during C-Mod’s FY2015 and FY2016 experimental campaigns with no evidence of damage, despite being regularly subjected to heat flux densities parallel to the magnetic field exceeding ~1 GW m -2 for short periods of time. A comparison between rail and proud probe data indicates that sheath expansion effects were successfully mitigated by the rail design, extending the use of these Langmuir probes to incident magnetic field line angles as low as 0.5°.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Triplet-pair spin signatures from macroscopically aligned heteroacenes in an oriented single crystal

The photo-driven process of singlet fission generates coupled triplet pairs (TT) with fundamentally intriguing and potentially useful properties. The quintet 5 TT 0 sublevel is particularly interesting for quantum information because it is highly entangled, is addressable with microwave pulses, and could be detected using optical techniques. Previous theoretical work on a model Hamiltonian and nonadiabatic transition theory, called the JDE model, has determined that this sublevel can be selectively populated if certain conditions are met. Among the most challenging, the molecules within the dimer undergoing singlet fission must have their principal magnetic axes parallel to one another and to an applied Zeeman field. Here, we present time-resolved electron paramagnetic resonance (TR-EPR) spectroscopy of a single crystal sample of a tetracenethiophene compound featuring arrays of dimers aligned in this manner, which were mounted so that the orientation of the field relative to the molecular axes could be controlled. The observed spin sublevel populations in the paired TT and unpaired (T+T) triplets are consistent with predictions from the JDE model, including preferential 5 TT 0 formation at z ‖ B 0 , with one caveat—two 5 TT spin sublevels have little to no population. This may be due to crossings between the 5 TT and 3 TT manifolds in the field range investigated by TR-EPR, consistent with the intertriplet exchange energy determined by monitoring photoluminescence at varying magnetic fields.

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

SAIGE-GPU: accelerating genome- and phenome-wide association studies using GPUs

Genome-wide association studies (GWAS) at biobank scale are computationally intensive, especially for admixed populations requiring robust statistical models. SAIGE is a widely used method for generalized linear mixed-model GWAS but is limited by its CPU-based implementation, making phenome-wide association studies impractical for many research groups. We developed SAIGE-GPU, a GPU-accelerated version of SAIGE that replaces CPU-intensive matrix operations with GPU-optimized kernels. The core innovation is distributing genetic relationship matrix calculations across GPUs and communication layers. Applied to 2068 phenotypes from 635 969 participants in the Million Veteran Program, including diverse and admixed populations, SAIGE-GPU achieved a 5-fold speedup in mixed model fitting on supercomputing infrastructure and cloud platforms. We further optimized the variant association testing step through multi-core and multi-trait parallelization. Deployed on Google Cloud Platform and Azure, the method provided substantial cost and time savings. Source code and binaries are available for download at https://github.com/saigegit/SAIGE/tree/SAIGE-GPU-1.3.3. A code snapshot is archived at Zenodo for reproducibility (DOI: [10.5281/zenodo.17642591]). SAIGE-GPU is available in a containerized format for use across HPC and cloud environments and is implemented in R/C++ and runs on Linux systems.

Rodriguez, Alex [Argonne National Laboratory (ANL)↗
DOE OSTI2025

Entanglement-enhanced ac magnetometry in the presence of Markovian noise

Entanglement is a resource to improve the sensitivity of quantum sensors. In an ideal case, using an entangled state as a probe to detect target fields, we can beat the standard quantum limit by which all classical sensors are bounded. However, since entanglement is fragile against decoherence, it is unclear whether entanglement-enhanced metrology is useful in a noisy environment. Its benefit is indeed limited when estimating the amplitude of dc magnetic fields under the effect of parallel Markovian decoherence, where the noise operator is parallel to the target field. In this paper, on the contrary, we show an advantage to using an entanglement over the classical strategy under the effect of parallel Markovian decoherence when we try to detect ac magnetic fields. We consider a scenario to induce a Rabi oscillation of the qubits with the target ac magnetic fields. Although we can, in principle, estimate the amplitude of the ac magnetic fields from the Rabi oscillation, the signal becomes weak if the qubit frequency is significantly detuned from the frequency of the ac magnetic field. We show that, by using the Greenberger-Horne-Zeilinger (GHZ) states, we can significantly enhance the signal of the detuned Rabi oscillation even under the effect of parallel Markovian decoherence. Further, our method is based on the fact that the interaction time between the GHZ states and ac magnetic fields scales as 1/L to mitigate the decoherence effect, where L is the number of qubits, which contributes to improving the bandwidth of the detectable frequencies of the ac magnetic fields. Our results pave the way for new applications of entanglement-enhanced ac magnetometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Whistler wave generation by anisotropic tail electrons during asymmetric magnetic reconnection in space and laboratory

Whistler wave generation near the magnetospheric separatrix during reconnection at the dayside magnetopause is studied with data from the Magnetospheric Multiscale (MMS) mission. The dispersion relation of the whistler mode is measured for the first time near the reconnection region in space, which shows that whistler waves propagate nearly parallel to the magnetic field line. A linear analysis indicates that the whistler waves are generated by temperature anisotropy in the electron tail population. This is caused by loss of electrons with a high velocity parallel to the magnetic field to the exhaust region. There is a positive correlation between activities of whistler waves and the lower-hybrid drift instability (LHDI) both in laboratory and space, indicating the enhanced transport by LHDI may be responsible for the loss of electrons with a high parallel velocity.

Dispersion relation↗
DOE OSTI2021

HYPPO v1.0.0

HYPPO is a hyperparameter optimization (HPO) software that uses surrogate modeling and uncertainty quantification to provide reliable neural network architectures. Similar state-of-the-art HPO technologies do not use such features and mostly perform random searches to explore the parameter space, which can easily become computationally expensive. Using HYPPO, the user is able to reduce by an order of magnitude the number of evaluations necessary to identify the most optimal region in the hyperparameter space. Finally, using asynchronous nested parallelism, we are able to decrease by two orders of magnitude the amount of time required to complete the HPO process.

Dumont, Vincent↗
DOE OSTI2021

Multiplexer System for the SPEAR3 Booster BPM Upgrade

BPM measurements in booster synchrotrons are often only critical during accelerator commissioning or when a problem occurs. As a result, many facilities do not make large investments in booster BPM signal processors; they either have very few BPMs and/or use older generation processors. The SPEAR3 booster BPM processor system, for instance, has operated since 1990 with commercial multiplexers to switch between BPM button signals into a single dated analog BPM processor that was developed at SLAC*. This system has reached its end-of-life so we are in the process of upgrading to modern multiplexers that feed a pair of turn-by-turn Libera SPARK-ERXR processors. This low-cost solution gives us the ability to arbitrarily multiplex between BPM signals during the energy ramp with modern BPM processors. The system can either measure 2 BPMs turn-by-turn in parallel during the entire energy ramp, or sequentially measure all BPMs (2 at a time) at different time slices within the ramp. Here we show measurements of the MiniCircuits switch we chose as well as our architecture for the upgrade.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Slycat™ User Manual

Slycat™ is a web-based system for performing data analysis and visualization of potentially large quantities of remote, high-dimensional data. Slycat™ specializes in working with ensemble data. An ensemble is a group of related data sets, which typically consists of a set of simulation runs exploring the same problem space. An ensemble can be thought of as a set of samples within a multi-variate domain, where each sample is a vector whose value defines a point in high-dimensional space. To understand and describe the underlying problem being modeled in the simulations, ensemble analysis looks for shared behaviors and common features across the group of runs. Additionally, ensemble analysis tries to quantify differences found in any members that deviate from the rest of the group. The Slycat™ system integrates data management, scalable analysis, and visualization. Results are viewed remotely on a user’s desktop via commodity web clients using a multi-tiered hierarchy of computation and data storage, as shown in Figure 1. Our goal is to operate on data as close to the source as possible, thereby reducing time and storage costs associated with data movement. Consequently, we are working to develop parallel analysis capabilities that operate on High Performance Computing (HPC) platforms, to explore approaches for reducing data size, and to implement strategies for staging computation across the Slycat™ hierarchy. Within Slycat™, data and visual analysis are organized around projects, which are shared by a project team. Project members are explicitly added, each with a designated set of permissions. Although users sign-in to access Slycat™, individual accounts are not maintained. Instead, authentication is used to determine project access. Within projects, Slycat™ models capture analysis results and enable data exploration through various visual representations. Although for scientists each simulation run is a model of real-world phenomena given certain conditions, we use the term model to refer to our modeling of the ensemble data, not the physics. Different model types often provide complementary perspectives on data features when analyzing the same data set. Each model visualizes data at several levels of abstraction, allowing the user to range from viewing the ensemble holistically to accessing numeric parameter values for a single run. Bookmarks provide a mechanism for sharing results, enabling interesting model states to be labeled and saved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

SUNDIALS Multiphysics+MPIManyVector Performance Testing

In this report we document performance test results on a SUNDIALS-based multiphysics demonstration application. We aim to assess the large-scale parallel performance of new capabilities that have been added to the SUNDIALS suite of time integrators and nonlinear solvers in recent years under funding from both the Exascale Computing Project (ECP) and the Scientific Discovery through Advanced Scientific (SciDAC) program, specifically: 1. SUNDIALS’ new MPIManyVector module, that allows extreme flexibility in how a solution "vector" is staged on computational resources. 2. ARKode’s new multirate integration module, MRIStep, allowing high-order accurate calculations that subcycle "fast" processes within "slow" ones. 3. SUNDIALS’ new flexible linear solver interfaces, that allow streamlined specification of problem-specific linear solvers. 4. SUNDIALS’ new N_Vector additions of "fused" vector operations (to increase arithmetic intensity) and separation of reduction operations into "local" and "global" versions (to reduce latency by combining multiple reductions into a single MPI_Allreduce call). We anticipate that subsequent reports will extend this work to investigate a variety of other new features, including SUNDIALS’ generic SUNNonlinearSolver interface and accelerator-enabled N Vector modules, and upcoming MRIStep extensions to support custom "fast" integrators (that leverage problem structure) and IMEX integration of the "slow" time scale (to add diffusion).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Xyce™ Parallel Electronic Simulator (V.7.9) Release Notes

The Xyce™ Parallel Electronic Simulator has been written to support the simulation needs of Sandia National Laboratories’ electrical designers. Xyce™ is a SPICE-compatible simulator with the ability to solve extremely large circuit problems on large-scale parallel computing platforms, but also includes support for most popular parallel and serial computers. For up-to-date information not available at the time these notes were produced, please visit the Xyce™ web page at http://xyce.sandia.gov.

42 ENGINEERING↗
DOE OSTI2025

Angle-resolved polarized Raman spectroscopy study of phosphorene nanoribbons

We present a systematic angle-resolved polarized Raman spectroscopy (ARPRS) study of black phosphorus (BP) nanostructures formed via electrochemical sodium- and lithium-ion intercalation. Sodium intercalation leads to bundles of densely packed, highly uniform phosphorene nanoribbons (PNRs) separated by parallel amorphous channels, whereas lithium intercalation results in shorter, irregular nanoribbon-like segments with lower aspect ratios. In both cases, six additional Raman peaks (P1–P6) appear alongside the three primary Raman-active modes of BP (A$^{1}_{g}$, B 2g , and A$^{2}_{g}$). These peaks are attributed to the amorphous regions, as confirmed by their isotropic angular dependence in ARPRS measurements. The three BP modes show pronounced angular variations that differ significantly between the two intercalated samples. In sodium-intercalated BP, A$^{1}_{g}$ and A$^{2}_{g}$ modes retain a dumbbell-like angular dependence under parallel polarization with enhanced anisotropy and reduced symmetry under crossed polarization. At the same time, the B 2g mode transitions from four-lobed (cloverleaf) polar plot to a butterfly-like one. In contrast, lithium-intercalated BP exhibits weaker anisotropy and less distinct angular polar plots for all three modes. These differences reflect the sensitivity of phonon behavior to underlying nanostructure morphology. The vibrational frequencies density of states (FDOS) calculations attribute the B 2g mode transformation to phonon band folding and mode mixing in PNRs. This study demonstrates the power of ARPRS in probing phonon-structure relationships and highlights the influence of edge geometry and quantum confinement on phonon dispersion in PNRs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Sub-Alfvénic Turbulence: Magnetic-to-kinetic Energy Ratio, Modification of Weak Cascade, and Implications for Magnetic Field Strength Measurements

Abstract We study the properties of sub-Alfvénic magnetohydrodynamic (MHD) turbulence, i.e., turbulence with Alfvén mach number M A = V L / V A < 1, where V L is the velocity at the injection scale and V A is the Alfvén velocity. We demonstrate that MHD turbulence can have different properties, depending on whether it is driven by velocity or magnetic fluctuations. If the turbulence is driven by isotropic bulk forces acting upon the fluid, i.e., is velocity driven, in an incompressible conducting fluid we predict that the kinetic energy is M A − 2 times larger than the energy of magnetic fluctuations. This effect arises from the long parallel wavelength tail of the forcing, which excites modes with k ∥ / k ⊥ < M A . We also predict that as the MHD turbulent cascade reaches the strong regime, the energy of slow modes exceeds the energy of Alfvén modes by a factor M A − 1 . These effects are absent if the turbulence is driven through magnetic fluctuations at the injection scale. We confirm our predictions with numerical simulations. Since the assumption of magnetic and kinetic energy equipartition is at the core of the Davis–Chandrasekhar–Fermi (DCF) approach to measuring magnetic field strength in sub-Alfvénic turbulence, we conclude that the DCF technique is not universally applicable. In particular, we suggest that the dynamical excitation of long azimuthal wavelength modes in the galactic disk may compromise the use of the DCF technique. We discuss alternative expressions that can be used to obtain magnetic field strength from observations and consider ways of distinguishing the cases of velocity and magnetically driven turbulence using observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Simulations of Magnetic Fields in Tidally Disrupted Stars

We perform the first magnetohydrodynamical simulations of tidal disruptions of stars by supermassive black holes. We consider stars with both tangled and ordered magnetic fields, for both grazing and deeply disruptive encounters. When the star survives disruption, we find its magnetic field amplifies by a factor of up to 20, but see no evidence for a self-sustaining dynamo that would yield arbitrary field growth. For stars that do not survive, and within the tidal debris streams produced in partial disruptions, we find that the component of the magnetic field parallel to the direction of stretching along the debris stream only decreases slightly with time, eventually resulting in a stream where the magnetic pressure is in equipartition with the gas. Our results suggest that the returning gas in most (if not all) stellar tidal disruptions is already highly magnetized by the time it returns to the black hole.

79 ASTRONOMY AND ASTROPHYSICS↗