Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective 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 253 records · Page 14

Effects of irradiation temperature on the microstructure and deuterium retention in γ-LiAlO 2 pellets

This report presents the experimental findings obtained from November 2023 to September 2024. The study aims to investigate the effects of temperature on the microstructure, deuterium (D) retention, and lithium (Li) loss in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation to a high dose. Sequential ion irradiation has been a key method in our previous studies to simulate the behavior of γ-LiAlO 2 under neutron irradiation. In addition to dose and dose rate, irradiation temperature is a critical factor influencing microstructural and compositional changes. LiAl 5 O 8 precipitates have been observed in γ-LiAlO 2 pellets irradiated with reactor neutrons at ~300°C. These precipitates also form during ion irradiation at an elevated temperature of 500°C, but not at 300°C. In our ion irradiation experiments, the dose rate is typically three orders of magnitude higher than that of neutron irradiation, leading to a more rapid damage production. To better emulate the microstructural features in neutron-irradiated pellets using ion irradiation, a higher irradiation temperature is needed to accelerate the diffusion of point defects and enhance defect recovery rates, thereby compensating for the effects of the higher dose rate. The microstructural changes observed are the result of competing processes occurring during ion irradiation.

36 MATERIALS SCIENCE↗

NOCI-F Electronic Couplings in Assemblies of Indolonaphthyridine Molecules: From Dimers to the Full Stack

Key electronic processes related to molecular excitonic states of finite stacks of indolonaphthyridine molecules are analyzed via the non-orthogonal configuration interaction with fragments (NOCI-F) method. Indolonaphthyridine is an organic chromophore that can undergo several electronic photoexcitation-related intermolecular processes, such as exciton and electron transfer. The structures studied here are noncrystalline arrangements built as either ordered stacks of indolonaphthyridine or stacks extracted from molecular dynamics simulations including thermal disorder. Taking dimers or trimers from either model, we performed CASSCF and NOCI-F calculations to quantify the intermolecular electronic couplings governing singlet fission, excited singlet and triplet diffusion, and hole and electron diffusion processes. Also, comparing the results for the different models, we studied the effect of structural disorder and distortion on these couplings. Finally, we present a newly developed, advanced postanalysis tool. It takes the NOCI-F data as input to carry out a multifragment full Hamiltonian procedure that involves the complete stack, providing physical information not available from the dimer/trimer models, hence giving access to additional insight into the material’s properties.

coupling reactions↗

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗

In-situ two-dimensional temperature measurements using x-ray fluorescence spectroscopy in laminar flames with high silica particle concentrations

X-ray fluorescence spectroscopy (XRF) was used to measure temperatures and study mixing, for the first time in silica particle synthesis flames. Hexamethyldisiloxane (HMDSO) and trimethylsilanol (TMSO) were the particle precursors. A multi-element diffusion burner was used to produce a flat methane flame, and the precursors, dilute in inert gas, were injected via a central jet. Krypton was the fluorescent medium at 3.2 % concentration by volume. Scans with Kr in the central flow and not in the main flow were made to assess the mixing effects between the central and main gas flows. When HMDSO or TMSO were added, a secondary diffusion flame formed between the jet and the main methane flame. The results revealed a dramatic change in the centerline temperature profile of the jet gases when HMDSO or TMSO were added. The main methane flame stoichiometry also affected the temperature profiles. The results show HMDSO and TMSO reactions are initiated in a lowtemperature and low-oxygen concentration region of the jet where thermal decomposition is not expected to be significant. Reaction of the particle precursors is therefore attributed to radical transport from the main methane flame. In the current work, particle number densities of up to 270 g/m 3 locally are estimated. Thus, the study also demonstrates the capability of the XRF technique for high spatial fidelity measurements in flames with high concentrations of condensed-phase particles, leveraging the attribute that the XRF signal is generally not impacted by condensed-phase interferences. In conclusion, the observations and data obtained in this study inform likely reaction pathways for this important class of siloxane compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Charge Carrier and Spin Diffusion in a Polycrystalline and Single Crystal Lead Halide Perovskite Semiconductor

