Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Findings”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mitigation of resistive drift wave and ion temperature gradient instabilities by velocity shear

The effects of velocity shear on the resistive drift wave instability in the non-adiabatic limit and the toroidal ion temperature gradient instability are investigated for a plasma of inhomogeneous density/ion temperature, respectively. For the resistive drift wave, we find that the instability growth rate decreases monotonically with increasing magnitude of shear, but we find that complete stabilization is impossible. For the ion temperature gradient instability, we find that the standard WKB approximation is insufficient to describe the full behavior of the instability, and that an analysis of the localized eigenmode problem reveals two separate unstable solutions, which the WKB approximation does not predict. The impact of flow shear on these two new unstable solutions is discussed. In both resistive drift wave and ion temperature gradient instabilities, the sheared flow causes a shifting, tilting, and sharpening of the electrostatic potential eddies.

Physics↗
DOE OSTI2024

Filamentary velocity scaling validation and spin dynamics in the DIII-D tokamak

Measured filament velocities in the DIII-D tokamak are compared against theoretical scalings, finding that the latter often represents an upper limit on experimental velocity distributions with most filaments possessing lower velocity. Filament spin from internal E × B drift is experimentally demonstrated to alter filament radial velocity. A critical spin velocity, where filament radial velocity peaks, is observed and corresponds to approximately 5 km/s. This transition is corroborated using a less direct measure of filament spin in the form of a temperature ratio. These techniques are combined to find that the critical spin velocity closely aligns with transport times along and across filaments becoming comparable. The normalized filament size distribution is consistent with the most stable size as dictated by Kelvin–Helmholtz and curvature-driven instabilities. Overall, the findings suggest filament stability and spin alter filamentary transport that may threaten the integrity of first walls in fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2020

Self-interaction error overbinds water clusters but cancels in structural energy differences

Here, we gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi–Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew–Burke–Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction–corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/ H 2 O from about 42 meV/ H 2 O for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction–corrected SCAN. These findings suggest that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

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

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

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

A qualitative study of marginalized students’ academic, physical, and social self-efficacy in a multiweek geoscience field program

Undergraduate summer field programs are valuable experiences that can foster or reduce students’ self-efficacy, an important factor in students’ success and retention in geoscience. Growing research findings show that science field experiences can be hostile and unwelcoming to students with marginalized identities, which may negatively impact their self-efficacy in geoscience, a discipline with a dearth of students from underrepresented, marginalized identities. We conducted an interpretive qualitative study examining how summer geoscience field programs affected two undergraduate, marginalized students’ self-efficacy. Adding to existing theoretical explanations of self-efficacy, we identified three types of self-efficacy impacted positively and negatively by geoscience field experiences: academic, physical, and social self-efficacy. We developed a nuanced understanding of the specific field experiences that influenced the ‘ups and downs’ of students’ self-efficacy and, ultimately, their intent in continuing to pursue a geoscience education or career. Despite negative experiences, including gender discrimination, crude sexual jokes, and a lack of belonging, the students described their intent to persist in geoscience. Here, our findings can assist geoscience educators (and others in field-based sciences) to consider experiences that support and hinder marginalized students’ self-efficacy. Also, our findings can guide efforts to improve geoscience field programs to create more inclusive environments.

case study↗
DOE OSTI2025

The CanBikeCO Full Pilot: Long-Term Results and Analysis From an E-Bike Program in Colorado, USA

Personal micromobility devices like bicycles, e-bikes, and scooters are low- or zero-energy alternatives to single-occupancy vehicles. However, a lack of data has led to a dearth of data-driven research on personally owned e-bike usage. We present longitudinal findings from the CanBikeCO program, focused on e-bike adoption and use across demographics, trip characteristics, and geographies in the state of Colorado. CanBikeCO recorded travel survey data from low-income individuals provided with personal e-bikes by the Colorado Energy Office in six communities across Colorado from July 2021 to December 2022. The data were collected using a custom instance of the National Renewable Energy Laboratory OpenPATH platform, which combines passive data collection with semantic information such as trip mode and purpose labels. To our knowledge, there are no prior travel survey data on personally owned e-bikes with this range and scope. Insights from this unique dataset include: (i) work trips were 17% more likely than average trips to be taken on an e-bike, (ii) e-bikes were most often reported to replace cars (34% of e-bike trips) and other personal micromobility devices (22%), and (iii) participants favored walking for trips less than 1 mile, e-bikes for trips of 1-3 miles, and e-bikes, cars, or shared rides for trips of 3-20 miles. The data used to generate these results have been made available in the Transportation Secure Data Center. We find e-bike use is appealing across age groups and may be related to characteristics of land use, urban form, occupation, income, and car ownership. We conclude for this population that the energy demand added by e-bike use (induced demand and replacing non-motorized modes) is outweighed by the reduction in energy demand from replacement of single-occupancy vehicle trips with e-bike trips. Our findings suggest considerable potential for energy savings from personal e-bike ownership.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Multi-Fidelity Bayesian Optimization with Gaussian Processes for Double Shell Inertial Confinement Fusion Target Design

