Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Drift kinetic electrostatic simulations of the edge localized mode heat pulse

In the present work, electrostatic drift kinetic simulations of parallel plasma transport within the tokamak scrape-off layer (SOL) are conducted using the COGENT code. The SOL configuration is represented in one-dimensional slab geometry, incorporating a heat source localized in the midplane. The heat source parameters correspond to those characterizing edge-localized modes observed in the Joint European Torus (JET) tokamak. The numerical model includes kinetic treatment of both ions and electrons, a simplified model for the gyrokinetic Poisson equation that allows one to step over short time scales associated with fast electrostatic shear Alfvèn waves, and the logical sheath boundary condition (LSBC) that enforces global system quasineutrality. A third-order accurate LSBC is derived to be consistent with the third-order accurate upwind advection scheme utilized in the code, and it was shown to noticeably impact the simulation results, especially parallel heat flux at the target plate. The findings of this study are in agreement with results from preceding fluid and kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling performance of data collection systems for high-energy physics

Exponential increases in scientific experimental data are outpacing silicon technology progress, necessitating heterogeneous computing systems—particularly those utilizing machine learning (ML)—to meet future scientific computing demands. The growing importance and complexity of heterogeneous computing systems require systematic modeling to understand and predict the effective roles for ML. We present a model that addresses this need by framing the key aspects of data collection pipelines and constraints and combining them with the important vectors of technology that shape alternatives, computing metrics that allow complex alternatives to be compared. For instance, a data collection pipeline may be characterized by parameters such as sensor sampling rates and the overall relevancy of retrieved samples. Alternatives to this pipeline are enabled by development vectors including ML, parallelization, advancing CMOS, and neuromorphic computing. By calculating metrics for each alternative such as overall F1 score, power, hardware cost, and energy expended per relevant sample, our model allows alternative data collection systems to be rigorously compared. We apply this model to the Compact Muon Solenoid experiment and its planned high luminosity-large hadron collider upgrade, evaluating novel technologies for the data acquisition system (DAQ), including ML-based filtering and parallelized software. The results demonstrate that improvements to early DAQ stages significantly reduce resources required later, with a power reduction of 60% and increased relevant data retrieval per unit power (from 0.065 to 0.31 samples/kJ). However, we predict that further advances will be required in order to meet overall power and cost constraints for the DAQ.

Olin-Ammentorp, Wilkie (ORCID:0000000224729862)↗
DOE OSTI2025

Off-state magnetoresistance in long-channel germanium Schottky-barrier MOSFETs

An increasing magnetic field perpendicular to an undoped semiconductor surface at low temperature is known to strengthen the binding of localized electrons to stationary ions, as the wavefunction's tails evolve from exponential to Gaussian. It is also known that application of a high bias voltage to a depleted semiconductor can liberate bound charge and induce a large drop in electrical resistance. We connect these established results to experimental electrical transport measurements on off-state germanium Schottky-barrier metal–oxide–semiconductor field-effect transistor (MOSFETs) with an aluminum oxide insulating dielectric and platinum germanide contacts. We make measurements at the three distinct orientations of the magnetic field with respect to the substrate and the current. At 6 K, we observe sharp attenuation of current by more than 2 orders of magnitude, within 60 mT, at a crossover magnetic field perpendicular to the substrate. A 1 T magnetic field attenuates the current by more than 4 orders of magnitude. The strength of the attenuation and the value of the crossover field are controlled by both the gate–source and drain–source voltages. The attenuation is much weaker when the magnetic field is parallel to the current. Finally, we orient the magnetic field parallel to the substrate, but perpendicular to the current, allowing us to distinguish charge hopping at the oxide interface from charge hopping in the bulk. In conclusion, this large off-state magnetoresistance can be exploited for cryogenic magnetic- and photo-detection, and for high-bias, low-leakage MOSFETs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quadrupolar density structures in driven magnetic reconnection experiments with a guide field

