Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel systems”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Flexible User-Defined Domain Decomposition in Kilometer-Scale E3SM Land Model Simulation

The Energy Exascale Earth System Model (E3SM) Land Model (ELM) has been extended to kilometer-scale (km-ELM) resolutions, enabling high-fidelity simulations of terrestrial processes at 1 km x 1 km grid spacing. In ELM, domain decomposition partitions the computational domain across processors, ensuring efficient parallel execution. Currently, round-robin decomposition is applied, providing a straightforward way to distribute computational workload. As ELM continues evolving at the kilometer-scale (km-scale), particularly with integrating lateral flow modeling, decomposition strategies must also account for the increased workload and data movement. This paper introduces a flexible user-defined domain decomposition framework, allowing users to customize domain partitioning based on application requirements. The impact of different decomposition strategies is evaluated across various applications concerning computation, communication, and I/O. Results demonstrate that while 1D partitioning yields superior I/O performance, k-nearest neighbors (KNN) clustering effectively reduces inter-process communication overhead. This study lays the groundwork for scalable partitioning in large-scale land surface simulations, enhancing next-generation Earth system modeling.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗
DOE OSTI2022

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

LW3D was used to explore the expected CSR from an undulator-based diagnostic using a first-principles approach. It was also used for comparison with an analytic Green-function approach. LW3D was extended to model transport systems involving dipoles and quadrupoles. It was used to explore the radiation produced in a 4-dipole bunch compressor. Radiation ahead of the electron bunch at the beamline exit was shown to be associated with wavefronts produced in upstream dipoles overtaking the bunch.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

ECP Software Technology Capability Assessment Report

The Exascale Computing Project Software Technology (ECP ST) focus area represents the key bridge between Exascale systems and the scientists developing applications that will run on those platforms. ECP offers a unique opportunity to build a coherent set of software (often referred to as the "software stack") that will allow application developers to maximize their ability to write highly parallel applications, targeting multiple Exascale architectures with runtime environments that will provide high performance and resilience. But applications are only useful if they can provide scientific insight, and the unprecedented data produced by these applications require a complete analysis work ow that includes new technology to scalably collect, reduce, organize, curate, and analyze the data into actionable decisions. This requires approaching scientific computing in a holistic manner, encompassing the entire user workflow - from conception of a problem, setting up the problem with validated inputs, performing high-fidelity simulations, to the application of uncertainty quantification to the final analysis. The software stack plan defined here aims to address all of these needs by extending current technologies to Exascale where possible, by performing the research required to conceive of new approaches necessary to address unique problems where current approaches will not suffice, and by deploying high-quality and robust software products on the platforms developed in the Exascale systems project. The ECP ST portfolio has established a set of interdependent projects that will allow for the research, development, and delivery of a comprehensive software stack,

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toward a Fully Integrated Multiphysics Simulation Framework for Fusion Blanket Design

Fusion is an attractive clean-energy solution, thanks to its various advantages, such as reduced radioactivity, little high-level nuclear waste, ample fuel supplies, and increased safety. However, the harsh operating environment introduced by a complex fusion plasma system makes design and integration of fusion blankets incredibly challenging and time-consuming. This work focuses on developing a fully integrated multiphysics simulation framework based on an advanced open-source platform—the Multiphysics Object-Oriented Simulation Environment (MOOSE)—to alleviate the difficulties in fusion blanket design and integration. MOOSE is a massively parallel finite element/volume multiphysics simulation platform that has been widely adopted within the nuclear fission community. Even though fission and fusion are fundamentally different, they involve similar multiphysics phenomena. A fully integrated open-source multiphysics simulation framework tailored for the fusion blanket design will be implemented by leveraging the well-established multiphysics capabilities in MOOSE. Once successfully developed, this fully integrated framework will rapidly evaluate a blanket design concept and offer insights for subsequent iterations. As the first step, we will mainly aim to integrate neutronics analysis, system thermal hydraulics simulation, and full 3-D heat transfer calculations. The efficacy of the integrated framework will be verified using an innovative solid ceramic blanket design. While the project’s final goal is to enable a fully integrated multiphysics simulation platform for various fusion blanket concepts, here this work, as a preliminary step, will mainly focus on a solid ceramic breeder helium-cooled blanket.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimal interconnection of three-terminal tandem solar cells

