Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SLOPE”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil water percolation and nutrient fluxes as a function of topographical, seasonal and soil texture variation in Central Amazonia, Brazil

Abstract Understanding soil water dynamics and transport of nutrients is challenging in tropical rainforests due to the uniqueness of several properties related to soils, vegetation and seasonality that make relating patterns found in temperate environments to tropical sites difficult. We address the need for better edaphic characterization in tropical environments by investigating soil water percolation rates and chemistry across topographic, soil texture and seasonal gradients in a mature tropical rainforest in Central Amazonia, Brazil. We utilized a passive wick flux meter (e.g., drainage lysimeter) to directly measure real‐time percolation fluxes at 60‐cm depth, and to sample a suite of chemical species across plateau, slope and valley topographic positions. We found percolation flux volume and chemical exports generally increase with decreasing elevation and clay content, which was lowest in the valley. Daily percolation flux was observed to be 2.39 ± 0.44 in plateau, 3.01 ± 0.50 in slope and 6.16 ± 0.83 mm in valley. Most solutes were present in small amounts of <1 mg L −1 , such as PO₄ 3− , Fe 2+ /Fe 3+ and Mn 2+ ; however, NO 3 − concentrations were >20 mg L −1 , even exceeding 100 mg L −1 in the valley. Based on additional isotopic analysis, we speculate high NO 3 − concentrations are partially an artefact of root decomposition following installation of the flux meters. The empirical relationships we show among percolation volume and nutrient exports under varying topographies and soil textures can improve Earth System Model performance by better constraining ecohydrological relationships to nutrient fluxes, which can in‐turn better illuminate the important factors that govern their behaviour.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Utility of near‐surface phenology in estimating productivity and evapotranspiration across diverse ecosystems

Abstract Agroecosystems, which include row crops, pasture, and grass and shrub grazing lands, are sensitive to changes in management, weather, and genetics. To better understand how these systems are responding to changes, we need to improve monitoring and modeling carbon and water dynamics. Vegetation Indices (VIs) are commonly used to estimate gross primary productivity (GPP) and evapotranspiration (ET), but these empirical relationships are often location and crop specific. There is a need to evaluate if VIs can be effective and, more general, predictors of ecosystem processes through time and across different agroecosystems. Near‐surface photographic (red‐green‐blue) images from PhenoCam can be used to calculate the VI green chromatic coordinate (G CC ) and offer a pathway to improve understanding of field‐scale relationships between VIs and GPP and ET. We synthesized observations spanning 76 site‐years across 15 agroecosystem sites with PhenoCam G CC and GPP or ET estimates from eddy covariance (EC) to quantify interannual variability (IAV) in the relationship between GPP and ET and G CC across. We uncovered a high degree of variability in the strength and slopes of the G CC ∼ GPP and ET relationships (R 2 = 0.1 ‐ 0.9) within and across production systems. Overall, G CC is a better predictor of GPP than ET (R 2 = 0.64 and 0.54, respectively), performing best in croplands (R 2 = 0.91). Shrub‐dominated systems exhibit the lowest predictive power of G CC for GPP and ET but have less IAV in slope. We propose that PhenoCam estimates of G CC could provide an alternative approach for predictions of ecosystem processes.

Environmental Sciences & Ecology↗
DOE OSTI2020

Plant‐mediated methane transport in emergent and floating‐leaved species of a temperate freshwater mineral‐soil wetland

Abstract Methane flux from freshwater mineral‐soil (FWMS) wetlands and its variability among sites is largely modulated by plant‐mediated transport. However, plant‐mediated transport processes are rarely resolved in land surface models and are poorly parametrized for plants commonly found in FWMS wetlands. Here, relationships between methane flux and CO 2 uptake, as well as plant conductance of methane were evaluated for three plant species and two characteristic functional types: emergent (narrow‐leaved cattail) and floating‐leaved (American lotus and water lily). We found significant but contrasting correlations between methane flux and CO 2 uptake in cattails ( r 2 = 0.51, slope = −0.16, during morning) and water lily ( r 2 = 0.32, slope = 0.064, after midday). This relationship was not significant in American lotus, showing that stomata regulation of methane fluxes is species‐specific and not generalizable across the floating‐leaved plant functional type. Conductance of methane per leaf area showed distinct seasonal dynamics across species. Conductance was similar among the floating‐leaved species (6.2 × 10 −3 m d −1 in lotus and 7.2 × 10 −3 m d −1 in water lily) and higher than conductance in the emergent species (2.7 × 10 −3 m d −1 ). Our results provide direct observations of plant conductance rates and identify the vegetation parameters (leaf area, stomatal conductance) that modify them. Our results further suggest that models of methane emissions from FWMS should parameterize plant‐mediated transport in different plant functional types, scaled by leaf area and with variable seasonal phenological dynamics, and consider possible species‐specific mechanisms that control methane transport through plants.

