Search NASA⌕ Search

SEARCH · Search NASA

Results for “rotational diffusion”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Integrable symplectic maps with a polygon tessellation

Identifying integrable dynamics remains a formidable challenge, and despite centuries of research, only a handful of examples are known to date. In this article, we explore a distinct form of area-preserving (symplectic) mappings derived from the stroboscopic Poincaré cross section of a kicked rotator—an oscillator subjected to an external force periodically switched on in short pulses. The significance of this class of problems extends to various applications in physics and mathematics, including particle accelerators, crystallography, and studies of chaos. Notably, Suris's theorem constrains the integrability within this category of mappings, outlining potential scenarios with analytic invariants of motion. In this paper, we challenge the assumption of the analyticity of the invariant by exploring piecewise linear transformations on a torus ( T 2 ) and associated systems on the plane ( R 2 ), incorporating arithmetic quasiperiodicity and discontinuities. Introducing a new automated technique, we discovered previously unknown scenarios featuring polygonal invariants that form perfect tessellations and, moreover, fibrations of the plane or torus. This work reveals a novel category of planar tilings characterized by discrete symmetries that emerge from the invertibility of transformations and are intrinsically linked to the presence of integrability. Our algorithm relies on the analysis of the Poincaré rotation number and its piecewise monotonic nature for integrable cases, contrasting with the noisy behavior in the case of chaos, thereby allowing for clear separation. Some of the newly discovered systems exhibit the peculiar behavior of “integrable diffusion,” characterized by infinite and quasirandom hopping between tiles while being confined to a set of invariant segments. Finally, through the implementation of a smoothening procedure, all mappings can be generalized to quasi-integrable scenarios with suppressed volume occupied by chaotic trajectories, thereby opening doors to potential practical applications. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2022

Adatom-Driven Oxygen Intermixing during the Deposition of Oxide Thin Films by Molecular Beam Epitaxy

Thin film deposition from the vapor phase is a complex process involving adatom adsorption, movement, and incorporation into the growing film. While simulations have been able to capture key aspects of these processes to reproduce film morphology, these models typically ignore any subsurface diffusion that occurs. The results of these diffusion and intermixing processes are often observed experimentally, but cannot be explained by invoking bulk lattice diffusion at the deposition temperature. In this work, we present quantitative experimental data that reveals anion and cation intermixing over long length scales during the deposition of epitaxial Fe 2 O 3 and Cr 2 O 3 films and heterostructures by oxygen-plasma-assisted molecular beam epitaxy. We track this diffusion by incorporating well-defined tracer layers containing 18 O and/or 57 Fe, and measure their redistribution on the nanometer scale with atom probe tomography. Molecular dynamics simulations suggest potential intermixing events, which are then examined via nudged elastic band calculations. These reveal that adatoms on the film surface act to “pull up” subsurface O and Fe. Subsequent ring-like rotation mechanisms involving both adatom and subsurface anions then facilitate their mixing. A simple model incorporating these events qualitatively agrees with the experimental intermixing results. In addition to film deposition, these intermixing mechanisms may to be operant during other surface-mediated processes such as heterogeneous catalysis and corrosion.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Realistic galaxy image simulation via score-based generative models

ABSTRACT We show that a denoising diffusion probabilistic model (DDPM), a class of score-based generative model, can be used to produce realistic mock images that mimic observations of galaxies. Our method is tested with Dark Energy Spectroscopic Instrument (DESI) grz imaging of galaxies from the Photometry and Rotation curve OBservations from Extragalactic Surveys (PROBES) sample and galaxies selected from the Sloan Digital Sky Survey. Subjectively, the generated galaxies are highly realistic when compared with samples from the real data set. We quantify the similarity by borrowing from the deep generative learning literature, using the ‘Fréchet inception distance’ to test for subjective and morphological similarity. We also introduce the ‘synthetic galaxy distance’ metric to compare the emergent physical properties (such as total magnitude, colour, and half-light radius) of a ground truth parent and synthesized child data set. We argue that the DDPM approach produces sharper and more realistic images than other generative methods such as adversarial networks (with the downside of more costly inference), and could be used to produce large samples of synthetic observations tailored to a specific imaging survey. We demonstrate two potential uses of the DDPM: (1) accurate inpainting of occluded data, such as satellite trails, and (2) domain transfer, where new input images can be processed to mimic the properties of the DDPM training set. Here we ‘DESI-fy’ cartoon images as a proof of concept for domain transfer. Finally, we suggest potential applications for score-based approaches that could motivate further research on this topic within the astronomical community.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A high-order, localized-artificial-diffusivity method for Eulerian simulation of multi-material elastic-plastic deformation with strain hardening

