Search NASA⌕ Search

SEARCH · Search NASA

Results for “Strongly inhomogeneous diffusivity”

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 OSTI2024

A finite difference informed random walker (FDiRW) solver for strongly inhomogeneous diffusion problems

In nature, many complex multi-physics coupling problems exhibit strong diffusivity inhomogeneity. For instance, in the context of radionuclide absorption by porous wasteform materials within a flowing waste stream, the difference of species’ diffusivity in solid and liquid phases spans by 3~8 orders of magnitude. To solve the diffusion equations with strongly inhomogeneous diffusivity, traditional discretization-based methods, such as the Finite Difference Method (FDM), require infinitesimally small time steps (<10 -10 ) as high spatial resolutions are employed in most microstructure evolution processes, leading to prohibitively high computational costs. Here, this work developed an integrated numerical approach (FDiRW: Finite Difference informed Random Walk) to tackle this challenge. The idea is that utilizing the Random Walk concept, the fast diffusion is modeled as a superposition of point source’s solution for a concentration distribution while FDM is used to obtain the point source’s solution at each node. A mesh-coarsening algorithm is developed to generate an exclusive coarse mesh for FDiRW approach to maximize its efficiency. The effectiveness of the coarse mesh-based FDiRW approach is validated by benchmarking Finite Difference solutions. Numerical results demonstrated that FDiRW achieves a remarkable 1000x computational efficiency improvement over FDM while preserving desired accuracy for a medium-sized model of 192 × 192 × 192 grids. Finally, as models scale up, a floating-point operations (PLOPs) analysis of the FDiRW algorithm reveals that its computational complexity grows quadratically in terms of the number of nodes employed in computation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Finite Difference informed Random Walk solver for simulating radiation defect evolution in polycrystalline structures with strongly inhomogeneous diffusivity

Diffusivity of species and defects on grain boundaries is usually several orders of magnitude larger than that inside grains. Such strongly inhomogeneous diffusivity requires prohibitively high computational demands for modeling microstructural evolution. Here, this paper presents a highly-efficient numerical solver, combining the Finite Difference method and Random Walk model, designed for accurately modeling strongly inhomogeneous diffusion within polycrystalline structures. The proposed solver, termed Finite Difference informed Random Walk (FDiRW), integrates a customized Finite Difference (cFD) scheme tailored for fast diffusion along thin grain boundaries represented by a single-layer of nodes. Numerical experiments demonstrate that the FDiRW solver achieves an impressive efficiency gain of 1560x compared to traditional Finite Difference methods while maintaining accuracy, making it feasible for personal computer machines to handle diffusional systems with strongly inhomogeneous diffusivity across static polycrystalline microstructures. The model has been successfully applied to simulate radiation defect evolution, showcasing its scalability to engineering scales in both length and time dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Efficient lasing in mixtures of helium and fluorine in diffuse discharges formed by runaway electrons

The parameters of stimulated lasing in diffuse discharges formed in mixtures of helium and fluorine in a strongly inhomogeneous electric field are investigated. Lasing is obtained in the visible and VUV spectral regions on the transitions of fluorine atoms and molecules. It is shown that lasing in He – F{sub 2} mixtures at a wavelength of 157 nm continues for several half-periods of the discharge current. Due to the homogeneity of the diffuse discharge, the maximum lasing efficiency of the F{sub 2} laser is 0.15 %, which corresponds to the efficiency of this type of lasers pumped by pre-ionised transverse volume discharges. (paper)

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Solute diffusion behavior during heat treatment and its impact in Sc-microalloyed Al-Cu system

Scandium (Sc) is a promising microalloying element that enhances the strength and thermal stability of aluminum (Al) alloys. However, these benefits are not fully realized in aluminum-copper (Al-Cu) systems, and corrosion resistance often declines due to the complex phase evolution and diffusion behavior of Cu and Sc. Here, to clarify this, solute diffusion behavior and Cu-Sc interaction during heat treatment (HT) were investigated using in situ synchrotron-based transmission x-ray microscopy (TXM), wide-angle x-ray scattering (WAXS), and x-ray absorption near-edge structure (XANES) spectroscopy. The results, supported by electron microscopy, reveal strong Cu-Sc bonding that significantly impedes solute diffusion, leading to inhomogeneous solute distribution and non-uniform precipitation. Moreover, the Al-Cu-Sc eutectic phase exhibits high thermal stability, resisting dissolution even near the matrix liquidus. These findings quantitatively elucidate the sluggish diffusion kinetics of Cu and Sc, which can help redesign HT schedules to improve both mechanical properties and corrosion resistance in Sc-microalloyed Al-Cu alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Highly skewed current–phase relation in superconductor–topological insulator–superconductor Josephson junctions