The directed movement of spin-bearing excitations in semiconductors can enable many potential schemes for spintronic applications. Understanding the mechanism of spin motion, as opposed to charge carrier or exciton motion, may require specialized techniques and sample conditions. Here we employ the noncontact and time-resolved technique of light-induced transient grating spectroscopy (LITG) to measure both carrier and spin motion in a perovskite semiconductor with varying crystallinity. The carrier motion aligns with expectations from past studies undertaken at low fluence, revealing an ambipolar diffusion coefficient on the order of 1 cm 2 /sec and diffusion length of roughly 1.5 μm for single crystals that is strongly attenuated as crystallite size decreases. The spin diffusion is measured with cross-polarized excitation beams and uncovers an intensity-dependent coefficient rising above 100 cm 2 /sec. The fast room-temperature spin relaxation limits the spin diffusion length, but the remarkable speedup of spin over carrier diffusion suggests a mechanism involving a combination of exchange-mediated and doping effects that enhance spin transport.

14 SOLAR ENERGY↗

Physically constrained 3D diffusion for inverse design of fiber-reinforced polymer composite materials

Designing fiber-reinforced polymer composites (FRPCs) with a tailored nonlinear stress-strain response is crucial for applications such as energy absorption in crash structures, flexible robotics, and impact-resistant protective gear. However, the inherent complexities of composite materials and the multitude of parameters involved, render traditional design and optimization methods inadequate for achieving effective inverse design of composites. In this paper, we present an AI-based inverse design framework that effectively and efficiently generates FRPCs with targeted nonlinear stress-strain responses. We introduce a physically constrained diffusion model (PC3D_Diffusion) capable of managing the complexities of composite materials and producing detailed, high-quality designs. We propose a loss-guided, learning-free approach to generate physically feasible microstructure designs by explicitly enforcing physical constraints during the generation process. For training purposes, 1.35 million FRPC samples were created, and their corresponding stress-strain curves were computed using established physics-based computational models. The results show that PC3D_Diffusion consistently generates high-quality designs with tailored mechanical behaviors, while guaranteeing compliance with the physical constraints. PC3D_Diffusion advances FRPC inverse design and may facilitate the inverse design of other 3D materials, offering potential applications in industries reliant on materials with custom mechanical properties.

Xu, Pei [Clemson Univ., SC (United States)]↗

Rooftop Photovoltaics and Electric Vehicle Co-Adoption: Attitudes, Norms, Diffusion, and Economics

The overall objectives of the project were to (1) Identify personal and social norms, peer effects, demographics and contextual factors (such as household energy expenditures, experience with outages, charging infrastructure availability) that may facilitate or hinder Solar-EV co-adoption; (2) outline strategies for the design of programs and practices that will enable reasoned RPV and EV adoption and their co-adoption; and (3) assemble a dataset comprised of survey responses of RPV-EV co-adopters, RPV-only adopters, EV-only adopters, and non-adopters. These objectives were only partially fulfilled given the decision to terminate the project. The research team was able to complete the semi-structured interviews, the core BY1 task and summarize the results below. As well, once the team learned of the research project’s fate, it sought and obtained the DOE SETO Technical Manager’s approval to repurpose funds for BY2 tasks to the extent feasible. Rather than invest further resources in expanding the pool of potential survey respondents and having no budget to survey, the research team repurposed funds subject and undertook a more limited (geographically) survey, but also expanded to it to include contingent valuation questions. The analysis of the survey data will necessarily be delayed given the lack of budget support, although some information on the respondent sample is in the results. Given a lack of attention to co-adoption previously, the insights generated in this project should assist decision-makers and industry actors in the development of programs and practices that enable reasoned RPV adoption and EV adoption and their co-adoption.

14 SOLAR ENERGY↗

Anharmonic phonons, thermal expansion, and nuclear quantum effects in Zn

The anharmonic behavior of phonons and thermal expansion of hexagonal zinc were studied from 15 to 690 K by inelastic neutron scattering (INS) and ab initio simulations. Phonon spectra were measured for Q-points covering the full Brillouin zone, giving the phonon density of states (DOS) and dispersions along high-symmetry directions. The dispersions were sharp at 15 K, but diffuse intensity was observed at energies above them. The dispersions broadened with temperature T and the diffuse intensity became stronger. This diffuse intensity appeared in all INS measurements and simulations, except for classical molecular dynamics at 15 K. The temperature-dependent effective potential (TDEP) method, which included the nuclear quantum effect from zero-point vibrational dynamics, was used to calculate the free energy and thermal expansion. For T < 100 K nuclear quantum effects were important for obtaining the correct negative thermal expansion, and path integral molecular dynamics (PIMD) was particularly effective for obtaining the negative thermal expansion in the basal plane. A Heisenberg-Langevin model for interacting phonons coupled to a thermal bath was able to reproduce the shape and intensity of the diffuse spectral features.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Superionic Surface Li-Ion Transport in Carbonaceous Materials