Three-terminal (3T) tandem solar cells require an adapted module integration scheme in order to explore their full efficiency potential. The three terminals allow to extract the power of the top and bottom cell separately. In a cell string, the wide bandgap top cells are interconnected in parallel to multiple bottom cells resulting in a parallel/series interconnection. This interconnection scheme affects the operation of the subcells, the resulting current path between the subcells, the layout of the cell interconnects, and the system level. Here, we analyze by simulations and experiments the aspects of the module integration of series- and reverse-connected 3T cells with their practical impact on module processes and performance as well as the effect of varying voltage ratios on the string-end losses. If the subcells are connected in series, the module integration requires insulation layers and significantly longer interconnects compared to devices with reverse-connected subcells. Tandem devices with a reverse connection and a voltage ratio between top and bottom cell of 2:1 allow a lean interconnection design and low integration losses. Here, we present an approach for the integration of bypass diodes for the protection against shading effects that allow to minimize string-end and shading losses for a system of modules featuring 3T cells.

14 SOLAR ENERGY↗
DOE OSTI2022

Multiscale Crystalline Structure of Confined Polypeptoid Films: The Effect of Alkyl Side Chain Branching

Here we report the effect of alkyl side chain branching on melt-recrystallization of nanoconfined polypeptoid films using poly(N-octyl glycine) (PNOG) and poly(N-2-ethyl-1-hexyl glycine) (PNEHG) as model systems. Upon cooling from the isotropic melt, confined PNOG molecules recrystallize into a near-perfect orthorhombic crystal structure with the board-like molecules stacked face-to-face in the substrate-parallel direction, resulting in long-range ordered wormlike lamellae that occupy the entire film. By contrast, rod-like PNEHG molecules bearing branched N-2-ethyl-1-hexyl side chains stack into a columnar hexagonal mesophase with their backbones oriented parallel to the substrates, forming micron-sized sheaf-like superstructures under confinement, exposing large areas of empty spaces in the film. These findings highlight the effect of alkyl side chain branching on the packing motif and multiscale crystalline structure of polypeptoids under a nanoconfined geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Intrinsic Ion Transport Properties of Block Copolymer Electrolytes

Knowledge of intrinsic properties is of central importance for materials design and assessing suitability for specific applications. Self-assembling block copolymer electrolytes (BCEs) are of great interest for applications in solid-state energy storage devices. A fundamental understanding of ion transport properties, however, is hindered by the difficulty in deconvoluting extrinsic factors, such as defects, from intrinsic factors, such as the presence of interfaces between the domains. Here, we quantify the intrinsic ion transport properties of a model BCE system consisting of poly(styrene-block-ethylene oxide) (SEO) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) salt using a generalizable strategy of depositing thin films on interdigitated electrodes and self-assembling fully connected parallel lamellar structures throughout the films. Comparison between conductivity in homopolymer poly(ethylene oxide) (PEO)-LiTFSI electrolytes and the analogous conducting material in SEO over a range of salt concentrations (r, molar ratio of lithium ion to ethylene oxide repeat units) and temperatures reveals that between 20% and 50% of the PEO in SEO is inactive. Using mean-field theory calculations of the domain structure and monomer concentration profiles at domain interfaces-both of which vary substantially with salt concentration-the fraction of inactive PEO in the SEO, as derived from conductivity measurements, can be quantitatively reconciled with the fraction of PEO that is mixed with greater than a few volume percent of polystyrene. Despite the detrimental interfacial effects for ion transport in BCEs, the intrinsic conductivity of the SEO studied here (ca. 10 -3 S/cm at 90 degrees C, r = 0.085) is an order of magnitude higher than reported values from bulk samples of similar molecular weight SEO (ca. 10 -4 S/cm at 90 degrees C, r = 0.085). Finally, this work provides motivation and methods for pursuing improved BCE chemical design, interfacial engineering, and processing.

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