Three-dimensional topological insulators (TIs) in proximity with superconductors are expected to exhibit exotic phenomena, such as topological superconductivity (TSC) and Majorana-bound states (MBS), which may have applications in topological quantum computation. In superconductor–TI–superconductor Josephson junctions, the supercurrent versus the phase difference between the superconductors, referred to as the current–phase relation (CPR), reveals important information including the nature of the superconducting transport. Here, we study the induced superconductivity in gate-tunable Josephson junctions (JJs) made from topological insulator BiSbTeSe 2 with superconducting Nb electrodes. We observe highly skewed (non-sinusoidal) CPR in these junctions. The critical current, or the magnitude of the CPR, increases with decreasing temperature down to the lowest accessible temperature ( T ~ 20 mK), revealing the existence of low-energy modes in our junctions. The gate dependence shows that close to the Dirac point the CPR becomes less skewed, indicating the transport is more diffusive, most likely due to the presence of electron/hole puddles and charge inhomogeneity. Our experiments provide strong evidence that superconductivity is induced in the highly ballistic topological surface states (TSS) in our gate-tunable TI-based JJs. Furthermore, the measured CPR is in good agreement with the prediction of a model which calculates the phase-dependent eigenstate energies in our system, considering the finite width of the electrodes, as well as the TSS wave functions extending over the entire circumference of the TI.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Substrate-Dependent Exciton Diffusion and Annihilation in Chemically Treated MoS 2 and WS 2

Atomically thin semiconductors such as monolayer MoS 2 and WS 2 exhibit nonlinear exciton–exciton annihilation at notably low excitation densities (below ~10 excitons/μm 2 in exfoliated MoS 2 ). In this work, we show that the density threshold at which annihilation occurs can be tuned by changing the underlying substrate. When the supporting substrate is changed from SiO 2 to Al 2 O 3 or SrTiO 3 , the rate constant for second-order exciton–exciton annihilation, k XX [cm 2 /s], is reduced by 1 or 2 orders of magnitude, respectively. Using transient photoluminescence microscopy, we measure the effective room-temperature exciton diffusion coefficient in bis(trifluoromethane)sulfonimide-treated MoS 2 to be in the range D = 0.03–0.06 cm 2 /s, corresponding to a diffusion length of L D = 350 nm for an exciton lifetime of τ = 18 ns, which does not depend strongly on the substrate. We discuss possible mechanisms for the observed behavior, including substrate permittivity, long-range exciton–exciton or exciton–charge interactions, defect-mediated Auger recombination, and spatially inhomogeneous exciton populations arising from substrate-induced disorder. Exciton annihilation limits the overall efficiency of 2D semiconductor devices operating at high exciton densities; the ability to tune these interactions via the underlying substrate is an important step toward more efficient optoelectronic technologies featuring atomically thin materials.

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

Probing exciton diffusion dynamics in photosynthetic supercomplexes via exciton–exciton annihilation

Photosynthesis converts solar energy into chemical energy through coordinated energy transfer between light-harvesting complexes and reaction centers (RCs). Understanding exciton motion, particularly the exciton diffusion length, is essential for optimizing energy efficiency in photosystems. In this work, we combine intensity-cycling transient absorption spectroscopy with kinetic Monte Carlo (kMC) simulation to investigate exciton motion in the C2S2 photosystem II supercomplex of spinach. Using exciton–exciton annihilation, revealed in the fifth-order response, we experimentally estimate an exciton diffusion length of 10.9 nm based on a 3D normal diffusion model, suggesting the ability of excitons to traverse the supercomplex. However, kMC simulations reveal that exciton motion is sub-diffusive because of spatial constraints and the strong RC traps. An anomalous diffusion model analysis of the experimental data yields a diffusion length of 9.7 nm, while the simulated diffusion length is 7.4 nm. The variable exciton residence time across subunits, partly influenced by their connectivity to the trap, indicates inhomogeneous annihilation probability and suggests how plants balance efficient light harvesting with photoprotection. We also explore the influence of specific assumptions in the annihilation simulation, which are challenging to access in more complex environments, such as the thylakoid membrane. Our study provides a framework for studying exciton dynamics using exciton–exciton annihilation, which can be extended to understand the light-harvesting efficiencies of larger, more complex photosynthetic assemblies.

