Search NASA⌕ Search

SEARCH · Search NASA

Results for “Anisotropic transport”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesoscopic finite-size effects of unconventional electron transport in PdCoO 2

A wide range of unconventional transport phenomena has recently been observed in single-crystal delafossite metals. Here, we present a theoretical framework to elucidate electron transport using a combination of first-principles calculations and numerical modeling of the anisotropic Boltzmann transport equation. Using PdCoO 2 as a model system, we study different microscopic electron and phonon scattering mechanisms and establish the mean free path hierarchy of quasiparticles at different temperatures. We treat the anisotropic Fermi surface explicitly to numerically obtain experimentally-accessible transport observables, which bridge between the “diffusive,” “ballistic,” and “hydrodynamic” transport regime limits. Further, we illustrate that the distinction between the “quasiballistic” and “quasihydrodynamic” regimes is challenging and often needs to be quantitative in nature. From first-principles calculations, we populate the resulting transport regime plots and demonstrate how the Fermi surface orientation adds complexity to the observed transport signatures in micrometer-scale devices. Our work provides key insights into microscopic interaction mechanisms on open hexagonal Fermi surfaces and establishes their connection to the macroscopic electron transport in finite-size channels.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Origin of Phase Separation in Ni-Rich Layered Oxide Cathode Materials During Electrochemical Cycling

In intercalation materials, the kinetics and uniformity of mass transport across the nanocrystalline domains dictate the structural reversibility and transport capability at the macroscopic level. (De)intercalation-induced interlayer disintegrations exhibit anisotropic crystallite size change. Due to the anisotropic mass transport mechanism, separated phases are inherently crystallographically oriented. One such material is LiNi $1–x–y$ Mn $x$ Co $y$ O 2 , which plays a pivotal role in advanced Li-ion batteries but suffers from severe phase inhomogeneities under fast charge or electrochemical aging. Here, using operando synchrotron techniques, we probe the mechanistic origins of the compositional and orientational-dependent phase separations during the electrochemical cycling of LiNi 0.8 Mn 0.1 Co 0.1 O 2 by comprehensive analysis of both in-plane and out-of-plane reflections. In the H2/H3 phase regime, in-plane domain propagation occurs due to increased covalency despite the severe decay of interlayer crystallographic order, resulting in the change of crystalline domain shape from 3D spheres to two-dimensional nanosheets. The crystallographically selective XRD line splitting is linked to the geometry of the facets as mass transfers along the ab plane. In conclusion, this work provides mechanistic insights into crystallographic orientation-dependent phase inhomogeneity under fast charge and extended cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Directed Self-Assembly of the Organic Semiconductor C8-BTBT-C8 in Anodic Aluminum Oxide Nanopores

Controlling the self-assembly of organic semiconductors at the nanoscale is critical for advancing high-performance electronic and photonic devices, yet it remains challenging due to their intrinsic anisotropic crystallization and sensitivity to processing conditions. Here, in this study, we demonstrate that cylindrical nanoconfinement within anodic aluminum oxide membranes provides a versatile platform to precisely tune the molecular orientation and phase behavior of the prototypical organic semiconductor 2,7-dioctyl[1]benzothieno[3,2-b][1]benzothiophene (C8-BTBT-C8). Combining temperature-dependent high-resolution synchrotron X-ray scattering with optical birefringence measurements, we uncover that confinement geometries (pore diameters 25–180 nm) and surface chemistry govern the emergence of distinct smectic A textures, featuring molecular layers either parallel or perpendicular to the pore axis. The competition between axial and radial smectic layering is modulated by pore size, surface hydrophilicity, and thermal history, enabling reversible control over domain orientations and transitions between liquid crystalline and crystalline states. Notably, nanoconfinement stabilizes the smectic phase over an expanded temperature range compared to bulk, while inducing complex multidomain configurations owing to geometric constraints and anchoring conditions. Our results elucidate fundamental mechanisms by which anisotropic nanoscale confinement directs the self-organization of highly conjugated organic molecules, with implications for optimizing directional charge transport and anisotropic optical responses in organic–inorganic hybrid nanoarchitectures. This study establishes nanoconfinement as a powerful strategy to engineer morphology and functional properties in organic semiconducting materials with nanoscale precision.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Anisotropic Thermal Conductivity of Inkjet-Printed 2D Crystal Films: Role of the Microstructure and Interfaces

