Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Effects of drift on the transport of cosmic rays. VI - A three-dimensional model including diffusion

The first results are presented from a series of computer simulations of the solar modulation of galactic cosmic rays using a full three-dimensional model which incorporates all known important effects on particle transport, particle drifts, convection with the solar wind, energy loss, and anisotropic diffusion. The model is time-independent in the coordinate frame rotating with the sun, so corotating effects can be studied. Modulation in an interplanetary magnetic field model is considered in which the current sheet separating the northern and southern solar hemispheres is warped and corotating with the sun. The amplitude of the warp is varied to simulate possible solar cycle variation of the magnetic field. Substantial effects due to the warp of the current sheet are found. Comparison of the model results with various data is presented. Among other things, it is shown that the intensity may decrease away from the current sheet for both signs of the magnetic field, as suggested by recent observations, and in contrast with inferences from earlier, more approximate calculations.

Kota, J.↗
DOE OSTI2020

Reduced Artifact Approach for Determining Diffusion Coefficients in Time-Resolved Microscopy

Ultrafast microscopy methods traditionally assume a Gaussian profile to extract excited state diffusivities from transport measurements. Although this fitting method recovers accurate diffusion coefficients when the point spread function is well-represented by a Gaussian, even minor spatial aberrations introduced by the imaging system cause significant errors in the determined value. To provide a more accurate measure of excited state transport in nano- and microscale materials systems, in this work an alternative analysis protocol is proposed that numerically convolves the Green’s function solution to the diffusion equation with the experimentally measured point spread function. In contrast to the Gaussian fitting approach, the numerical convolution is shown to be robust against artifacts caused by nonideal point spread functions. Furthermore, the numerical convolution approach is highly effective at resolving anisotropic diffusion in modeled data.

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

On Thermodynamic and Kinetic Mechanisms for Stabilizing Surface Solid Solutions

Many processes for energy storage rely on transformations between phases with strong separation tendencies. In these systems, performance limitations can arise from undesirable chemical and mechanical factors associated with the phase separation behavior. Solid solutions represent a desirable alternative, provided the conditions for their formation are known. In the present work, we invoke linear stability theory and diffuse-interface mesoscopic simulations to demonstrate that solid solutions can be stabilized near surface layers of phase-separating systems. Two factors are found to drive surface solid-solution formation: surface relaxation of solution self-strain energy and anisotropy of diffusion mobility. Using a strongly phase-separating LiXFePO4 particle as a model system, we show that the relaxation of the solution self-strain energy competes against the relaxation of the coherency strain energy to stabilize surface solid solutions. Our theoretical understanding also suggests that highly anisotropic diffusion mobility can provide an alternative kinetic route to achieve the same aim, with stabilizing behavior strongly dependent on the specific alignment of the surface orientation. Our findings provide fundamental guidance for manipulating solid-solution behavior in nanoscale structures, in which surface effects become especially significant. Beyond energy storage materials, our findings have important implications for understanding solid-solution formation in other phase-separating systems from metal alloys to ceramics.

25 ENERGY STORAGE↗
DOE OSTI2025