Zhang, Kunyan (ORCID:000000026830409X)↗
DOE OSTI2021

Enhanced superconductivity and ferroelectric quantum criticality in plastically deformed strontium titanate

The properties of quantum materials are commonly tuned using experimental variables such as pressure, magnetic field and doping. Here we explore a different approach using irreversible, plastic deformation of single crystals. We show that compressive plastic deformation induces low-dimensional superconductivity well above the superconducting transition temperature (Tc) of undeformed SrTiO3, with evidence of possible superconducting correlations at temperatures two orders of magnitude above the bulk Tc. The enhanced superconductivity is correlated with the appearance of self-organized dislocation structures, as revealed by diffuse neutron and X-ray scattering. We also observe deformation-induced signatures of quantum-critical ferroelectric fluctuations and inhomogeneous ferroelectric order using Raman scattering. Our results suggest that strain surrounding the self-organized dislocation structures induces local ferroelectricity and quantum-critical dynamics that strongly influence Tc, consistent with a theory of superconductivity enhanced by soft polar fluctuations. Our results demonstrate the potential of plastic deformation and dislocation engineering for the manipulation of electronic properties of quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Combined microscopy studies of complex electronic materials. Final report

The aim of this project was to study the electronic properties of exotic electronic materials which are only available as micron-scale samples or have microscopic inhomogeneity, using a special combination of multiple submicron scanning microscopy techniques allied to control of temperature, magnetic and electric fields and mechanical strain. In nanobeams of VO 2 , with its still mysterious metal-insulator phase transition, we determined the phase diagram and demonstrated that hydrogen diffusion is strongly channelled along one crystal axis. Using graphene, we developed and applied the scanning photo-Nernst effect, in which a current is generated when a laser is focused on the edge a sample in a magnetic field. We demonstrated that CrI 3 , a two-dimensional layered magnet, retains true Ising ferromagnetization down the monolayer level due to its out-of-plane anisotropy. In exfoliated WTe 2 , a layered topological semimetal, we discovered that monolayers are the first monolayer topological insulator and we imaged the helical edge states with scanning microwave microscopy. We also discovered that the same material when stacked in a few layers is a ferroelectric metal, and that the monolayer can be converted to a superconducting state by electron doping with a gate a temperatures below 1 K. Additionally, we developed the capability to perform momentum-resolved electronic spectroscopy on operating nanoscale devices by applying micron-scale angle resolved photoemission spectroscopy to two-dimensional heterostructure devices incorporating graphene and transition metal dichalcogenide semiconductors. Using it, we measured the band gaps of a series of monolayer semiconductors and observed band renormalization as a function of electron doping in monolayer WSe 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multifaceted Water Dynamics in Spherical Nanocages

We present here a new method to study position-dependent, anisotropic diffusion tensors inside spherically confined systems—a geometry that is common to many chemical nanoreactors. We use this method to elucidate the surprisingly rich solvent dynamics of confined water. The spatial variation of the strongly anisotropic diffusion predicted by the model agrees with the results of explicit molecular dynamics simulations. The same approach can be directly transferred to the transport of solutes to and from reaction sites located at nanoreactor interfaces. We complement our study by a detailed analysis of water hydrogen bond kinetics, which is intimately coupled to diffusion. Despite the inhomogeneity in structure and translational dynamics inside our nanocages, a single set of well-defined rate constants is sufficient to accurately describe the kinetics of hydrogen bond breaking and formation. We find that once system size effects have been eliminated, the residence times of water molecules inside the coordination shell of a hydrogen bond partner are well correlated to average diffusion constants obtained from the procedure above.

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

Magnetic Reconnection in Three Dimensions: Modeling and Analysis of Electromagnetic Drift Waves in the Adjacent Current Sheet