Two-dimensional (2D) materials are uniquely suited for highly anisotropic thermal transport, which is important in thermoelectrics, thermal barrier coatings, and heat spreaders. Solution-processed 2D materials are attractive for simple, low-cost, and large-scale fabrication of devices on, virtually, any substrate. However, to date, there are only few reports with contrasting results on the thermal conductivity of graphene films, while thermal transport has been hardly measured for other types of solution-processed 2D material films. In this work, inkjet-printed graphene, h-BN and MoS2 films are demonstrated with thermal conductivities of ∼10 Wm−1K−1 and ∼0.3 Wm−1K−1 along and across the basal plane, respectively, giving rise to an anisotropy of ∼30, hardly dependent on the material type and annealing treatment. First-principles calculations indicate that portion of the phonon spectrum is cut-off by the quality of the thermal contact for transport along the plane, yet the ultra-low conductivity across the plane is associated with high-transmissivity interfaces. These findings can drive the design of highly anisotropic 2D material films for heat management applications.

2D materials↗
DOE OSTI2020

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗
DOE OSTI2022

Complex structure due to As bonding and interplay with electronic structure in superconducting BaNi 2 As 2

BaNi 2 As 2 is a superconductor chemically related to the Fe-based superconductors, with a complex and poorly understood structural phase transition. We show, based on first-principles calculations, that, in fact, there are two distinct competing structures. These structures are different from the electronic, transport, and bonding points of view but are close in energy. These arise due to complex As bonding patterns and drive distortions of the Ni layers. This is supported by photoemission experiments. This is very distinct from views of the distortion as being primarily driven by electrons at the Fermi surface as in a classical charge density wave, or by correlated electron physics associated with the Ni d electrons and the Ni d orbitals. The structural distortion, although connected with As does lead to an interplay of electronic and structural behavior including induced anisotropic electronic transport. Here, the local bonding nature of the instabilities and the competition between distortions is discussed in the context of the complex behavior observed in BaNi 2 As 2 samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Aligning thermal and current quenches with a high density low-Z injection

The conventional approach for thermal quench (TQ) mitigation in a tokamak disruption is through a high- Z impurity injection that radiates away the plasma’s thermal energy before it reaches the wall. The downside is a robust Ohmic-to-runaway current conversion due to the radiatively clamped low post-thermal-quench electron temperature. An alternative approach is to deploy a low- Z (either deuterium or hydrogen) injection that aims to slow down the TQ, and ideally aligns it with the current quench (CQ). This approach has been investigated here via 3D MHD simulations using the PIXIE3D code. By boosting the hydrogen density, a fusion-grade plasma is dilutionally cooled at approximately the original pressure. Energy loss to the wall is controlled by a Bohm outflow condition at the boundary where the magnetic field intercepts a thin plasma sheath at the wall, in addition to Bremsstrahlung bulk losses. Robust MHD instabilities proceed as usual, while the collisionality of the plasma has been greatly increased and parallel transport is now in the Braginskii regime. The main conclusion of this study is that the decreased transport loss along open field lines due to a sufficient low- Z injection slows down the TQ rate to the order of 20 ms, aligned with the CQ timescale for a 15 MA ITER plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetic field dependence of low-energy magnons, anisotropic heat conduction, and spontaneous relaxation of magnetic domains in the cubic helimagnet ZnCr 2 Se 4

Here, anisotropic low-temperature properties of the cubic spinel helimagnet ZnCr 2 Se 4 in the single-domain spin-spiral state are investigated by a combination of neutron scattering, thermal conductivity, ultrasound velocity, and dilatometry measurements. In an applied magnetic field, neutron spectroscopy shows a complex and nonmonotonic evolution of the spin-wave spectrum across the quantum-critical point that separates the spin-spiral phase from the field-polarized ferromagnetic phase at high fields. A tiny spin gap of the pseudo-Goldstone magnon mode, observed at wave vectors that are structurally equivalent but orthogonal to the propagation vector of the spin helix, vanishes at this quantum critical point, restoring the cubic symmetry in the magnetic subsystem. The anisotropy imposed by the spin helix has only a minor influence on the lattice structure and sound velocity but has a much stronger effect on the heat conductivities measured parallel and perpendicular to the magnetic propagation vector. The thermal transport is anisotropic at T≲ 2K, highly sensitive to an external magnetic field, and likely results directly from magnonic heat conduction. We also report long-time thermal relaxation phenomena, revealed by capacitive dilatometry, which are due to magnetic domain motion related to the destruction of the single-domain magnetic state, initially stabilized in the sample by the application and removal of magnetic field. Our results can be generalized to a broad class of helimagnetic materials in which a discrete lattice symmetry is spontaneously broken by the magnetic order.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Light–ferroelectric interaction in two-dimensional lead iodide perovskites