Magnetic reconnection is a ubiquitous process in plasma physics, driving rapid and energetic events such as coronal mass ejections. Reconnection between magnetic fields with arbitrary shear can be decomposed into an anti-parallel reconnecting component and a non-reconnecting guide-field component, which is parallel to the reconnecting electric field. This guide field modifies the structure of the reconnection layer and the reconnection rate. We present results from experiments on the MAIZE pulsed-power generator (500 kA peak current, 200 ns rise time), which use two exploding wire arrays, tilted in opposite directions, to embed a guide field in the plasma flows with a relative strength b≡B g /B rec =0, 0.4, or 1. The reconnection layers in these experiments have widths that are less than the ion skin depth, d i =c/ω pi , indicating the importance of the Hall term, which generates a distinctive quadrupolar magnetic field structure along the separatrices of the reconnection layer. Using laser imaging interferometry, we observe quadrupolar structures in the line-integrated electron density, consistent with the interaction of the embedded guide field with the quadrupolar Hall field. Our measurements extend over much larger length scales (40d i ) at higher β (∼1) than previous experiments, providing an insight into the global structure of the reconnection layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On transport of electrons in a magnetized plasma sheath

Electron trajectories are analyzed in the plasma sheath region in a magnetic field that intersects the surface at a shallow angle. Depending on the size of the Debye width relative to the electron gyroradius, λ D /ρ e , an electron either drifts in the sheath electric field or follows a sequence of reflections (“bouncing”) on the surface before reversing its parallel velocity. For thin sheath, λ D /ρ e ≪ 1, the conservation of the electron magnetic moment can be violated, depending on the ratio of the electron perpendicular and parallel velocities V ⊥ /V || and the magnetic field pitch angle with respect to the surface. Small-scale waviness of the surface can strongly amplify this effect, leading to violation of the magnetic moment conservation even for λ D /ρ e ≳ 1, which opens a channel of anomalous radial transport for electrons in plasma devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗
DOE OSTI2026

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗
DOE OSTI2023

Structural insights into plasmalemma vesicle-associated protein (PLVAP): Implications for vascular endothelial diaphragms and fenestrae

In many organs, small openings across capillary endothelial cells (ECs) allow the diffusion of low–molecular weight compounds and small proteins between the blood and tissue spaces. These openings contain a diaphragm composed of radially arranged fibers, and current evidence suggests that a single-span type II transmembrane protein, plasmalemma vesicle-associated protein-1 (PLVAP), constitutes these fibers. Here, we present the three-dimensional crystal structure of an 89-amino acid segment of the PLVAP extracellular domain (ECD) and show that it adopts a parallel dimeric alpha-helical coiled-coil configuration with five interchain disulfide bonds. The structure was solved using single-wavelength anomalous diffraction from sulfur-containing residues (sulfur SAD) to generate phase information. Biochemical and circular dichroism (CD) experiments show that a second PLVAP ECD segment also has a parallel dimeric alpha-helical configuration—presumably a coiled coil—held together with interchain disulfide bonds. Overall, ~2/3 of the ~390 amino acids within the PLVAP ECD adopt a helical configuration, as determined by CD. We also determined the sequence and epitope of MECA-32, an anti-PLVAP antibody. Taken together, these data lend strong support to the model of capillary diaphragms formulated by Tse and Stan in which approximately ten PLVAP dimers are arranged within each 60- to 80-nm-diameter opening like the spokes of a bicycle wheel. Passage of molecules through the wedge-shaped pores is presumably determined both by the length of PLVAP—i.e., the long dimension of the pore—and by the chemical properties of amino acid side chains and N-linked glycans on the solvent-accessible faces of PLVAP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

A twist grain boundary phase in aqueous solutions of the nucleic acid tetramer GTAC