Reliable, secure access to energy is a major focus for national security efforts. One potential route to such energy is through fusion reactions in inertial confinement fusion (ICF) experiments. Such experiments are carried out at facilities such as the National Ignition Facility (NIF) in Livermore, California, where high powered lasers are used to compress a DT fuel-containing target to the necessary high temperature, high pressure conditions. These experiments are limited in number, which creates a heavy dependence on high fidelity predictive physics simulations and analysis performed “pre shot,” or before the experiment occurs. Many of these simulations in higher dimensions (2D and 3D) are computationally expensive, so finding optimal simulation-based designs presents its own challenges. In this work, we present our multi-fidelity Bayesian optimization with Gaussian processes (GPs) for ICF double shell targets, where a 1D surrogate model is used to help find a 2D surrogate model, enabling us to find optimal targets in the higher fidelity (2D), while saving computational cost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Late-time post-merger modeling of a compact binary: effects of relativity, r-process heating, and treatment of transport

Detectable electromagnetic counterparts to gravitational waves from compact binary mergers can be produced by outflows from the black hole-accretion disk remnant during the first 10 s after the merger. Two-dimensional axisymmetric simulations with effective viscosity remain an efficient and informative way to model this late-time post-merger evolution. In addition to the inherent approximations of axisymmetry and modeling turbulent angular momentum transport by a viscosity, previous simulations often make other simplifications related to the treatment of the equation of state and turbulent transport effects. In this paper, we test the effect of these modeling choices. By evolving with the same viscosity the exact post-merger initial configuration previously evolved in Newtonian viscous hydrodynamics, we find that the Newtonian treatment provides a good estimate of the disk ejecta mass but underestimates the outflow velocity. We find that the inclusion of heavy nuclei causes a notable increase in ejecta mass. An approximate inclusion of r-process effects has a comparatively smaller effect, except for its designed effect on the composition. Diffusion of composition and entropy, modeling turbulent transport effects, has the overall effect of reducing ejecta mass and giving it a speed with lower average and more tightly-peaked distribution. Also, we find significant acceleration of outflow even at distances beyond 10 000 km, so that thermal wind velocities only asymptote beyond this radius and at higher values than often reported.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Skyrmion transport and annihilation in funnel geometries

Using atomistic simulations, we have investigated the transport and annihilation of skyrmions interacting with a funnel array under a current applied perpendicular to the funnel axis. We find that transport without annihilation is possible at low currents, when the motion is dominated by skyrmion–skyrmion interactions and skyrmions push each other through the funnel opening. Skyrmion annihilation occurs for higher currents when skyrmions in the upper half of the sample exert pressure on skyrmions in the bottom half of the sample due to the external current. Upon interacting with the funnel wall, the skyrmions undergo a size reduction that makes it easier for them to pass through the funnel opening. We find five phases as a function of the applied current and the size of the funnel opening: (i) pinned, (ii) transport without annihilation, (iii) transport with annihilation, (iv) complete annihilation, and (v) a reentrant pinning phase that only occurs for very narrow openings. Finally, our findings provide insight into how to control skyrmion transport using funnel arrays by delineating regimes in which transport of skyrmions is possible as well as the conditions under which annihilation occurs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spin–orbit torque in a three-fold-symmetric bilayer and its effect on magnetization dynamics

Field-free switching of perpendicular magnetization has previously been observed in an epitaxial L1 1 -ordered CoPt/CuPt bilayer and attributed to spin–orbit torque (SOT) arising from the crystallographic 3m point group of the interface. Using a first-principles nonequilibrium Green's function formalism combined with the Anderson disorder model, we calculate the angular dependence of the SOT in a CoPt/CuPt bilayer and find that the magnitude of the 3m SOT is about 20% of the conventional dampinglike SOT. We further study the magnetization dynamics in perpendicularly magnetized films in the presence of 3m SOT and Dzyaloshinskii–Moriya interaction, using the equations of motion for domain wall dynamics and micromagnetic simulations. For systems where strong interfacial DMI results in the Néel character of domain walls, we find that a very large current density is required to achieve deterministic switching, because reorientation of the magnetization inside the domain wall is necessary to induce the switching asymmetry. For thicker films with relatively weak interfacial DMI and the Bloch character of domain walls, deterministic switching is possible at much smaller current densities, which may explain the recent experimental findings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Multiferroic ground states in free standing perovskite-based nanodots: a density functional theory study