Frustrated charge density wave and quasi-long-range bond-orientational order in the magnetic kagome FeGe

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and the multiple-q charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and obscured by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-q CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector q = q M at T CDW < T < T*(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic B-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Subires, D. [Donostia International Physics Center↗
DOE OSTI2014

Large-scale anisotropy of TeV-band cosmic rays

The expected anisotropy in the 1 to 10{sup 4} TeV energy range is calculated for Galactic cosmic rays with both anisotropy in the diffusion tensor and source discreteness taken into account. We find that if the sources are distributed radially (but with azimuthal symmetry) in proportion to Galactic pulsars, the expected anisotropy almost always exceeds the observational limits by one order of magnitude in the case of isotropic diffusion. If the radial diffusion is more than an order of magnitude smaller than the azimuthal diffusion rate, the radial gradient of the sources can be accommodated about 5% of the time. If the sources are concentrated in the spiral arms, then the anisotropy depends on our location between them, but in some spatial window, roughly equidistant from adjacent spiral arms, the observational constraints on anisotropy are obeyed roughly 20%-30% of the time for extremely anisotropic diffusion. The solar system is in that window less than 10% of the time, but it may be there now. Under the assumption of isotropic diffusion, nearby supernovae are found to produce a discreteness anisotropy that is nearly two orders of magnitude in excess of the observational limit if all supernovae are assumed to contribute equally with a source rate 1 in every 100 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Atomic-Scale Mechanism of Unidirectional Oxide Growth

A basic knowledge of the unidirectional growth mechanisms is required for precise control on size, shape, and thereby functionalities of nanostructures. The oxidation of many metals results in oxide nanowire growth with a bicrystal grain boundary along the axial direction. Using transmission electron microscopy that spatially and temporally resolves CuO nanowire growth during the oxidation of copper, here we provide direct evidence of the correlation between unidirectional crystal growth and bicrystal grain boundary diffusion. Based on atomic scale observations of the upward growth at the nanowire tip, oscillatory downward growth of atomic layers on the nanowire sidewall and the parabolic kinetics of lengthening, bicrystal grain boundary diffusion is the mechanism by which Cu ions are delivered from the nanowire root to the tip. Together with density-functional theory calculations, we further show that the asymmetry in the corner-crossing barriers promotes the unidirectional oxide growth by hindering the transport of Cu ions from the nanowire tip to the sidewall facets. We anticipate the broader applicability of these results in manipulating the growth of nanostructured oxides by controlling the bicrystal grain boundary structure that favors anisotropic diffusion for unidirectional, one-dimensional crystal growth for nanowires or isotropic diffusion for two-dimensional platelet growth.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Leith diffusion model for homogeneous anisotropic turbulence

Here, a proposal for a spectral closure model for homogeneous anisotropic turbulence. The systematic development begins by closing the third-order correlation describing nonlinear interactions by an anisotropic generalization of the Leith diffusion model for isotropic turbulence. The correlation tensor is then decomposed into a tensorially isotropic part, or directional anisotropy, and a trace-free remainder, or polarization anisotropy. The directional and polarization components are then decomposed using irreducible representations of the SO(3) symmetry group. Under the ansatz that the decomposition is truncated at quadratic order, evolution equations are derived for the directional and polarization pieces of the correlation tensor. Here, numerical simulation of the model equations for a freely decaying anisotropic flow illustrate the non-trivial effects of spectral dependencies on the different return-to-isotropy rates of the directional and polarization contributions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Diffuse scattering from dynamically compressed single-crystal zirconium following the pressure-induced α → ω phase transition

The prototypical α → ω phase transition in zirconium is an ideal test bed for our understanding of polymorphism under extreme loading conditions. After half a century of study, a consensus had emerged that the transition is realized via one of two distinct displacive mechanisms, depending on the nature of the compression path. However, recent dynamic-compression experiments equipped with diffraction diagnostics performed in the past few years have revealed new transition mechanisms, demonstrating that our understanding of the underlying atomistic dynamics and transition kinetics is in fact far from complete. We present classical molecular dynamics simulations of the α → ω phase transition in single-crystal zirconium shock compressed along the [0001] axis using a machine-learning-class potential. The transition is predicted to proceed primarily via a modified version of the two-stage Usikov-Zilberstein mechanism, whereby the high-pressure ω phase heterogeneously nucleates at boundaries between grains of an intermediate β phase. We further observe the fomentation of atomistic disorder at the junctions between β grains, leading to the formation of highly defective interstitial material between the ω grains. We directly compare synthetic x-ray diffraction patterns generated from our simulations with those obtained using femtosecond diffraction in recent dynamic-compression experiments, and show that the simulations produce the same unique, anisotropic diffuse scattering signal unlike any previously seen from an elemental metal. Our simulations suggest that the diffuse signal arises from a combination of thermal diffuse scattering, nanoparticlelike scattering from residual kinetically stabilized α and β grains, and scattering from interstitial defective structures. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2014

Effects of anisotropies in turbulent magnetic diffusion in mean-field solar dynamo models

We study how anisotropies of turbulent diffusion affect the evolution of large-scale magnetic fields and the dynamo process on the Sun. The effect of anisotropy is calculated in a mean-field magnetohydrodynamics framework assuming that triple correlations provide relaxation to the turbulent electromotive force (so-called the 'minimal τ-approximation'). We examine two types of mean-field dynamo models: the well-known benchmark flux-transport model and a distributed-dynamo model with a subsurface rotational shear layer. For both models, we investigate effects of the double- and triple-cell meridional circulation, recently suggested by helioseismology and numerical simulations. To characterize the anisotropy effects, we introduce a parameter of anisotropy as a ratio of the radial and horizontal intensities of turbulent mixing. It is found that the anisotropy affects the distribution of magnetic fields inside the convection zone. The concentration of the magnetic flux near the bottom and top boundaries of the convection zone is greater when the anisotropy is stronger. It is shown that the critical dynamo number and the dynamo period approach to constant values for large values of the anisotropy parameter. The anisotropy reduces the overlap of toroidal magnetic fields generated in subsequent dynamo cycles, in the time-latitude 'butterfly' diagram. If we assume that sunspots are formed in the vicinity of the subsurface shear layer, then the distributed dynamo model with the anisotropic diffusivity satisfies the observational constraints from helioseismology and is consistent with the value of effective turbulent diffusion estimated from the dynamics of surface magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Robust Implicit Adaptive Low Rank Time-Stepping Methods for Matrix Differential Equations

In this work, we develop implicit rank-adaptive schemes for time-dependent matrix differential equations. The dynamic low rank approximation (DLRA) is a well-known technique to capture the dynamic low rank structure based on Dirac–Frenkel time-dependent variational principle. In recent years, it has attracted a lot of attention due to its wide applicability. Our schemes are inspired by the three-step procedure used in the rank adaptive version of the unconventional robust integrator (the so called BUG integrator) (Ceruti et al. in BIT Numer Math 62(4):1149–1174, 2022) for DLRA. First, a prediction (basis update) step is made computing the approximate column and row spaces at the next time level. Second, a Galerkin evolution step is invoked using an implicit solves for the small core matrix. Finally, a truncation is made according to a prescribed error threshold. Since the DLRA is evolving the differential equation projected on to the tangent space of the low rank manifold, the error estimate of the BUG integrator contains the tangent projection (modeling) error which cannot be easily controlled by mesh refinement. This can cause convergence issue for equations with cross terms. To address this issue, we propose a simple modification, consisting of merging the row and column spaces from the explicit step truncation method together with the BUG spaces in the prediction step. In addition, we propose an adaptive strategy where the BUG spaces are only computed if the residual for the solution obtained from the prediction space by explicit step truncation method, is too large. Here, we prove stability and estimate the local truncation error of the schemes under assumptions. We benchmark the schemes in several tests, such as anisotropic diffusion, solid body rotation and the combination of the two, to show robust convergence properties.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Phase transformation mechanism in irradiation-induced superlattice formation

Atomic kinetic Monte Carlo simulations were used to model void superlattice formation under irradiation in molybdenum, driven by anisotropic diffusion of self-interstitial atoms. A change in the phase transformation mechanism from nucleation and growth to spinodal decomposition occurred with increasing dose rate, with both mechanisms leading to superlattice formation. Analysis of a rate-theory based analytical model showed that an observed change in the kinetics of vacancy accumulation, the appearance of a region of positive second derivative in the plot of average vacancy concentration versus time, was caused by the onset of spinodal instability. Further, the analytical model showed that for molybdenum and several other metals where void superlattice formation is commonly observed, the phase transformation likely occurs by nucleation and growth. However, nickel may offer the possibility of experimental observation of the transition between phase transformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unveiling the lithium-ion transport mechanism in Li{sub 2}ZrCl{sub 6} Solid-State Electrolyte {ital via} deep learning-accelerated molecular dynamics simulations.

Lithium zirconium chlorides (LZCs) present a promising class of cost-effective solid electrolytes for next-generation all-solid-state batteries. The unique crystal structure of LZCs plays a crucial role in facilitating lithium-ion mobility, which further affects the electrochemical performance. To understand the underlying mechanism governing ion transport, we employed deep learning-accelerated molecular dynamics simulation on Li2ZrCl6 (trigonal alpha- and monoclinic beta-LZC), focusing specifically on the zirconium coordination environment. Our results reveal that disordered alpha-LZC exhibits the highest ionic conductivity, while beta-LZC demonstrates significantly lower conductivity, closely aligning with experimental findings. The study confirms that across all phases, lithium migration proceeds via the site-to-site hopping mechanism, where variations in site residence times critically impact the overall ionic conductivity. In alpha-LZCs, lithium ions prefer to anisotropically diffuse across interlayers as the result of a lower energy barrier, driven primarily by collective diffusion. In contrast, lithium ions in beta-LZC primarily isotropically diffuse within the intralayer, hindered by higher energy barriers and determined by individual diffusion. The variation in ZrCl6 2- octahedral unit softening, induced by the specific layered arrangement of zirconium atoms, emerges as a critical determinant of the energy barriers across the LZC phases. These atomic-scale insights into the transport processes provide valuable guidance for the rational design and optimization of LZCs-based electrolytes, accelerating their practical application in advanced energy storage technologies.

Guo, Hanzeng↗
DOE OSTI2015

Thermodynamic calculations of Fe–Mg interdiffusion in (Mg,Fe){sub 2}SiO{sub 4} polymorphs and perovskite

In this study, we show that the temperature and pressure dependence of Fe–Mg interdiffusion in (Fe,Mg){sub 2}SiO{sub 4} polymorphs (olivine, wadsleyite, and ringwoodite) and perovskite can be successfully reproduced in terms of bulk elastic and expansivity data through a thermodynamic model (so-called cBΩ model) that interconnects point defect parameters with bulk properties. Under dry and wet conditions, our calculated Fe–Mg interdiffusion coefficients D{sub calc}{sup Fe-Mg} (particularly for anisotropic diffusivity in olivine), activation enthalpy h{sup act}, and activation volume υ{sup act} over a wide range of geologically relevant temperatures (1000–2400 K) and pressures (0–100 GPa) are consistent with the experimental ones when the uncertainties are considered.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Energy-conserving contact dynamics of nonspherical rigid-body particles

Understanding the contact dynamics of nonspherical particles is crucial for accurately modeling colloidal and granular systems where shape anisotropy dictates structural organization and transport properties. We here introduce an energy-conserving contact dynamics framework for arbitrary convex rigid-body particles by implementing vertex–boundary interactions in 2D and vertex–surface and edge–edge detection in 3D. The established formulation enables continuous force evaluation and prevents particle overlap while conserving total energy during translational and rotational motion. We demonstrate the framework’s stability and its utility to capture packing behavior, anisotropic diffusion, and equations of state of polygonal and polyhedral particles as examples. The framework establishes a robust and extensible foundation for investigating nonequilibrium dynamics of complex nonspherical particulate systems, enabling enhanced understanding of phenomena across spatiotemporal scales in self- and directed-assembly, granular flows, and hydrodynamics, potentially coupled with interactions that represent underlying physical mechanisms.

Discrete element method↗
DOE OSTI2017

Cosmic-ray propagation with DRAGON2: I. numerical solver and astrophysical ingredients

We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions. The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios. The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well. This new version facilitates the users to include their own physical models by means of a modular C++ structure.

79 ASTRONOMY AND ASTROPHYSICS↗