Unlike Li-ion transport in the bulk of carbonaceous materials, little is known about Li-ion diffusion on their surface. Here, in this study, we have discovered an ultrafast Li-ion transport phenomenon on the surface of carbonaceous materials with limited reversible Li insertion capacity and high surface area. An ionic conductivity of 18.1 mS cm –1 at room temperature is observed in lithiated Ketjen black (KB), far exceeding those of most solid-state ion conductors. Theoretical calculations reveal low diffusion barriers for the surface Li species. As a result, lithiated KB functions effectively as an interlayer between Li and solid-state electrolytes (SSEs) to mitigate dendrite growth. Further, lithiated KB acts as a high-performance mixed ionic–electronic conductor and replaces solid electrolytes to enhance graphite anode performance, demonstrating full utilization with ∼85% capacity retention over 300 cycles. The discovery of this surface-mediated ultrafast Li-ion transport mechanism provides new directions for the design of solid-state ion conductors and solid-state batteries.

Li metal batteries↗

The Impact of the 8–10 March 2012 Geomagnetic Storm on Inner Zone Protons as Measured by Van Allen Probes

The Relativistic Electron Proton Telescope (REPT) instrument on the Van Allen Probes observed a double-peaked inner zone proton population throughout the 7 year lifetime of the mission. M. Hudson et al. (2023) showed that a strong SEP event accompanied by a CME-shock in early March 2012 provided the Solar Energetic Proton (SEP) source for the higher L trapped proton population, which then diffused radially inward to be observed by REPT at L = 2. The study followed trajectories of SEP protons launched isotropically from a sphere at 7 Re for 2.5 hr in fields calculated by the LFM-RCM global MHD model, which includes electric fields needed to model the transport and trapping of the protons by the shock, and then a radial diffusion simulation was run for 2 years using the result from the test-particle simulation as the initial condition. The simulation result was compared with REPT measurement in November 2013 and showed reasonable agreement. However, the simulation overestimated the Phase Space Density by a factor of four due to lack of field line curvature scattering during the storm in the model. In this study, a test-particle simulation is performed for 2 days following the injection and trapping of protons in March 2012 using TS05 fields to simulate the field line curvature scattering of the trapped SEP due to the buildup of the ring current during the geomagnetic storm. The resulting sample distribution was then weighted using the flux at the end of the two-hour MHD-test particle simulation. A radial diffusion simulation is then run using the initial profile that included the loss effect, with improved comparison with REPT measurements after 2 years.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing the influence of ion-pairing on ligand-field excited-state dynamics

Exploration of the photophysical and photochemical properties of transition metal complexes has driven ground-breaking advancements in solar energy conversion technologies, including photoredox catalysis. While significant research has been devoted to understanding excited state properties of second- and third-row transition metal complexes, earth-abundant first-row metal complexes have received comparatively little attention in this context until very recently. In particular, the role of ion-pairing – which has been identified as a potentially significant factor for Ir(III)-based photosensitizers – has not been examined with regard to its influence on the ligand-field excited states that dominate much of first-row photophysics. A key challenge in studying ion-pair interactions lies in quantifying the extent and nature of ion-pairing, particularly in non-aqueous media where the vast majority of photophysical studies are performed. Cobalt(III) polypyridyl complexes provide an attractive platform to address such questions due to their demonstrated potential for applications in photoredox catalysis involving ligand-field excited states. In the present study, we prepared a cobalt(III) polypyridyl complex, [Co(4,4′-OMebpy) 3 ](BAr F 4 ) 3 (where 4,4′-OMebpy is 4,4′-dimethoxy-2,2′-bipyridine and BAr F 4 is tetrakis(3,5-bis(trifluoromethyl)-phenyl)borate) to probe ion-pairing in non-aqueous solutions. Specifically, analysis of data acquired from both variable-temperature diffusion ordered spectroscopy (DOSY) NMR and 1-D rotating-frame nuclear Overhauser effect (ROE) experiments allowed us to identify and differentiate between solvent-separated ion pairs in high-dielectric media and contact ion pairs in a low-dielectric solvent. Time-resolved absorption spectroscopy was then used to measure ground-state recovery dynamics under these varying conditions of ion-pairing, the results of which revealed an increase in excited-state lifetime for contact ion-pairs that we suggest arises from a reduction in outer-sphere reorganization energy relative to conditions which favored solvent-separated ion pairs. We believe this study demonstrates that one can leverage broadly available NMR-based methods to understand ion-pairing in non-aqueous solutions, which in turn can provide a microscopic picture of intermolecular interactions that can impact the photophysical properties of transition metal-based chromophores.