The unique physical properties of two-dimensional (2D) metal halide perovskites (MHPs) such as nonlinear optics, anisotropic charge transport, and ferroelectricity have made these materials promising candidates for multifunctional applications. Recently, fluorine derivatives such as 4,4-difluoropiperidinium lead iodide perovskite or (4,4-DFPD, C 5 H 10 F 2 N) 2 PbI 4 have shown strong ferroelectricity as compared to other 2D MHPs. Although it was previously addressed that the ferroelectricity in MHPs can be affected by illumination, the underlying physical mechanisms of light–ferroelectricity interaction in 2D MHPs are still lacking. Here, we explore the electromechanical responses in 4,4-(DFPD) 2 PbI 4 thin films using advanced scanning probe microscopy techniques revealing ferroelectric domain structures. In this work, hysteretic ferroelectric loops measured by contact-Kelvin probe force microscopy are dependent on domain structures under dark conditions, while ferroelectricity weakens under illumination. The X-ray diffraction patterns exhibit significant changes in preferential orientation of individual lattice planes under illumination. Particularly, the reduced intensity of the (1 1 1) lattice plane under illumination leads to transitioning from a ferroelectric to a paraelectric phase. The instability of positive ions, especially molecular organic cations, is observed under illumination by time-of-flight secondary ion mass spectrometry. The combination of crystallographic orientation and chemical changes under illumination clearly contributes to the origin of light–ferroelectricity interaction in 2D (4,4-DFPD, C 5 H 10 F 2 N) 2 PbI 4 .

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

Structural and physical properties of two distinct 2D lead halides with intercalated Cu( II )