In this work, we use density functional theory to investigate the possibility of polar and multiferroic states in free-standing, perovskite-based nanodots at the atomic limit of miniaturization: single unit cells with terminations which allow centro-symmetry. We consider both A-O and B-O 2 terminations for three families of nanodots: (i) A = Ba with B = Ti, Zr, and Hf; (ii) A = Ca and Sr with B = Ti; and (iii) A = Na and K with B = Nb. We find all A–O terminated dots to be non-polar and to exhibit cubic symmetry (except for K 8 NbO 6 ), regardless of the presence of ferroelectricity in the bulk. In contrast, all the B–O 2 terminated nanodots considered relax to a non-cubic ground state. Rather surprisingly, all of these structures exhibit polar ground states (except NaNb 8 O 12 ). We propose a new structural parameter, the cluster tolerance factor (CTF), to determine whether a particular chemistry will result in a polar ground state nanodot, analogous to the Goldschmidt factor for bulk ferroelectrics. In addition, we find that all A–O terminated (except Ca 8 TiO 6 ) and all polar B–O 2 terminated nanodots are magnetic, where none show magnetism in the bulk. As with bulk systems, multiferroicity in the B–O 2 terminated dots originates from separation between spin density in peripheral B atoms and polarity primarily caused by the off-center central A atom. Our findings stress that surface termination plays a crucial role in determining whether ferroelectricity is completely suppressed in perovskite-based materials at their limit of miniaturization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2025

Assessing the numerical stability of physics models to equilibrium variation through database comparisons on DIII-D

High fidelity kinetic equilibria are crucial for tokamak modeling and analysis. Manual workflows for constructing kinetic equilibria are time consuming and subject to user error, motivating development of automated equilibrium reconstruction tools to provide accurate and consistent reconstructions for downstream physics analysis. These automated tools also provide access to kinetic equilibria at large database scales, which enables the quantification of general uncertainties arising from equilibrium reconstruction techniques. In this paper, we compare a large database of DIII-D kinetic equilibria generated manually by physics experts to equilibria from automated kinetic reconstruction tools, assessing the impact of reconstruction method on equilibrium parameters and resulting magnetohydrodynamic stability calculations. We find agreement among scalar parameters, whereas profile quantities, such as the bootstrap current, show larger disagreements. We analyze ideal kink and classical tearing stability with DCON and STRIDE respectively, finding that the kink stability calculation is generally more robust than the tearing index Δ' calculation. We find that in 90% of cases, both kink stability classifications are unchanged between the manual expert and automated kinetic equilibria.

CAKE↗
DOE OSTI2020

Skyrmion dynamics and topological sorting on periodic obstacle arrays

Abstract We examine skyrmions under a dc drive interacting with a square array of obstacles for varied obstacle size and damping. When the drive is applied in a fixed direction, we find that the skyrmions are initially guided in the drive direction but also move transverse to the drive due to the Magnus force. The skyrmion Hall angle, which indicates the difference between the skyrmion direction of motion and the drive direction, increases with drive in a series of quantized steps as a result of the locking of the skyrmion motion to specific symmetry directions of the obstacle array. On these steps, the skyrmions collide with an integer number of obstacles to create a periodic motion. The transitions between the different locking steps are associated with jumps or dips in the velocity–force curves. In some regimes, the skyrmion Hall angle is actually higher than the intrinsic skyrmion Hall angle that would appear in the absence of obstacles. In the limit of zero damping, the skyrmion Hall angle is 90°, and we find that it decreases as the damping increases. For multiple interacting skyrmion species in the collective regime, we find jammed behavior at low drives where the different skyrmion species are strongly coupled and move in the same direction. As the drive increases, the species decouple and each can lock to a different symmetry direction of the obstacle lattice, making it possible to perform topological sorting in analogy to the particle sorting methods used to fractionate different species of colloidal particles moving over two-dimensional obstacle arrays.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Polarized QED cascades

By taking the spin and polarization of the electrons, positrons and photons into account in the strong-field QED processes of nonlinear Compton emission and pair production, we find that the growth rate of QED cascades in ultra-intense laser fields can be substantially reduced. While this means that fewer particles are produced, we also found them to be highly polarized. We further find that the high-energy tail of the particle spectra is polarized opposite to that expected from Sokolov–Ternov theory, which cannot be explained by just taking into account spin-asymmetries in the pair production process, but results significantly from 'spin-straggling'. We employ a kinetic equation approach for the electron, positron and photon distributions, each of them spin/polarization-resolved, with the QED effects of photon emission and pair production modelled by a spin/polarization dependent Boltzmann-type collision operator. For photon-seeded cascades, depending on the photon polarization, we find an excess or a shortage of particle production in the early stages of cascade development, which provides a path towards a controlled experiment. Throughout this paper we focus on rotating electric field configuration, which represent an idealized model and allows for a straightforward interpretation of the observed effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