Magnetic reconnection at the subsolar magnetopause is persistently accompanied by strong fluctuations of the magnetic field (B), plasma density (n), and all components of the electric field (E) and current (J). The strongest fluctuations are at frequencies below the lower hybrid frequency and observed in a thin, intense current sheet adjacent to the electron diffusion region. In this current sheet, the background magnitudes of B and n are changing considerably, creating an inhomogeneous plasma environment. In this paper, we show that the fluctuations in B and n are consistent with an oscillatory displacement of the current sheet in the surface normal direction. The displacement is propagating parallel to the magnetic reconnection X-line. Wavelengths are on the order of or longer than the thickness of the current sheet to which they are confined, so we label these waves electromagnetic drift waves. E and J fluctuations are more complex than a simple displacement. They have significant variations in the component along B, which suggest that the drift waves also may be confined along B. The current sheet is supported by an electron drift driven by normal electric field, which, in turn, is balanced by an ion pressure gradient. We suggest that wave growth comes from an instability related to the drift between the electrons and ions. We discuss the possibility that drift waves can displace or penetrate into the electron diffusion region giving magnetic reconnection three-dimensional structure. Drift waves may corrugate the X-line, possibly breaking the X-line and generating turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sectoral and growth rate control on elemental uptake by individual calcite crystals

Heterogeneous distribution of trace elements (impurities) within individual calcite crystals is a phenomenon commonly observed in natural and laboratory systems. Changes in thermodynamic intensive parameters (mostly chemical potential and temperature) cannot always explain the inhomogeneous impurity patterns in calcite crystals and it has been suggested that growth rate and crystallographic orientation may exert strong effects on elemental incorporation into calcite. In addition, there are a number of experimental studies on micro-scale element (E) distribution between non-equivalent pairs of calcite vicinal faces (known as sectoral zoning); however, the variability of partition coefficients (e.g., K E = (E/Ca) calcite /(E/Ca) fluid ) within individual crystals remains undetermined. In this study, we have extended the work on elemental distribution between crystal sectors to evaluation of partition coefficients of trace and minor elements (Li, B, Mg, and Sr) in calcite crystal faces (10–14) and (01-12), whose growth rates were assessed. Growth entrapment model (GEM) and lattice strain theory were applied to explain K E heterogeneity by varying near-surface diffusivity of Mg and Sr and by varying surface enrichment factor for Li, B, Mg, and Sr. Decoupling of sectoral and growth rate effects reveals that sectoral zoning plays a key role in elemental distribution. More specifically, K Li and K B vary by more than one order of magnitude and K Mg varies by a factor of two within individual crystal faces. Strontium sectoral distribution is different from those of Li, B, and Mg and K Sr varies by up to a factor of two. Finally, these behaviors likely reflect different mechanisms of incorporation of Li, B, Mg, and Sr into calcite.

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

Magnetic anisotropy and two-dimensional short-range chemical ordering in Ba 1−𝑥 ⁢Na 𝑥 ⁢Fe 2 ⁢As 2

A true understanding of the properties of pnictide superconductors requires the development of high-quality materials and performing measurements designed to unravel their intrinsic properties and short-range nematic correlations which are often obscured by extrinsic effects such as poor crystallinity, inhomogeneity, domain formation, and twinning. Here, in this paper, we report the systematic growth of high-quality Na-substituted BaFe 2 As 2 single crystals and their characterization using pulsed magnetic fields x-ray diffraction and x-ray diffuse scattering. Analysis of the properties and compositions of the highest-quality crystals shows that their actual Na stoichiometry is about 50-60% of the nominal content and that the targeted production of crystals with specific compositions is accessible. We derived a reliable equation to estimate the Na stoichiometry based on the measured superconducting T c of these materials. Attempting to force spin reorientation and induce tetragonality, orthorhombic Ba 1-x Na x Fe 2 As 2 single crystals subjected to out-of-plane magnetic fields up to 31.4T are found to exhibit strong in-plane magnetic anisotropy demonstrated by the insufficiency of such high fields in manipulating the relative population of their twinned domains or in suppressing the orthorhombic order. Broad x-ray diffuse-intensity rods observed at temperatures between 30 and 300 K uncover short-range structural correlations. Local structure modeling together with 3D-Δ pair-distribution function mapping of real-space interatomic vectors show that the diffuse scattering arises from in-plane short-range chemical correlations of the Ba and Na atoms coupled with short-range atomic displacements within the same plane due to an effective size difference between the two atomic species.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A DECam view of the diffuse dwarf galaxy Crater II: the colour–magnitude diagram