Ghosh, Atanu [Michigan State University, East Lans↗

Nanoscale structural correlations in a model cuprate superconductor

Understanding the extent and role of inhomogeneity is a pivotal challenge in the physics of cuprate superconductors. While it is known that structural and electronic inhomogeneity is prevalent in the cuprates, it has proven difficult to disentangle compound-specific features from universally relevant effects. Here, in this study, we combine advanced neutron and x-ray diffuse scattering with numerical modeling to obtain insight into bulk structural correlations in HgBa 2 ⁢ CuO 4+δ . This cuprate exhibits a high optimal transition temperature of nearly 100 K, pristine charge-transport behavior, and a simple average crystal structure without long-range structural instabilities, and is therefore uniquely suited for investigations of intrinsic inhomogeneity. We uncover diffuse reciprocal-space patterns that correspond to prominent nanoscale correlations of atomic displacements perpendicular to the CuO 2 planes. The real-space nature of the correlations is revealed through three-dimensional pair distribution function analysis and complementary numerical refinement. We find that relative displacements of ionic and CuO 2 layers play a crucial role, and that the structural inhomogeneity is not directly caused by the presence of conventional point defects. The observed correlations are therefore intrinsic to HgBa 2 ⁢ CuO 4+δ , and thus likely important for the physics of cuprates more broadly. It is possible that the structural correlations are closely related to the unusual superconducting fluctuations and Mott-localization in these complex oxides. As advances in scattering techniques yield increasingly comprehensive data, the experimental and analysis tools developed here for large volumes of diffuse scattering data can be expected to aid future investigations of a wide range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modification of Ni-20Cr corrosion dealloying behavior in molten fluorides via cold work induced plastic deformation

The corrosion dealloying behavior of cold-worked (CW) Ni20Cr alloy (wt%) was studied in molten LiF-NaF-KF (or FLiNaK) salts at 600 °C, equal to a homologous temperature (TH) of 0.52. Alloys were cold-rolled to achieve reductions of thickness of 10%, 30%, and 50% introducing plastic deformation and a high density of dislocations. Potentiostatic holds (Eapplied) were applied in two different electrode potential regimes. At 1.75VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, Cr dealloying to Cr(II) and Cr(III) is predominant, while at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, both Ni and Cr are oxidized in molten FLiNaK at 600 °C. In these potential regimes, dealloyed Ni20Cr displayed bicontinuous porosity within the grain interior and at grain boundaries, driven by the high driving force for Cr dissolution and sustained by defect mediated outward solid state diffusion of Cr in parallel with surface diffusion of Ni. The bicontinuous porous structure developed was observed to undergo further coarsening and densification of the Ni-rich ligaments at a higher electrode potential. The main effect of CW observed is the introduction of plastic deformation and dislocation substructures that serve as short-circuit paths for Cr solid state diffusion to surfaces exposed to FLiNaK. This modified the evolution of the bicontinuous porous structure which increased with CW substantially. Kinetic analysis reveals that the Cr dealloying at +1.75 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ and 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$ is initially charge transfer controlled, except for the 50% CW condition at 1.90 VK+/K$${\text{V}}_{{\text{K}}^{+}/{\text{K}}}$$, where the process becomes limited by slow Cr defect mediated bulk diffusion. The rate determining factors are explored and compared to experimental results.

Chan, Ho Lun↗

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗

Topotactic Phase Transformation of Lithiated Spinel to Layered LiMn0.5Ni0.5O2: The Interaction of 3-D and 2-D Li-ion Diffusion

This study investigates the structural evolution of LiMn0.5Ni0.5O2 cathode materials for Li-ion batteries as a function of synthesis temperature and its effect on electrochemical performance. It is demonstrated that, as the synthesis temperature increases from 400 to 900 ?C, a gradual topotactic transformation occurs between a lithiated spinel structure, denoted herein as “lithium-excess spinel” LxS-LiMn0.5Ni0.5O2 (or LxS-LMNO), and the well-known layered LiMn0.5Ni0.5O2 structure prepared at high temperature, HT-LiMn0.5Ni0.5O2 (HT-LMNO). The electrochemical capacity of the LiMn0.5Ni0.5O2 electrodes follows a parabolic trend with increasing synthesis temperature, which is attributed primarily to the gradual transformation of 3-dimensional (3-D) to 2-dimensional (2-D) diffusion pathways for the Li ions. When synthesized at 400 °C, LxS-LiMn0.5Ni0.5O2 electrodes perform well, benefitting from the 3-D network of channels within the LxS structure. By contrast, when prepared at 500-700 °C, LiMn0.5Ni0.5O2 electrodes operate poorly, which is attributed to the formation of locally disordered structural arrangements that impede Li-ion diffusion. Such an increase in local disorder in the mid-temperature synthesis range is attributed to the structural frustration between the lithium-excess spinal and layered end-members. The transformation from the locally disordered to more ordered layered components between 700 °C and 900 °C enhances electrochemical performance. The study opens new avenues for designing next-generation Mn-rich cathode materials by fine-tuning the synthesis conditions as well as the composition and structure of LxS-LMNO electrodes.

energy storage↗

Applying a Phase-Separation Parameterization in Modeling Secondary Organic Aerosol Formation from Acid-Driven Reactive Uptake of Isoprene Epoxydiols under Humid Conditions

Secondary organic aerosol (SOA) from acid-driven reactive uptake of isoprene epoxydiols (IEPOX) contributes up to 40% of organic aerosol (OA) mass in fine particulate matter. Our previous work showed substantial conversions of particulate inorganic sulfates to surface-active organosulfates (OSs) by IEPOX decreases aerosol acidity and creates a viscous organic-rich shell that poses as a diffusion barrier, inhibiting additional reactive uptake of IEPOX. In this study, to account for this “self-limiting” effect, a phase-separation box model was developed to evaluate parameterizations of IEPOX reactive uptake against time-resolved chamber measurements of IEPOX-SOA tracers, including 2-methyltetrols (2-MT) and methyltetrol sulfates (MTS), at ~ 50% relative humidity. The phase-separation model was most sensitive to mass accommodation coefficient, IEPOX diffusivity in the organic shell, and ratio of the third-order reaction rate constants forming 2-MT and MTS (k MT /k MTS ). In particular, k MT /k MTS had to be lower than 0.1 to bring model predictions of 2-MT and MTS in closer agreement with chamber measurements, while prior studies reported values larger than 0.71. The model-derived rate constants favor more particulate MTS formation due to 2-MT likely off-gassing at ambient-relevant OA loadings. Incorporating this parametrization into chemical transport models is expected to predict lower IEPOX-SOA mass and volatility due to the predominance of OSs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Minimized aging of isocyanurate-based rigid cellular foams for buildings through tailored barrier facers and optimized formulation

The thermal resistivity (h·ft 2 ·°F/Btu/in.) of closed-cell rigid foam insulation materials significantly decreases over time. Diffusion-tight facers are designed to significantly enhance initial thermal resistivity and long-term thermal performance by preventing gas diffusion into the foam cells and inhibiting the escape of low thermal conductivity blowing agents. In addition to the gas diffusion property of the facer film, adhesion between the foam and the facer is crucial for achieving diffusion-tight bonding. Here, this study addresses the challenge of thermal aging by investigating how facer film properties and foam formulation influence the durability of thermal performance. A systematic evaluation was conducted to understand the effects of polymeric barrier films, surface treatments, facer coverage, and foam matrix rigidity on thermal resistivity over time. Key findings reveal that diffusion-tight facers, particularly those with metallized layers and compatible heat seal layers, significantly reduce gas exchange and improve foam-facer adhesion. The optimized system, incorporating barrier facers and a tailored polyisocyanurate foam formulation, achieved initial and aged thermal resistivity values after 200 days of approximately 8.3 and 7.4 h ft 2 ·°F/Btu/in., respectively representing only a 10 % reduction compared to a 17 % reduction observed in control samples without facers. Notably, polyurethane spray foams with facers exhibited only a 4 % reduction in thermal resistivity, compared to a 23 % decrease in control samples, demonstrating nearly six times better retention of thermal performance. This innovative facer technology presents a promising solution for reducing energy costs and represents a significant advancement in optimizing energy management for building envelopes in future technologies. This technology can also be adapted for other applications that necessitate the preservation of long-term thermal performance.

Wanasinghe, Shiwanka Vidarshi [Oak Ridge National ↗

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