Transition metal cation intercalation between the layers of two-dimensional (2D) metal halides is an underexplored research area. In this work we focus on the synthesis and physical property characterizations of two layered hybrid lead halides: a new compound [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 and the previously reported [Cu(O 2 C–(CH 2 ) 3 –NH 3 ) 2 ]PbBr 4 . These compounds exhibit 2D layered crystal structures with incorporated Cu 2+ between the metal halide layers, which is achieved by combining Cu(II) and lead bromide with suitable amino acid precursors. The resultant [Cu(O 2 C–(CH 2 ) 3 –NH 3 ) 2 ]PbBr 4 adopts a 2D layered perovskite structure, whereas the new compound [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 crystallizes with a new structure type based on edge-sharing dodecahedral PbBr5O3 building blocks. [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 is a semiconductor with a bandgap of 3.25 eV. It shows anisotropic charge transport properties with a semiconductor resistivity of 1.44 × 10 10 Ω cm (measured along the a-axis) and 2.17 × 10 10 Ω cm (along the bc-plane), respectively. The fabricated prototype detector based on this material showed response to soft low-energy X-rays at 8 keV with a detector sensitivity of 1462.7 μCGy –1 cm –2 , indicating its potential application for ionizing radiation detection. Finally, these encouraging results are discussed together with the results from density functional theory calculations, optical, magnetic, and thermal property characterization experiments.

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

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗
DOE OSTI2021

Modification of favorable average curvature effect by changing parallel sound wave behavior in tokamak plasmas

Abstract The favorable average curvature effect, also known as the GGJ effect (Glasser et al 1975 Phys. Fluids 18 875), is intrinsically associated with parallel sound wave propagation in a tokamak plasma. This work investigates how the GGJ effect is modified by changing the parallel sound wave behavior. Two physics models beyond the standard single fluid theory, i.e. an anisotropic thermal transport model and a parallel sound wave damping model, are employed to change parallel sound waves in a toroidal plasma, and the consequence on the GGJ effect is demonstrated for two important classes of problems, i.e. the resistive plasma response to the applied resonant magnetic perturbation and the stability of the tearing mode. Toroidal modeling reveals that the GGJ effect is significantly altered by both of the aforementioned physics effects. Compared to the thermal transport physics, which completely removes the GGJ effect, the sound wave damping effect only offers partial mitigation. The differences between these two models are further illustrated in terms of the radial structure of the shielding current and the eigenfunction of the tearing instability. In particular, a fundamental reason for complete suppression of the GGJ effect by the thermal transport is identified as an extra toroidal coupling of the poloidal harmonics.

Physics↗
DOE OSTI2022

Toroidal modeling of energetic passing particle drift kinetic effects on tearing mode stability

Abstract Drift kinetic effects of the neutral beam injection induced passing energetic particles (EPs) on the linear stability of the n = 1 tearing mode (TM) (with the dominant poloidal harmonic of m = 2) are numerically investigated utilizing the MARS-K code (Liu et al 2008 Phys. Plasmas 15 112503), in a tokamak plasma with finite equilibrium pressure and anisotropic thermal transport. In the low plasma pressure regime, it is found that co- (counter-) passing EPs stabilize (destabilize) the TM, agreeing with previous studies. However, as the plasma pressure increases beyond a critical value, it is found that co-passing EPs also destabilize the mode. An in-depth analysis reveals that the net effect of co-passing EPs is a result of competition between the stabilizing contribution from the non-adiabatic drift kinetic terms and the destabilizing contribution associated with adiabatic terms, with the latter becoming more dominant at higher equilibrium pressure. Non-perturbative magnetohydrodynamic-kinetic hybrid modeling also finds that co- and counter-passing EPs modify the TM eigenfunction differently, with the counter-passing EPs enhancing the sideband harmonics. Furthermore, effects of the plasma resistivity and toroidal rotation, as well as that of the equilibrium distribution of EPs in the particle pitch angle space, are also investigated, showing asymmetric results on the TM stability between the co- and counter-passing EPs. The first order finite orbit width correction is found to be stabilizing with co-passing EPs and destabilizing with counter-passing particles. Finally, drift resonances between passing EPs and the TM induce finite frequency to the mode and generate finite net torques inside the plasma, due to the neoclassical toroidal viscosity and the Reynolds stress associated with 3D perturbations.

Physics↗
DOE OSTI2023

Thermokinetics of point defects in α-Fe 2 O 3

Abstract Point defect formation and migration in oxides governs a wide range of phenomena from corrosion kinetics and radiation damage evolution to electronic properties. In this study, we examine the thermodynamics and kinetics of anion and cation point defects using density functional theory in hematite ( α -Fe 2 O 3 ), an important iron oxide highly relevant in both corrosion of steels and water-splitting applications. These calculations indicate that the migration barriers for point defects can vary significantly with charge state, particularly for cation interstitials. Additionally, we find multiple possible migration pathways for many of the point defects in this material, related to the low symmetry of the corundum crystal structure. The possible percolation paths are examined, using the barriers to determine the magnitude and anisotropy of long-range diffusion. Our findings suggest highly anisotropic mass transport in hematite, favoring diffusion along the c -axis of the crystal. In addition, we have considered the point defect formation energetics using the largest Fe 2 O 3 supercell reported to date.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the effect of scalar flux weighing of linearly anisotropic scattering matrices in few-group transport calculations

The majority of the codes available for homogenized group constant generation for deterministic transport calculations apply the approximation of scalar flux weighting during energy group condensation of higher-order anisotropic scattering matrices. In this paper, we discuss the effect of scalar flux weighting of the linearly anisotropic scattering matrices in the frame of S P{sub 3} and S{sub 12} calculations performed for a two-dimensional VVER-440 reactor benchmark. To compare group constants generated for two neutron energy groups, an infinite pin cell was homogenized with Serpent 2 and ERANOS ECCO. Serpent 2 applies scalar flux, while ERANOS ECCO performs current weighting of the linearly anisotropic scattering matrices during energy group condensation. For analyzing the effect of the various weighting options, three simple reactor models were built assuming different core sizes using standard rectangular assemblies with 15*15 fuel pins. Diffusion, SP{sub 3} and S{sub 12} calculations were performed for the 3 models using group constants generated with Serpent 2 and ERANOS ECCO. The effect of scalar flux weighting of linearly anisotropic scattering matrices is shown by comparing the decrease in reactivity due to the decreased reactor size, as well as the assembly power distribution to reference results obtained with Serpent 2 Monte Carlo calculations. Neglecting higher than linearly anisotropic scattering and indirect application of diffusion coefficients in higher-order transport calculations is advised if angular flux-moment spectra weighted higher-order scattering matrices cannot be generated. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Analysis of Anisotropic Nonlocal Diffusion Models: Well-posedness of Fractional Problems for Anomalous Transport

We analyze the well-posedness of an anisotropic, nonlocal diffusion equation. Establishing an equivalence between weighted and unweighted anisotropic nonlocal diffusion operators in the vein of unified nonlocal vector calculus, we apply our analysis to a class of fractional-order operators and present rigorous estimates for the solution of the corresponding anisotropic anomalous diffusion equation. Furthermore, we extend our analysis to the anisotropic diffusion-advection equation and prove well-posedness for fractional orders s ∊ [0.5, 1). We also present an application of the advection-diffusion equation to anomalous transport of solutes.

97 MATHEMATICS AND COMPUTING↗