Villa, Jorge A.↗
DOE OSTI2022

Gravitational waves from an inflation triggered first-order phase transition

Large excursion of the inflaton field can trigger interesting dynamics. One important example is a first-order phase transition in a spectator sector which couples to the inflaton. Gravitational waves (GWs) from such a first-order phase transition during inflation, an example of an instantaneous source, have an oscillatory feature. In this work, we show that this feature is generic for a source in an era of accelerated expansion. We also demonstrate that the shape of the GW signal contains information about the evolution of the early universe following the phase transition. In particular, the slope of the infrared part of the GW spectrum is sensitive to the evolution of the Hubble parameter when the GW modes reenter the horizon after inflation. The slope of the profile of the intermediate oscillatory part and the ultraviolet part of the GW spectrum depend on the evolution of the Hubble parameter when the modes exit horizon during the inflation and when they reenter the horizon during the reheating. The ultraviolet spectrum also depends on the details of the dynamics of the phase transition. We consider the GW signal in several models of evolution during and after inflation, and compare them with the minimal scenario of quasi- de Sitter inflation followed by radiation domination after a fast reheating, and demonstrate that the shape of the GW can be used to distinguish them. In this way, the GW signal considered in this paper offers a powerful probe to the dynamics of the early universe which is otherwise difficult to explore directly through CMB, large scale structure, big bang nucleosynthesis (BBN), and other well-studied cosmological observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Probing the chiral magnetic wave with charge-dependent flow measurements in Pb-Pb collisions at the LHC

The Chiral Magnetic Wave (CMW) phenomenon is essential to provide insights into the strong interaction in QCD, the properties of the quark-gluon plasma, and the topological characteristics of the early universe, offering a deeper understanding of fundamental physics in high-energy collisions. Measurements of the charge-dependent anisotropic flow coefficients are studied in Pb-Pb collisions at center-of-mass energy per nucleon-nucleon collision $\sqrt{s_{\textrm{NN}}}$ = 5.02 TeV to probe the CMW. In particular, the slope of the normalized difference in elliptic (v 2 ) and triangular (v 3 ) flow coefficients of positively and negatively charged particles as a function of their event-wise normalized number difference, is reported for inclusive and identified particles. The slope ${r}_3^{\textrm{Norm}}$ is found to be larger than zero and to have a magnitude similar to ${r}_2^{\textrm{Norm}}$, thus pointing to a large background contribution for these measurements. Furthermore, ${r}_2^{\textrm{Norm}}$ can be described by a blast wave model calculation that incorporates local charge conservation. In addition, using the event shape engineering technique yields a fraction of CMW (f CMW ) contribution to this measurement which is compatible with zero. This measurement provides the very first upper limit for f CMW , and in the 10–60% centrality interval it is found to be 26% (38%) at 95% (99.7%) confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Metal fluxes at the sediment–water interface in a free water surface constructed wetland

Metal biogeochemistry in the sediment, water, and the sediment–water interface (SWI) was studied in a free water surface constructed wetland. Concentrations of labile copper (Cu), zinc (Zn), sulfate, chloride, and dissolved organic carbon (DOC) were measured with the diffusive gradients in thin film (DGT) and peeper. A good agreement between peeper- and DGT-measured metals was observed for Cu (regression r 2 = 0.3, 95% CI of the slopes > 0), but not for Zn (95% CI of the slopes overlapped with 0), which was attributed to the different complexed compounds between Cu and Zn in porewater. The depth profile of labile metals in sediment porewater varied with time and was consistent with the solid-phase metal deposition, showing higher concentrations in the surface layer (3 to –3 cm) than in the bottom layer (–4 to –13 cm). The depth-averaged labile Cu and Zn concentrations measured by DGT were 1.0 and 3.1 µg/L, and labile sulfate, chloride, and DOC concentrations measured by peeper were 1.8, 3.6, and 2.1 mg/L, respectively. A sharp decrease in sulfate occurred in September when sulfate concentrations became the lowest among sampling months. This was caused by the seasonal sulfur cycles in the wetland, where the dominant sulfur reaction is sulfate reduction in warm seasons and sulfide oxidation in cold seasons. Different metal-removal mechanisms were observed in the two wetland cells; sulfur dynamics controlled the removal processes in the cell without frequent disturbance but failed to influence metal removal in the cell with frequent disturbance due to the interruption of anoxic layers. The flux ratios that compare labile element concentrations between the water column and the SWI (R-Cu, R-Zn, R-DOC, R-sulfate, and R-chloride) were generated to determine metal diffusive fluxes at the interface. Labile Zn was mostly transported from the water to the SWI during all seasons (R-Zn < 1 for all months except January). Labile Cu moved from the SWI to the water during the warm months (R-Cu < 1), which was explained by the bioturbation-induced transport of organic matter based on the positive correlations between R-Cu and R-DOC. In general, sediment can serve either as a sink or a source depending on the environmental conditions, metal speciation, and presence of living organisms. Here, metal flux at the SWI is a key component in the biogeochemical cycling of a constructed wetland.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Attributing differences of solar-induced chlorophyll fluorescence (SIF)-gross primary production (GPP) relationships between two C4 crops: corn and miscanthus

