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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Overlapping Schwarz Methods Are Not Anisotropy‐Robust Multigrid Smoothers

We analyze overlapping multiplicative Schwarz methods as smoothers in the geometric multigrid solution of two-dimensional anisotropic diffusion problems. For diffusion equations, it is well known that the smoothing properties of point-wise smoothers, such as Gauss Seidel, rapidly deteriorate as the strength of anisotropy increases. On the other hand, global smoothers based on line smoothing are known to generally provide good smoothing for diffusion problems, independent of the anisotropy strength. Here, a natural question is whether global methods are really necessary to achieve good smoothing in such problems, or whether it can be obtained with locally overlapping block smoothers using sufficiently large blocks and overlap. Through local Fourier analysis and careful numerical experimentation, we show that global methods are indeed necessary to achieve anisotropy-robust smoothing. Specifically, for any fixed block size bounded sufficiently far away from the global domain size, we find that the smoothing properties of overlapping multiplicative Schwarz rapidly deteriorate with increasing anisotropy, irrespective of the amount of overlap between blocks. Moreover, our results indicate that anisotropy-robust smoothing requires blocks of diameter 𝒪⁡(𝜖 −1/2 ) for anisotropy ratio 𝜖 ∈(0,1] .

97 MATHEMATICS AND COMPUTING↗
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 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 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 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↗
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↗
DOE OSTI2026

Effect of magneto-mechanical synergism in the process-structure correlation in Fe–C alloys: A phase-field modeling approach

Applied magnetic fields can alter phase equilibria and kinetics in steels; however, quantitatively resolving how magnetic, chemical, and elastic driving forces jointly influence the microstructure remains challenging. We develop a quantitative magneto-mechanically coupled phase-field model for the Fe–C system that couples a CALPHAD-based chemical free energy with demagnetization-field magnetostatics and microelasticity. Here, the model reproduces single- and multi-particle evolution during the α → γ inverse transformation at 1023 K under external fields up to 20 T, including ellipsoidal morphologies observed experimentally at 8 T. Chemically driven growth is isotropic; a magnetic interaction introduces an anisotropic driving force that elongates γ precipitates along the field into ellipsoids, while elastic coherency promotes faceting, yielding elongated cuboidal or “brick-like” particles under combined magneto-elastic coupling. Growth kinetics increase with C content, and decrease with field strength and misfit strain. Multi-particle simulations reveal dipolar interaction-mediated coalescence for field-parallel neighbors and ripening for field-perpendicular neighbors. Incorporating field-dependent diffusivity from experiment slows kinetics as expected; a first-principles-motivated anisotropic diffusivity correction is estimated to be small (<2%). These results establish a process-structure link for magnetically assisted heat treatments of Fe–C alloys and provide guidance for microstructure control via chemo-magneto-mechanical synergism.

Magnetic field↗
DOE OSTI2025

Parallel diffusion operator for magnetized plasmas with improved spectral fidelity

Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.

Anisotropic diffusion↗
DOE OSTI2026

Origins of anisotropic linear magnetoresistance with isotropic mobility in Cd 3 As 2 films on GaAs$[110]$

Thin film synthesis allows for the potential to orient crystals in different orientations, permitting measurement of orientation-dependent material aspects such as band structures and transport anisotropy. Here, Dirac semimetal Cd 3 As 2 films are epitaxially grown on GaAs$[110]$ substrates, which has a $[001]$ orientation in-plane. Films contain domains of two different c-axis orientations resulting from an aligned a-axis in-plane. Magnetoresistance measurements performed along both $[1\bar{1}0]$ and $[001]$ substrate directions reveal similar mobility and carrier concentration, but much larger magnetoresistance along the $[001]$ direction, which can be explained by the guiding center diffusion model as arising from anisotropic disorder and different atomic spacings.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exploring nuclear structure with multiparticle azimuthal correlations at the LHC

Details of the nuclear structure of Xe, such as the quadrupole deformation and the nuclear diffuseness, are studied by extensive measurements of anisotropic-flow-related observables in Xe–Xe collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.44 TeV with the ALICE detector at the LHC. The results are compared with those from Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for a baseline, given that the 208 Pb nucleus exhibits a very weak deformation. Furthermore, comprehensive comparisons are performed with a state-of-the-art hybrid model using IP-Glasma+MUSIC+UrQMD. It is found that among various IP-Glasma+MUSIC+UrQMD calculations with different values of nuclear parameters, the one using a nuclear diffuseness parameter of α 0 = 0.492 and a nuclear quadrupole deformation parameter of β 2 = 0.207 provides a better description of the presented flow measurements. These studies represent the first systematic exploration of nuclear structure at TeV energies, utilizing a comprehensive set of anisotropic flow observables. The measurements serve as a critical experimental benchmark for rigorously testing the interplay between nuclear structure inputs and heavy-ion theoretical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

In-Plane Anisotropy of Charge Density Wave Fluctuations in 1⁢𝑇−TiSe 2

We report measurements of anisotropic triple-𝑞 charge density wave (CDW) fluctuations in the transition metal dichalcogenide 1⁢𝑇−TiSe 2 over a large volume of reciprocal space with x-ray diffuse scattering. Above the transition temperature, 𝑇 CDW , the in-plane diffuse scattering is marked by ellipses which reveal that the in-plane fluctuations are anisotropic. In addition, the out-of-plane diffuse scattering is characterized by rodlike structures which indicate that the CDW fluctuations in neighboring layers are largely decoupled. Furthermore, our analysis of the diffuse scattering line shapes and orientations suggests that the three charge density wave components contain independent phase fluctuations with a hierarchy of length scales, leading to intricate fluctuation patterns that go beyond the conventional 2D-to-3D crossover picture.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