At high concentration, long Watson/Crick (WC) double-helixed DNA forms columnar crystal or liquid crystal phases of linear, parallel duplex chains packed on periodic lattices. This can also be a structural motif of short NA oligomers such as the 5’-GTAC-3’ studied here, which makes four-base WC duplexes having hydrophobic blunt ends. End-to-end aggregation then assembles these duplexes into columns and columnar phases are stabilized, in spite of breaks in the double helix every four bases. But the new degrees of freedom introduced by such breaks also enable opportunities for a more diverse palette of self-assembly modes, producing striking self-assemblies of DNA that would not be achievable with contiguous polymers. These include recently reported three-dimensional (3D) periodic low-density nanoscale networks of GCCG, and the twist grain boundary (TGB) phase presented here. In the TGB, columns of GTAC pairs assemble into monolayer sheets in which the duplex columns are mutually parallel. However, unlike in the columnar crystals, these sheets stack in helical fashion into lamellar arrays in which the column axis of each layer is rotated through a 60° angle with respect to the columns in neighboring layers. This assembly of DNA is unique in that it the fills a 3D volume wherein the major grooves of columns in each layer mutually enter and interlock with the major grooves of columns in neighboring layers. This locking is optimized by small adjustments in structure enabled by the breaks in the duplex backbones.

DNA↗
DOE OSTI2023

Numerical Simulations of Flow-Induced Deflections in MITR LEU Fuel Plate Due to Channel Size Disparity

The hydromechanical stability of the fuel plates in parallel coolant channels of a Materials Testing Reactor (MTR) fuel element design is of great importance to the safety of research and test reactors. Previous analytical, experimental, and numerical efforts focused on parallel channels with the same or similar size; also, in the prior numerical simulations, the fuel plate was often assumed to be perfectly flat. This work presents the results of a fluid-structure interaction simulation performed to evaluate the flow-induced deflections of the fuel plates in the low-enriched uranium (LEU, <20 wt% 235 U) fuel element design for the conversion (from highly enriched uranium) of the Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR). Various manufacturing and assembly tolerances of the MITR LEU elements are considered in the analysis, and the effects of channel size disparity, nonideal plate shape, and flow rate uncertainty are investigated. Results show that, for all cases analyzed, the deflection occurs toward the larger channel, and the change in any channel stripe remains small (less than 0.021 mm) compared to fabrication tolerances. In addition to simulation work, a hydraulic performance test of the MITR LEU fuel element is currently planned to support conversion to the use of LEU fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Diffusion Synthetic Acceleration for Heterogeneous Domains, Compatible with Voids

A standard approach to solving the S N transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (a) preconditioned source iteration converges rapidly and (b) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using O(1) AMG iterations, convergence of the transport iteration typically requires 2 to 5× fewer iterations than current state-of-the-art “full” DSA, and the proposed method is trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3% of the mesh elements with five to ten AMG iterations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Soft particles and infinite-dimensional geometry

Abstract In the sigma model, soft insertions of moduli scalars enact parallel transport ofS-matrix elements about the finite-dimensional moduli space of vacua, and the antisymmetric double-soft theorem calculates the curvature of the vacuum manifold. We explore the analogs of these statements in gauge theory and gravity in asymptotically flat spacetimes, where the relevant moduli spaces are infinite-dimensional. These models have spaces of vacua parameterized by (trivial) flat connections on the celestial sphere, and soft insertions of photons, gluons, and gravitons parallel transportS-matrix elements about these infinite-dimensional manifolds. We argue that the antisymmetric double-soft gluon theorem ind + 2 bulk dimensions computes the curvature of a connection on the infinite-dimensional space Map ( S d , G ) / G , whereGis the global part of the gauge group. The analogous metrics in abelian gauge theory and gravity are flat, as indicated by the vanishing of the antisymmetric double-soft theorems in those models. In other words, Feynman diagram calculations not only know about the vacuum manifold of Yang–Mills theory, they can also be used to compute its curvature. The results have interesting implications for flat space holography.

Astronomy & Astrophysics↗
DOE OSTI2022

The theory of transport in helical spin-structure crystals