There remains limited information to characterize the solar-induced chlorophyll fluorescence (SIF)-gross primary production (GPP) relationship in C4 cropping systems. The annual C4 crop corn and perennial C4 crop miscanthus differ in phenology, canopy structure and leaf physiology. Investigating the SIF-GPP relationships in these species could deepen our understanding of SIF-GPP relationships within C4 crops. Using in situ canopy SIF and GPP measurements for both species along with leaf-level measurements, we found considerable differences in the SIF-GPP relationships between corn and miscanthus, with a stronger SIF-GPP relationship and higher slope of SIF-GPP observed in corn compared to miscanthus. These differences were mainly caused by leaf physiology. For miscanthus, high non-photochemical quenching (NPQ) under high light, temperature and water vapor deficit (VPD) conditions caused a large decline of fluorescence yield (Φ F ), which further led to a SIF midday depression and weakened the SIF-GPP relationship. The larger slope in corn than miscanthus was mainly due to its higher GPP in mid-summer, largely attributed to the higher leaf photosynthesis and less NPQ. Furthermore, our results demonstrated variation of the SIF-GPP relationship within C4 crops and highlighted the importance of leaf physiology in determining canopy SIF behaviors and SIF-GPP relationships.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mapping tree height in complex terrain of northern China using ultra-high-resolution images

Tree height is a key parameter for estimating forest biomass and carbon sequestration. In recent years, notable progress has been made in mapping tree height using satellite imagery. However, existing tree height products show low accuracy in mountainous and complex terrains, and few studies typically addressed tree height estimations in mountain areas. This study examines the Mentougou district of Beijing, China, characterized by complex terrain and mountainous landscapes. We analyzed two methods for estimating tree height: one using only spectral features and another combining spectral features with topographic factors (elevation, slope, aspect). We used 3-m resolution PlanetScope 8-band multispectral imagery, with 710 field-measured individual tree heights averaged to obtain 471 pixel-level tree height values as ground-truth, to develop tree height prediction models using eXtreme Gradient Boosting (XGBoost), Random Forest (RF), and Gradient Boosting Machine (GBM) models. The results show that the XGBoost model consistently presented the highest accuracy for both methods evaluated. Specifically, the XGBoost model that combined spectral data with elevation and slope variables with an R² of 0.75 and an RMSE of 2.69 m. Using the XGBoost model, we generated the tree height map for the Mentougou area at 3 m resolution, showing tree heights ranging from 0.5 to 30.4 m, and the model’s prediction error standard deviations ranged from 2.50 to 4.71 m, indicating reliable performance across varied terrain. Additionally, we compared and evaluated the global tree height products, identifying limitations in the accuracy within complex terrains. This study demonstrates the potential for accurately predicting tree heights by combining high-resolution multispectral satellites with a terrain factor modeling approach.

Complex terrain↗
DOE OSTI2022

A mass–momentum consistent coupling for mesh-adaptive two-phase flow simulations

Here, we present a novel mass-momentum consistent coupling between a geometric volume-of-fluid scheme and an incompressible flow solver with differing directional-splitting approaches. The advection of the volume fraction is performed using a direction-split algorithm, whereas the momentum advection algorithm uses a traditional unsplit, fractional-step approach. Both algorithms employ finite-volume discretizations based on Cartesian meshes. In solving the mass-momentum consistency problem, momentum fluxes at the cell faces are weighted by the density fluxes based on the already advected volume fraction. The success of our approach lies on introducing a Favre-averaged velocity interpolation at the liquid/gas interface along with a minmod slope limiter. Mesh-convergence studies show that when the minmod slope limiter is used, the two-phase solver retains an accuracy between first and second order, but when a purely upwind scheme is considered, its accuracy drops to first order. Finally, after considering several validation problems, the solver is shown to agree well with reference numerical and experimental data while retaining its robustness and efficiency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal expansion of liquid Fe-S alloy at high pressure

