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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 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 OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1971

Transport of solar flare protons: Comparison of a new analytic model with spacecraft measurements

An analytic solution was obtained to the complete Fokker-Planck equation including the effects of convection, interplanetary deceleration and acceleration, corotation, and anisotropic diffusion. The solution yields a time to maximum for the particle flux of about 10 h and an exponential decay time of about 5 h. Several solar flare particle events were observed with the solar and galactic cosmic ray experiment on OGO-6. Comparisons of the calculated time dependence of the fluxes with these observations of 1 to 70 MeV protons show that the model adequately describes both the rise and decay times.

Lupton, J. E.↗
NASA NTRS1971

The propagation of solar cosmic-ray bursts.

Description of a model showing analytically the three phases of anisotropy which occur during solar cosmic-ray events in the 7.5 to 21 MeV kinetic-energy interval and reported by McCracken et al. (1971): (1) a highly anisotropic, near field-aligned, initial phase, (2) a convective phase, and (3) a late-time phase in which the anisotropy is approximately perpendicular to the mean interplanetary magnetic field. The model is based on the cosmic-ray particles being convectively transported out from the sun, undergoing anisotropic diffusion along the interplanetary magnetic-field lines, and losing energy by adiabatic deceleration or by collision processes. The event is seen simply as a pulse moving outward from the sun after a cosmic-ray burst with a negative density-gradient in front of it and a positive gradient behind.

Ng, C. K.↗
NASA NTRS1981

Temporal Aperture Modulation

The two types of modulation techniques useful to X-ray imaging are reviewed. The use of optimum coded temporal aperature modulation is shown, in certain cases, to offer an advantage over a spatial aperture modulator. Example applications of a diffuse anisotropic X-ray background experiment and a wide field of view hard X-ray imager are discussed.

Proctor, R. J.↗
NASA NTRS1985

The expected cosmic ray density and stream distributions at the heliolatitudinal asymmetry of solar wind

The results of the spatial distribution of cosmic ray density, gradients, and anisotropy obtained on the basis of the numerical solution of the anisotropic diffusion equation with an account of solar wind velocity change depending on the latitudinal angle theta of the form U=u sub oe sup alpha theta and the diffusion coefficient depending on the spatial coordinates and the particle rigidity are presented. It is shown that the increase of the solar wind velocity and the diffusion coefficient with heliolatitude leads to gradient distributions that are in accord with experimental data observed in space. The results of the energetic spectrum of 11 and 22-year cosmic ray variations obtained with an account of direction of the general magnetic field of the Sun are presented are given.

Alania, M. V.↗
NASA NTRS1985

Spatial distribution of cosmic rays in magnetic cycles of the Sun

The space-time distribution of cosmic ray (CR) density is determined from the solution of the anisotropic diffusion equation. The boundary problem is solved in a quasistationary axisymmetric case, taking into account the delay of cosmic rays relative to the cyclic variations of solar activity. The specificty of the modulation model is the use of the empirical dependence of the free path of particles for scattering on the observed latitudinal distribution of solar activity. Another specific feature of the model is an account of an additional effect due to the drift of charged particles towards the helio-equator plane or backward, depending on the direction of the total magnetic field of the Sun.

Alania, M. V.↗
NASA NTRS1985

The modulation features of the long-period cosmic ray variations in connection with the sign change of the general magnetic field of the Sun

On the basis of the model and experimental investigations, the spatial distribution of cosmic ray anisotropy for the different epochs of solar activity is studied. A solution is offered to the anisotropic diffusion equation with regard to the electromagnetic conditions for the periods of minimum and maximum solar activity and cosmic ray particle drift in a regular interplanetary magnetic field. It is shown that the long period changes amplitude and phase of the diurnal variations of cosmic rays is limited not only by convection and diffusion of particles but also by the drift effect before and after the sing change of the general magnetic field of the Sun. The calculated model is compared with the results obtained on the basis of an analysis of the experimental data from the neutron super monitor station, and it is shown that for the periods when the lines of magnetic force of the Sun come from the Northern Hemisphere the phase of the first harmonic diurnal variation is shifted forwards to an earlier time.

Iskra, K.↗
NASA NTRS1985

The theoretical and experimental investigation of cosmic ray Forbush-effects

The theoretical results of analysis of the expected spatial distributions of density, gradients and anisotropy of cosmic rays, obtained on the basis of the numerical solution of anisotropic diffusion equation in the presence of the disturbances of shock wave type in the interplanetary space are presented. The theoretical calculations on the definition of the energy spectrum and anisotropy of galactic cosmic rays during Forbush effect are compared with the experimental data.

Alania, M. V.↗
NASA NTRS1995

Calculations of Diffuser Flows with an Anisotropic K-Epsilon Model