ABSTRACT We present a deep Blanco/DECam colour–magnitude diagram (CMD) for the large but very diffuse Milky Way satellite dwarf galaxy Crater II. The CMD shows only old stars with a clearly bifurcated subgiant branch (SGB) that feeds a narrow red giant branch. The horizontal branch (HB) shows many RR Lyrae and red HB stars. Comparing the CMD with [Fe/H] = −2.0 and [α/Fe] = + 0.3 alpha-enhanced BaSTI isochrones indicates a mean age of 12.5 Gyr for the main event and a mean age of 10.5 Gyr for the brighter SGB. With such multiple star formation events Crater II shows similarity to more massive dwarfs that have intermediate age populations, however for Crater II there was early quenching of the star formation and no intermediate age or younger stars are present. The spatial distribution of Crater II stars overall is elliptical in the plane of the sky, the detailed distribution shows a lack of strong central concentration, and some inhomogeneities. The 10.5 Gyr subgiant and upper main-sequence stars show a slightly higher central concentration when compared to the 12.5 Gyr population. Matching to Gaia DR2 we find the proper motion of Crater II: μαcos δ = −0.14 ± 0.07, μδ = −0.10 ± 0.04 mas yr−1, approximately perpendicular to the semimajor axis of Crater II. Our results provide constraints on the star formation and chemical enrichment history of Crater II, but cannot definitively determine whether or not substantial mass has been lost over its lifetime.

Walker, A. R.↗
DOE OSTI2025

Cosmic-ray transport in inhomogeneous media

A theory of cosmic-ray transport in multiphase diffusive media is developed, with the specific application to cases in which the cosmic-ray diffusion coefficient has large spatial fluctuations that may be inherently multiscale. We demonstrate that the resulting transport of cosmic rays is diffusive in the long-time limit, with an average diffusion coefficient equal to the harmonic mean of the spatially varying diffusion coefficient. Thus, cosmic-ray transport is dominated by areas of low diffusion even if these areas occupy a relatively small, but not infinitesimal, fraction of the volume. On intermediate time-scales, the cosmic rays experience transient effective subdiffusion, as a result of low-diffusion regions interrupting long flights through high-diffusion regions. In the simplified case of a two-phase medium, we show that the extent and extremity of the subdiffusivity of cosmic-ray transport is controlled by the spectral exponent of the distribution of patch sizes of each of the phases. We finally show that, despite strongly influencing the confinement times, the multiphase medium is only capable of altering the energy dependence of cosmic-ray transport when there is a moderate (but not excessive) level of perpendicular diffusion across magnetic-field lines.

Astronomy and Astrophysics↗
DOE OSTI2021

Probing large-scale UV background inhomogeneity associated with quasars using metal absorption

ABSTRACT We study large-scale UV background inhomogeneity in three-dimensions associated with the observed quasar populations at high redshift. We do this by measuring metal absorption through quasar absorption spectrum stacking as a function distance to closest quasar in SDSS-IV/eBOSS on 10 s of comoving megaparsec scales. We study both intergalactic medium absorbers and mixed circumgalactic medium absorbers and probe absorption in O vi, C iv, and Si iv, and Si iii. Overall, stronger high ionization species absorption is seen closer to quasars at 2.4 < z < 3.1. O vi absorption shows a particularly strong change, with effects in C iv evident in some cases, and more marginal effects in Si iii and Si iv. We further study 2.05 < z < 2.4 (with weak signs of increasing homogeneity with time) and explore the study of metal absorption as a function of integrated SDSS-r band flux quasar flux (yielding consistent but less significant results). While the metal absorption does show sensitivity to large-scale 3D quasar proximity, the current incomplete quasar samples limit detailed interpretation. This work does, however, demonstrate that UV background inhomogeneities exist on scales of several 10 s of comoving megaparsecs associated with quasars and that they can be measured with precision by examining metal absorption in the intergalactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