A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios

Abstract We formulate a generalized periodic boundary condition as a limit of the standard twist-and-shift parallel boundary condition that is suitable for simulations of plasmas with low magnetic shear. This is done by applying a phase shift in the binormal direction when crossing the parallel boundary. While this phase shift can be set to zero without loss of generality in the local flux-tube limit when employing the twist-and-shift boundary condition, we show that this is not the most general case when employing periodic parallel boundaries, and may not even be the most desirable. A non-zero phase shift can be used to avoid the convective cells that plague simulations of the three-dimensional Hasegawa–Wakatani system, and is shown to have measurable effects in periodic low-magnetic-shear gyrokinetic simulations. We propose a numerical program where a sampling of periodic simulations at random pseudo-irrational flux surfaces are used to determine physical observables in a statistical sense. This approach can serve as an alternative to applying the twist-and-shift boundary condition to low-magnetic-shear scenarios, which, while more straightforward, can be computationally demanding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Substantial reduction of the anisotropy in the critical current densities J c of Ni-doped CaKFe 4 As 4 single crystals by chemical and irradiation-induced disorder

Here, we demonstrate that the anisotropy in the critical current densities, J c , of CaK(Fe 1–x Ni x ) 4 As 4 (x = 0, 0.015, 0.025, and 0.030) single crystals reduces with increasing chemical and irradiation-induced disorder. The magnetic field dependences of J c are analyzed by performing magnetization measurements with H applied parallel and perpendicular to the crystallographic c-axis. The results show that undoped crystals display large anisotropies in Jc due to an enhancement of the vortex pinning with H applied parallel to the crystallographic ab-planes. This anisotropy reduces substantially as Ni addition increases. Moreover, we found that random disorder introduced by proton irradiation enhances mainly the vortex pinning for H parallel to the c-axis. Consequently, using adequate fluencies, the vortex pinning at low temperatures in both undoped and doped samples becomes close to isotropic. These results make the CaKFe 4 As 4 system promising for applications that require isotropic J c under magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Status and future plans for C 3 R&D

C 3 is an opportunity to realize an e + e - collider for the study of the Higgs boson at √s = 250 GeV, with a well defined upgrade path to 550 GeV while staying on the same short facility footprint. C 3 is based on a fundamentally new approach to normal conducting linear accelerators that achieves both high gradient and high efficiency at relatively low cost. Given the advanced state of linear collider designs, the key system that requires technical maturation for C 3 is the main linac. This paper presents the staged approach towards a facility to demonstrate C 3 technology with both Direct (source and main linac) and Parallel (beam delivery, damping ring, ancillary component) R&D. The primary goal of the C 3 Demonstration R&D Plan is to reduce technical and cost risk by building and operating the key components of C 3 at an adequate scale. This R&D plan starts with the engineering design, and demonstration of one cryomodule and will culminate in the construction of a 3 cryomodule linac with pre-production prototypes. This R&D program would also demonstrate the linac rf fundamentals including achievable gradient and gradient stability over a full electron bunch train and breakdown rates. It will also investigate beam dynamics including energy spread, wakefields, and emittance growth. This work will be critical to confirm the suitability of the C 3 beam parameters for the physics reach and detector performance in preparation for a Conceptual Design Report (CDR), as well as for follow-on technology development and industrialization. The C 3 Demonstration R&D Plan will open up significant new scientific and technical opportunities based on development of high-gradient and high-efficiency accelerator technology. It will push this technology to operate both at the GeV scale and mature the technology to be reliable and provide high-brightness electron beams. The timeline for progressing with C 3 technology development will be governed by practical limitations on both the technical progress and resource availability. It consists of four stages: Stage 0) Ongoing fundamental R&D on structure prototypes, damping and vibrations. Stage 1) Advancing the engineering maturity of the design and developing start-to-end simulations including space-charge and wakefield effects. This stage will include testing of strucutres operating at cryogenic temperatures. Beam tests would be performed with high beam current to test full beam loading. Stage 2) Production and testing of the first cryomodule at cryogenic temperatures. This would provide sufficient experimental data to compile a CDR and it is anticipated for Stage 2 to last 3 years and to culminate with the transport of photo-electrons through the first cryomodule. Stage 3) Updates to the engineering design of the cryomodules, production of the second and third cryomodule and their installation. Lower charge and lower emittance beams will be used to investigate emittance growth. The successful full demonstration of the 3 cryomodules to deliver up to a 3 GeV beam and achieve the C 3 five gradient will allow a comprehensive and robust evaluation of the technical design of C 3 as well as mitigate technical, schedule, and cost risks required to proceed with a Technical Design Report (TDR).