A newly developed anisotropic K-epsilon model is applied to calculate three axisymmetric diffuser flows with or without separation. The new model uses a quadratic stress-strain relation and satisfies the realizability conditions, i.e., it ensures both the positivity of the turbulent normal stresses and the Schwarz' inequality between any fluctuating velocities. Calculations are carried out with a finite-element method. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to ensure numerical credibility of the solutions. The standard K-epsilon model is also used in order to highlight the performance of the new model. Comparison with the experimental data shows that the anisotropic K-epsilon model performs consistently better than does the standard K-epsilon model in all of the three test cases.

Zhu, J.↗
NASA NTRS1992

Model atmospheres and radiation of magnetic neutron stars. I - The fully ionized case

Model neutron star atmospheres are calculated for typical cooling stars with a strong magnetic field and effective temperatures of 10 exp 5 to 10 exp 6 K. The effect of anisotropic photon diffusion in two normal modes are examined under the assumption that the opacity is due solely to the bremsstrahlung and Thomson scattering processes under conditions of LTE that are expected to prevail at the temperatures and densities obtained. The main aspects of anisotropic photon diffusion, and an original procedure for calculating model atmospheres and emitted spectra are discussed. Representative calculated spectra are given, and it is found that the hard spectral excess characterizing the nonmagnetic case, while still present, becomes less prominent in the presence of magnetic fields in the range of 10 exp 11 to 10 exp 13 G.

Shibanov, Iu. A.↗
DOE OSTI2025

Techniques for improved statistical convergence in quantification of eddy diffusivity moments

While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪⁡(100) simulations rather than 𝒪⁡(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.

general physics↗
DOE OSTI2026

Chemo-Mechanics of α-V 2 O 5 During Lithiation and Implications for Rechargeable Battery Cathodes

Chemo-mechanical degradation of layered oxide electrodes is strongly influenced by crystallographic anisotropy, local stress evolution, and ion insertion, yet the intrinsic mechanical response of layered materials remains incompletely understood. Indeed, most prior studies have focused on polycrystalline materials but single crystals enable direct observation of coupling between anisotropic ion diffusion and mechanical response. This study aims to determine how crystallographic anisotropy and lithiation affect deformation, fracture, and mechanical properties in single-crystal V 2 O 5 , and compares this behavior with polycrystalline counterparts. Polycrystalline V 2 O 5 thin films and single-crystal α-V 2 O 5 were studied using nanoindentation, scanning electron microscopy, focused ion beam cross-sectioning, and Raman spectroscopy. Single crystals were tested in pristine and chemically lithiated states, including experiments in which crystals were first plastically deformed via nanoindentation and subsequently lithiated. Polycrystalline films exhibited significantly higher hardness and elastic modulus than single crystals. Single crystals indented normal to the exposed (001) basal plane exhibited pronounced anisotropic deformation, including directional slip, crystallographically-guided cracking, anisotropic crack propagation, interlayer separation, and shear localization. Lithiation caused substantial softening, reduced hardness and modulus, and suppressed displacement bursts during nanoindentation, while previously indented regions showed crack formation and growth upon lithiation. Mechanical behavior of α-V 2 O 5 is strongly governed by crystallographic anisotropy and further altered by lithiation, with pre-existing deformation serving as a strong driver of fracture during ion insertion. These findings illuminate the coupling among ion insertion, deformation, and fracture in layered oxides and provide a basis for understanding and mitigating mechanical failure in electrochemical energy-storage materials.

Anisotropy↗
DOE OSTI2025

Anisotropic Thermal Transport in Quasi-2D Ruddlesden-Popper Hybrid Perovskite Superlattices

Two-dimensional hybrid metal-halide perovskites (2D-MHPs) have emerged as important solution-processed semiconductors with favorable optical and electronic properties for diverse applications in photovoltaics, optoelectronics, and spintronics. The quasi-2D layered structures, featuring large acoustic impedance mismatches between the organic and inorganic sublattices, are expected to result in distinct and anisotropic thermal transport properties along the cross-plane and in-plane directions. Here, we introduce transducer-free vibrational-pump-visible-probe (VPVP) approaches that enable accurate quantification of anisotropic thermal transport properties in various archetypical single-crystalline 2D-MHPs. Specifically, using VPVP spectroscopy and VPVP microscopy, we measure the anisotropic thermal diffusivities of 2D-MHPs with systematically varied Pb-I octahedral layer thicknesses, as well as organic spacer types and lengths, revealing how these structural parameters alter the cross-plane and in-plane thermal transport properties in distinct ways. While diffuse interface scattering plays an important role in dictating cross-plane thermal transport, in-plane thermal transport is primarily determined by phonon transport within interconnected inorganic layers. Density functional theory incorporating four-phonon scatterings provides further insight into the low thermal conductivity and modest thermal conduction anisotropy in 2D-MHPs. Our work demonstrates a new all-optical and noncontact method, which requires minimal sample preparation and allows direct visualization of cross-plane and in-plane thermal transport, potentially compatible with sample environments. The demonstrated VPVP approaches can advance understanding of thermal transport in 2D-MHPs as well as wide-ranging hybrid and polymeric semiconductors beyond 2D-MHPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion

We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