The local structure and density of liquid Fe-S alloys at high pressure have been determined in situ by combined angle and energy dispersive X-ray diffraction experiments in a multi-anvil apparatus, covering a large temperature and compositional range. Precise density measurements collected for increasing temperature allowed us to directly derive the thermal expansion coefficients for liquid Fe-S alloys as a function of composition. In turn, thermal expansion has been used to refine thermodynamic models and to address the crystallization regime of telluric planetary cores by comparing the adiabatic temperature gradient and the slope of the liquidus in the Fe-FeS system. For Fe-S cores of asteroids and small planetesimals, top-down solidification is the dominant scenario as the compositional domain for which the slope of the liquidus is greater than the adiabatic gradient is limited to a narrow portion on the Fe-rich side. However, bottom-up growth of the inner core is expected for S-poor cases, with this compositional domain expanding to more S-rich compositions with increasing pressure (size of the planetary body). In particular, bottom-up crystallization cannot be excluded for the Moon and Ganymede.

36 MATERIALS SCIENCE↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2021

A screening model for predicting injection well pressure buildup and plume extent in CO 2 geologic storage projects

We present a new screening model for predicting injection well pressure buildup and CO 2 plume migration for CO 2 geologic sequestration projects. The model requires only limited information and is quite accurate when compared to detailed simulation results. Such models can assist project developers during the early days of project planning (e.g., for 45Q related projects), and also help regulators perform simple checks against detailed numerical models. The screening model consists of two correlations: one for injectivity index in terms of the slope of the CO 2 fractional flow curve, and the second for total storage efficiency within the footprint of plume as a function of gravity number and slope of the CO 2 fractional flow curve. Using these two correlations, as well as a knowledge of some basic reservoir characteristics and estimates of fluid properties from standard correlations, the injection-well pressure buildup and CO 2 plume extent in the formation can be readily estimated. Finally, the new correlations show a good match with the results of the underlying simulations, and also provide good agreement with independent calculations for two example problems.

42 ENGINEERING↗
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Significance of slight ambient temperature variation on the electrocatalyst performance toward oxygen evolution reaction

Hundreds of oxygen evolution reaction (OER) electrocatalysts have been developed over the past few decades, and their performances are evaluated and compared at ambient temperature. However, the effect of ambient temperature variation on OER electrocatalyst performance has received less attention, which may play a remarkable role in the electrocatalytic activity. In this work, we systematically investigated the influence of ambient temperature variation on electrocatalyst performance toward OER. The results show that the slight ambient temperature variation has a significant effect on OER catalyst performance based on the changes of overpotential (10mAcm -2 ) and Tafel slope. Both remarkable chances are observed on transition metal (Ni) and noble metal (IrO 2 ) electrocatalysts, and the overpotentials decrease around 81 mV with a temperature increase by 20°C (from 10° to 30°C) for both Ni and IrO 2 electrocatalysts with the Tafel slope drops of 36.9 and 29.5 mV dec -1 , respectively. Further, a similar trend is also found in the electrochemically active surface (ECSA) normalized performance and the charge transfer resistance. This study demonstrates that reporting the actual operating temperature for OER is not only recommended but also necessary to evaluate and compare electrocatalyst activities from different materials systems properly, and neglecting the ambient temperature variation effect can highly mislead conclusions.

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

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗
DOE OSTI2025

High-pressure Cr 3+ luminescence and Raman spectroscopy of a natural MgAl 2 O 4 spinel to ∼60 GPa

Cr 3+ luminescence and Raman spectroscopy of a natural spinel were collected up to ∼60 GPa at ambient temperature using a Ne or a 4:1 methanol-ethanol mixture as the pressure transmitting media. Here, in the Ne experiments above ∼40 GPa, and in the 4:1 methanol–ethanol experiments above ∼20 GPa, all luminescence emission lines shift non-linearly, and by 59 GPa the spectra are broadened and are significantly changed. The Raman modes shift smoothly to higher frequencies with pressure, with a slight change in slope above 40 GPa, with no kinks or discontinuities. By ∼60 GPa, only the most intense Eg mode remains resolvable, appearing broadened and reduced in intensity. The change in slope in the Raman modes and luminescence bands may be due to a change in the compression mechanism of spinel. When spinel was decompressed from 40 GPa all changes were fully reversible; however, when the samples were decompressed from 54 to 60 GPa in Ne and 59 GPa in 4:1 methanol-ethanol the spectra were irreversibly changed. The change in the luminescence spectra and the Raman spectra are likely due to a pressure induced disorder of the anions or partial amorphization between 40 and 60 GPa.

Cr3+↗