radiation hardened magnets↗
DOE OSTI2021

Grand Unification of Quantum Algorithms

Quantum algorithms offer significant speed-ups over their classical counterparts for a variety of problems. The strongest arguments for this advantage are borne by algorithms for quantum search, quantum phase estimation, and Hamiltonian simulation, which appear as subroutines for large families of composite quantum algorithms. A number of these quantum algorithms have recently been tied together by a novel technique known as the quantum singular value transformation (QSVT), which enables one to perform a polynomial transformation of the singular values of a linear operator embedded in a unitary matrix. In the seminal GSLW’19 paper on the QSVT [Gilyén et al., ACM STOC 2019], many algorithms are encompassed, including amplitude amplification, methods for the quantum linear systems problem, and quantum simulation. Here, we provide a pedagogical tutorial through these developments, first illustrating how quantum signal processing may be generalized to the quantum eigenvalue transform, from which the QSVT naturally emerges. Paralleling GSLW’19, we then employ the QSVT to construct intuitive quantum algorithms for search, phase estimation, and Hamiltonian simulation, and also showcase algorithms for the eigenvalue threshold problem and matrix inversion. This overview illustrates how the QSVT is a single framework comprising the three major quantum algorithms, suggesting a grand unification of quantum algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗
DOE OSTI2021

Error-controlled, progressive, and adaptable retrieval of scientific data with multilevel decomposition

Extreme-scale simulations and high-resolution instruments have been generating an increasing amount of data, which poses significant challenges to not only data storage during the run, but also post-processing where data will be repeatedly retrieved and analyzed for a long period of time. The challenges in satisfying a wide range of post-hoc analysis needs while minimizing the I/O overhead caused by inappropriate and/or excessive data retrieval should never be left unmanaged. In this paper, we propose a data refactoring, compressing, and retrieval framework capable of 1) fine-grained data refactoring with regard to precision; 2) incrementally retrieving and recomposing the data in terms of various error bounds; and 3) adaptively retrieving data in multi-precision and multi-resolution with respect to different analysis. With the progressive data re-composition and the adaptable retrieval algorithms, our framework significantly reduces the amount of data retrieved when multiple incremental precision are requested and/or the downstream analysis time when coarse resolution is used. Experiments show that the amount of data retrieved under the same progressively requested error bound using our framework is 64% less than that using state-of-the-art single-error-bounded approaches. Parallel experiments with up to 1, 024 cores and ~ 600 GB data in total show that our approach yields 1.36× and 2.52× performance over existing approaches in writing to and reading from persistent storage systems, respectively.

Liang, Xin↗
DOE OSTI2023

Laser-Induced Ultrafast Demagnetization and Spin Switching (Final Report)

The information technology revolution demands bigger and faster magnetic storage devices. The funded research focuses on a new frontier - Laser-induced ultrafast all-optical spin switching (AOS), where a single laser pulse can switch spins from one direction to another permanently, without a magnetic field, on a time scale of 1-10 picoseconds, 2-3 orders of magnitude faster than the existing technology. The funded research has successfully resolved a long-standing issue in ultrafast demagnetization whether the lattice contributes to the demagnetization. We developed a high level first-principles theory and employed a massively parallel computer to pin down the details of the physical mechanisms that are inaccessible to other methods. This unique feature allows us to compute the force on each atom in both bulk and monolayer systems. We find that it is indeed true there is a nonzero force on each atom, but the idea of vibration contribution is not crucial to our basic understanding of physics. This paves the way to future experimental, theoretical and technological developments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Chemoselective, osmium-free, dihydroxylation/oxidative cleavage of heteroaryl isoprenes by a contemporary Malaprade reaction