Abstract We study helical structures in spin-spiral single crystals. In the continuum approach for the helicity potential energy the simple electronic band splits into two non-parabolic bands. For low exchange integrals, the lower band is described by a surface with a saddle shape in the direction of the helicity axis. Using the Boltzmann equation with the relaxation due to acoustic phonons, we discover the dependence of the current on the angle between the electric field and helicity axis leading to the both parallel and perpendicular to the electric field components in the electroconductivity. The latter can be interpreted as a planar Hall effect. In addition, we find that the transition rates depend on an electron spin allowing the transition between the bands. The electric conductivities exhibit nonlinear behaviors with respect to chemical potential µ . We explain this effect as the interference of the band anisotropy, spin conservation, and interband transitions. The proposed theory with the spherical model in the effective mass approximation for conduction electrons can elucidate nonlinear dependencies that can be identified in experiments. We find the excellent agreement between the theoretical and experimental data for parallel resistivity depending on temperature at the phase transition from helical to ferromagnetic state in a M n P single crystal. In addition, we predict that the perpendicular resistivity abruptly drops to zero in the ferromagnetic phase.

Physics↗
DOE OSTI2020

Stepped pressure equilibrium with relaxed flow and applications in reversed-field pinch plasmas

The multi-region relaxed magnetohydrodynamics (MRxMHD) has been successful in the construction of equilibria in three-dimensional (3D) configurations. In MRxMHD, the plasma is sliced into sub-volumes separated by ideal interfaces, each undergoing relaxation, allowing the formation of islands and chaos. The resulting equilibrium has a stepped pressure profile across sub-volumes. The stepped pressure equilibrium code (SPEC) was developed to calculate MRxMHD equilibria numerically. In this work, we have extended the SPEC code to compute MRxMHD equilibria with field-aligned flow and rotation, following the theoretical development to incorporate cross-helicity and angular momentum constraints. The code has been verified for convergence and compared to a Grad–Shafranov solver in 2D. We apply our new tool to study the flow profile change before and after the sawtooth crash of a reversed-field pinch discharge, in which data of the parallel flow is available. We find the promising result that under the constraints of cross-helicity and angular momentum, the parallel flow profile in post-crash SPEC equilibrium is flat in the plasma core and the amplitude of the flow matches experimental observations. Lastly, we provide an example equilibrium with a 3D helical field structure as the favored lower energy state. This will be the first 3D numerical equilibrium in which the flow effects are self-consistently calculated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer

Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B → . The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e. the temperatures are assumed to be isotropic. Additionally, in certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the two-dimensional (2D) transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Relevance of E × B drifts for particle and heat transport in divertors

Radial electric fields up to ~4 kV m –1 are observed in the boundary between the private flux region (PFR) and the scrape-off layer (SOL) driving E × B drifts between the inner and outer targets at speeds up to 2.8 km s –1 in the Tokamak à configuration variable divertor. The resulting E × B fluxes, located in a narrow region ($\Delta {\rho _\Psi } < 0.012$ in normalized radius or $\Delta $ R – R sep <4 mm mapped to the outer midplane) are equivalent to around 20% of the total heat and particle flux to the divertor targets (inner + outer). At the peak E r , the E × B poloidal transport is equivalent to parallel flows with M ∥ ~ 3. In the snowflake divertor with a second X-point in the outer SOL, the drifts in the PFR-SOL boundary were equivalent to around 30% of the total heat and particle flux to the divertor targets and cover a region ~50% wider than in the single null ($\Delta {\rho _\Psi }$ ~ 0.018, $\Delta $R – R sep ~ 6 mm). Furthermore, the location of the PFR-SOL boundary drift shifts radially in the E ∥ × B direction when reversing the toroidal field direction. Peaks in density and electron pressure have been identified near the primary X-point along with large gradients in density, temperature, and potential, the latter resulting in a local electric field ~2.7 kV m –1 which drives a drift (1.9 km s –1 ) upwards towards the closed flux surfaces. Floating potential (V f ) magnitudes up to 75 V (~2 kT e ) were measured, indicating that V f and parallel currents should not be neglected when estimating plasma potential.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