A high-order method for Eulerian simulation of material undergoing large elastic–plastic deformation is developed. Thermodynamically consistent hyperelastic constitutive relations are assumed, facilitating the treatment of solids, liquids, and gases in a unified manner. Here, the method enables the simulation of multi-material interactions using a diffuse interface approach. Numerical capturing of material interfaces, shock waves, contact surfaces, and elastic-plastic strain discontinuities using high-order compact-difference schemes is assisted by Localized Artificial Diffusivity (LAD). In the new setting involving elastic–plastic deformation, the previously established terms for the artificial properties are verified to effectively regularize normal shocks. Additional LAD terms are introduced to the elastic and plastic kinematic equations to regularize shear shocks and other strain discontinuities, improving solution stability. Other important features of the method that improve robustness include the numerical treatment of compatibility terms in the kinematic equations, and the treatment of rotation. Particular emphasis is focused toward new advancements of the methods for plastic-deformation integration and the associated strain hardening of the material, including rate-dependent plasticity. The method is demonstrated on a variety of test problems, including 1-D impacts, a variant of the Shu-Osher problem, a Taylor impact, and a Richtmyer-Meshkov instability between two elastic–plastic solids with strain hardening.

42 ENGINEERING↗
DOE OSTI2020

Critical Practices in Rigorously Assessing the Inherent Activity of Nanoparticle Electrocatalysts

Accurate activity measurements for electrocatalytic materials are the backbone of impactful electrocatalyst research. Inherent to accurate measurements is the correct accounting of the active catalyst sites and the design of experiments to ensure that these sites are participating in the measured reaction. Improvements in electrocatalyst activity for fuel cell reactions (hydrogen oxidation and oxygen reduction) are essential for the widespread adoption of this carbon-neutral energy source. Activities for these half-reactions can be reported in a variety of ways including specific activities, mass activities, volumetric activities, and half-wave potentials. These values are traditionally measured with a rotating disk electrode (RDE) in which the electrocatalyst is supported on a porous active layer film attached to the RDE. To accurately measure electrocatalyst activity with an RDE, two major sources of errors must be addressed: (i) accurate electrochemical surface area (ECSA) quantification and (ii) fast diffusion of reactants through the active layer film. In this contribution, we first aim to detail through RDE experiments and mass transport modeling the potential errors that can be observed when these factors are not properly addressed. We then present recommendations and discuss techniques for properly constructing electrocatalyst thin films and characterizing the ECSA on monometallic and novel alloy electrocatalysts. These practices, when adopted, ensure greater confidence and reliability in reports of electrocatalyst activity.

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

Pinning dislocations in colloidal crystals with active particles that seek stacking faults

There is growing interest in functional, adaptive devices built from colloidal subunits of micron size or smaller. A colloidal material with dynamic mechanical properties could facilitate such microrobotic machines. Here we study via computer simulation how active interstitial particles in small quantities can be used to modify the bulk mechanical properties of a colloidal crystal. Passive interstitial particles are known to pin dislocations in metals, thereby increasing resistance to plastic deformation. We extend this tactic by employing anisotropic active interstitials that travel super-diffusively and bind strongly to stacking faults associated with partial dislocations. We find that: (1) interstitials that are effective at reducing plasticity compromise between strong binding to stacking faults and high mobility in the crystal bulk. (2) Reorientation of active interstitials in the crystal depends upon rotational transitions between high-symmetry crystal directions. (3) The addition of certain active interstitial shapes at concentrations as low as 60 per million host particles (0.006%) can create a shear threshold for dislocation migration. This work demonstrates how active materials in a dense matrix can locally sense their environment and lead to bulk property changes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model

Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Virtual Volumetric Additive Manufacturing (VirtualVAM)

Abstract Tomographic volumetric additive manufacturing (VAM) produces arbitrary 3D geometries by exposure of a rotating volume of photopolymer resin to tomographically‐patterned illumination. This enables high speed, layer‐less printing of parts from a wide range of photopolymers not amenable to layer‐by‐layer processes. Since the entire geometry is produced at once over the course of a few seconds to minutes, molecular diffusion length scales become significant to the printing process. Understanding these molecular reaction and diffusion processes is imperative for advancing VAM to a usable technology. These processes are experimentally very difficult to monitor and measure. Herein, VirtualVAM ‐ a simulation framework for modeling the tomographic VAM process, is developed and experimentally validated. VirtualVAM simulates reaction, diffusion, and heat generation processes over the course of a print with single‐voxel resolution. From a few experimentally‐determined input parameters and a set of images for projection, VirtualVAM is able to generate a large spatio‐temporal data set for any given tomographic VAM print. Using VirtualVAM, a number of experimentally‐unattainable aspects of the VAM process are investigated such as single‐voxel conversion profiles, effect of molecular oxygen, and stopping time determination. VirtualVAM also enables the optimization of exposure patterns to further improve contrast between in‐part and out‐of‐part delivered dose.

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