The methyl ketone is a central synthetic building block for the construction of advanced heteroaryl scaffolds and systems. Reactions, including oxidative cyclization strategies, are often predicated on efficient access to this ubiquitous moiety. In the context of arenes, standard approaches leveraging Markovnikov hydration/oxidation or oxidative cleavage of the C–C π bond often afford satisfactory performance. However, when the substrate contains an electron-deficient heteroaryl core, the traditional Malaprade reaction, and related oxidative-cleavage strategies, frequently result in diminished performance over carbon-based arenes. In this work we present the development and application of an oxidative cleavage reaction of various pyridinyl isoprenes towards accessing the downstream methyl ketone for utilization in advanced cyclizations for the preparation of soft-N-donor complexant scaffolds. Further, this efficient protocol parallels the principles of Green chemistry by exchanging KMnO 4 for the toxic OsO 4 and offers the end-user an efficient, more environmentally friendly option for accessing heteroaryl methyl ketones in one hour of reaction time using potassium permanganate and sodium paraperiodate as a synergistically potent oxidative cleavage system. The wide substrate scope defined access to simple, as well as advanced heteroaryl methyl ketones. Method development, optimization, substrate scope, preliminary mechanistic observations, and a scale up reaction are delineated herein.

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

Memory access optimization for particle operations in computational fluid dynamics-discrete element method simulations

Computational Fluid Dynamics - Discrete Element Method is used to model gas-solid systems in several applications in energy, pharmaceutical and petrochemical industries. Computational performance bottlenecks often limit the problem sizes that can be simulated at industrial scale. The data structures used to store several millions of particles in such large-scale simulations have a large memory footprint that does not fit into the processor cache hierarchies on current high-performance-computing platforms, leading to reduced computational performance. This paper specifically addresses this aspect of memory access bottlenecks in industrial scale simulations. The use of space-filling curves to improve memory access patterns is described and their impact on computational performance is quantified in both shared and distributed memory parallelization paradigms. The Morton space filling curve applied to uniform grids and k-dimensional tree partitions are used to reorder the particle data-structure thus improving spatial and temporal locality in memory. The performance impact of these techniques when applied to two benchmark problems, namely the homogeneous-cooling-system and a fluidized-bed, are presented. We report these optimization techniques lead to approximately two-fold performance improvement in particle focused operations such as neighbor-list creation and data-exchange, with ~ 1.5 times overall improvement in a fluidization simulation with 1.27 million particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Influence of Supporting Electrolyte on Hydroxide Exchange Membrane Water Electrolysis Performance: Catholyte

Alkaline or hydroxide exchange membrane water electrolysis (HEMWE) is a promising technology for green hydrogen production using platinum group metal-free catalysts and stainless steel, an advantage of alkaline water electrolysis (AWE), and a gas-impermeable membrane, a parallel to proton exchange membrane electrolysis (PEMWE). However, the HEMWE requires supporting electrolytes and there is minimal understanding of their role on the respective reactions. Without SELs, HEMWE performance and durability are worse than PEMWE systems. Herein, consistently feeding potassium hydroxide anolyte, we systematically study the effects of catholyte SELs in HEMWEs including dry vs. wet operation, cation effects, anion effects, and cation/OH ratios on cell potential and stability. We report that (i) hydration of the cathode improves high current density operation by preventing dehydration of the hydroxide exchange membrane (HEM), (ii) there was no correlation between cation type and cell potential, (iii) cell potential and high frequency resistance did not correlate with SEL conductivity, (iv) cathodic carbonate SEL had a significant negative effect on cell performance, (v) increased cation/OH ratio also caused increased cell potentials. Overall, this study concludes that feeding water or potassium hydroxide solution is desirable to improve the AEMWE performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