Gravito-turbulence and dynamo in poorly ionized protostellar discs – I. Zero-net-flux case

In their early stages, protoplanetary discs are sufficiently massive to undergo gravitational instability (GI). This instability is thought to be involved in mass accretion, planet formation via gas fragmentation, the generation of spiral density waves, and outbursts. A key and very recent area of research is the interaction between the GI and magnetic fields in young protoplanetary discs, in particular whether this instability is able to sustain a magnetic field via a dynamo. We conduct 3D, stratified shearing-box simulations using two independent codes, PLUTO and Athena++, to characterize the GI dynamo in poorly ionized protostellar discs subject to ambipolar diffusion. We find that the dynamo operates across a large range of ambipolar Elssaser number Am (which characterizes the strength of ambipolar diffusion) and is particularly strong in the regime Am = 10–100, with typical magnetic to thermal energy ratios of order unity. The dynamo is only weakly dependent on resolution (at least for Am ≲ 100), box size, and cooling law. The magnetic field is produced by the combination of differential rotation and large-scale vertical roll motions associated with spiral density waves. Our results have direct implications for the dynamo process in young protoplanetary discs and possibly some regions of active galactic nucleus discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Understanding Large-Scale Dynamos in Unstratified Rotating Shear Flows

For this work, we combine simulations with new analyses that overcome previous pitfalls to explicate how nonhelical mean-field dynamos grow and saturate in unstratified, magnetorotationally driven turbulence. Shear of the mean radial magnetic field amplifies the azimuthal component. Radial fields are regenerated by velocity fluctuations that induce shear of radial magnetic fluctuations, followed by Lorentz and Coriolis forces that source a negative off-diagonal component in the turbulent diffusivity tensor. We present a simple schematic to illustrate this dynamo growth. A different part of the Lorentz force forms a third-order correlator in the mean electromotive force that saturates the dynamo.

accretion disk & black-hole plasma↗
DOE OSTI2024

Quenched lattice fluctuations in optically driven SrTiO3

Abstract Crystal lattice fluctuations, which are known to influence phase transitions of quantum materials in equilibrium, are also expected to determine the dynamics of light-induced phase changes. However, they have only rarely been explored in these dynamical settings. Here we study the time evolution of lattice fluctuations in the quantum paraelectric SrTiO 3 , in which mid-infrared drives have been shown to induce a metastable ferroelectric state. Crucial in these physics is the competition between polar instabilities and antiferrodistortive rotations, which in equilibrium frustrate the formation of long-range ferroelectricity. We make use of high-intensity mid-infrared optical pulses to resonantly drive the Ti–O-stretching mode at 17 THz, and we measure the resulting change in lattice fluctuations using time-resolved X-ray diffuse scattering at a free-electron laser. After a prompt increase, we observe a long-lived quench in R-point antiferrodistortive lattice fluctuations. Their enhancement and reduction are theoretically explained by considering the fourth-order nonlinear phononic interactions to the driven optical phonon and third-order coupling to lattice strain, respectively. These observations provide a number of testable hypotheses for the physics of light-induced ferroelectricity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fully quantum algorithm for mesoscale fluid simulations with application to partial differential equations

Fluid flow simulations marshal our most powerful computational resources. In many cases, even this is not enough. Quantum computers provide an opportunity to speed up traditional algorithms for flow simulations. We show that lattice-based mesoscale numerical methods can be executed as efficient quantum algorithms due to their statistical features. This approach revises a quantum algorithm for lattice gas automata to reduce classical computations and state preparation at every time step. For this, the algorithm approximates the qubit relative phases and subtracts them at the end of each time step. Phases are evaluated using the iterative phase estimation algorithm and subtracted using single-qubit rotation phase gates. Further, this method optimizes the quantum resource required and makes it more appropriate for near-term quantum hardware. We also demonstrate how the checkerboard deficiency that the D1Q2 scheme presents can be resolved using the D1Q3 scheme. The algorithm is validated by simulating two canonical partial differential equations: the diffusion and Burgers' equations on different quantum simulators. We find good agreement between quantum simulations and classical solutions for the presented algorithm.

97 MATHEMATICS AND COMPUTING↗