Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupling constant”

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

Interfacial Dzyaloshinskii-Moriya interaction of antiferromagnetic materials

The interface between a ferromagnet (FM) or antiferromagnet (AFM) and a heavy metal (HM) results in an antisymmetric exchange interaction known as the interfacial Dzyaloshinskii-Moriya interaction (iDMI) which favors non-collinear spin configurations. The iDMI is responsible for stabilizing noncollinear spin textures such as skyrmions in materials with bulk inversion symmetry. Interfacial DMI values have been previously determined theoretically and experimentally for FM/HM interfaces, and, in this work, they are calculated for the metallic AFM MnPt and the insulating AFM NiO. The heavy metals considered are W, Re, and Au. Values for the iDMI, exchange, and anisotropy constants are determined for different AFM and HM thicknesses. The iDMI values of the MnPt heterolayers are comparable to those of the common FM materials, and those of NiO are lower. In few-layer films of (001) MnPt, the high spin orbit coupling of the Pt layers can give rise to a small DMI in the absence of a HM layer.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Helimagnetism in MnBi 2 Se 4 Driven by Spin-Frustrating Interactions Between Antiferromagnetic Chains

We report the magnetic properties and magnetic structure determination for a linear-chain antiferromagnet, MnBi 2 Se 4 . The crystal structure of this material contains chains of edge-sharing MnSe 6 octahedra separated by Bi atoms. The magnetic behavior is dominated by intrachain antiferromagnetic (AFM) interactions, as demonstrated by the negative Weiss constant of –74 K obtained by the Curie–Weiss fit of the paramagnetic susceptibility measured along the easy-axis magnetization direction. The relative shift of adjacent chains by one-half of the chain period causes spin frustration due to interchain AFM coupling, which leads to AFM ordering at TN = 15 K. Neutron diffraction studies reveal that the AFM ordered state exhibits an incommensurate helimagnetic structure with the propagation vector k = (0, 0.356, 0). The Mn moments are arranged perpendicular to the chain propagation direction (the crystallographic b axis), and the turn angle around the helix is 128°. The magnetic properties of MnBi2Se4 are discussed in comparison to other linear-chain antiferromagnets based on ternary mixed-metal halides and chalcogenides.

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

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Abstract The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Efficient continuous Energy-Multigroup hybrid depletion scheme using the Shift Monte Carlo code. Part I: Energy condensation sensitivity analysis

Monte Carlo (MC) codes coupled to depletion solvers are increasingly used to provide high fidelity fuel cycle modeling capabilities. Here, these coupled depletion-MC tools produce accurate results in general but can experience nonphysical spatial oscillations when time steps are large or when a system’s dominance ratio approaches unity. Two substepping techniques have been developed previously to remedy and dampen these spatial oscillations without needing to reduce step sizes. The first approach relied on higher-order techniques to account for spectral changes within steps (extrapolation and interpolation techniques). The second approach used the first order perturbation (FOP) theory to account for the change in the one-group spatial flux distribution within steps. This paper develops a hybrid depletion methodology which, in a way, combines how the flux is handled in both substepping techniques. Specifically, the multigroup (MG) MC Shift code is used to update the flux distribution within steps rather than a one-group FOP solver. A fully reflected pincell is investigated, which is not spatially dependent in the MG representation. Thus, the analysis in this paper is an initial demonstration of hybrid depletion. An upcoming companion paper will focus on how the hybrid depletion dampens spatial oscillations. The hybrid depletion approach is verified to be consistent with previous constant extrapolation depletion (CED) methods. This paper finds that the hybrid CED exhibits some error in the eigenvalue and one group constants within macro steps. To address this discrepancy, a simple interpolation scheme (CELI) is investigated. This work found that CELI sufficiently addresses the discrepancy in spectrum for macro steps up to 100 days. Overall, this work demonstrates that the hybrid depletion method can significantly reduce the number of high fidelity MC executions in a MC-coupled depletion with an acceptable eigenvalue error.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Electron-phonon coupling induced intrinsic Floquet electronic structure

Floquet states are a topic of intense contemporary interest, which is often induced by coherent external oscillating perturbation (e.g., laser, or microwave) which breaks the continuous time translational symmetry of the systems. Usually, electron–phonon coupling modifies the electronic structure of a crystal as a non-coherent perturbation and seems difficult to form Floquet states. Surprisingly, we found that the thermal equilibrium electron–phonon coupling in M(MoS) 3 and M(MoSe) 3 (where M is a metallic element) exhibits a coherent behavior, and the electronic structure can be described by the Floquet theorem. Such a coherent Floquet state is caused by a selective giant electron–phonon coupling, with thermodynamic phonon oscillation serving as a driving force on the electronic part of the system. The quasi-1D Dirac cone at the Fermi energy has its band gap open and close regularly. Similarly, the electric current will oscillate even under a constant voltage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the fusion yield dependencies in OMEGA DT-layered implosion experiments using a physics-based statistical mapping model

Improving the performance of inertial confinement fusion implosions requires physics models that can accurately predict the response to changes in the experimental inputs. Good predictive capability has been demonstrated for the fusion yield using a statistical mapping of simulated outcomes to experimental data [Gopalaswamy et al., Nature 565(771), 581–586 (2019)]. In this paper, a physics-based statistical mapping approach is used to extract and quantify all the major sources of degradation of fusion yield for direct-drive implosions on the OMEGA laser. Here, the yield is found to be dependent on the age of the deuterium tritium fill, the ℓ = 1 asymmetry in the implosion core, the laser beam-to-target size ratio, and parameters related to the hydrodynamic stability. A controlled set of experiments were carried out where only the target fill age was varied while keeping all other parameters constant. The measurements were found to be in excellent agreement with the fill age dependency inferred using the mapping model. In addition, a new implosion design was created, guided by the statistical mapping model by optimizing the trade-offs between increased laser energy coupling at larger target size and the degradations caused by the laser beam-to-target size ratio and hydrodynamic instabilities. When experimentally performed, an increased fusion yield was demonstrated in targets with larger diameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Long-lived axionlike particles from tau decays

Axionlike particles (ALPs) are well-motivated examples of light, weakly coupled particles in theories beyond the Standard Model. In this work, we study long-lived ALPs coupled exclusively to leptons in the mass range between 2⁢𝑚 𝑒 and 𝑚 𝜏 − 𝑚 𝑒 . For anarchic flavor structure the leptophilic ALP production in tau decays or from ALP-tau bremsstrahlung is enhanced thanks to derivative couplings of the ALP and can surpass production from electron and muon channels, especially for ALPs heavier than 𝑚 𝜇 . Using past data from high-energy fixed-target experiments such as CHARM and BEBC we place new constraints on the ALP decay constant 𝑓 𝑎 , reaching scales as high as 𝒪⁡(10 8 ) GeV in lepton-flavor-violating channels and 𝑓 𝑎 ∼ 𝒪⁡(10 2 ) GeV in lepton-flavor-conserving ones. We also present projections for the event-rate sensitivity of current and future detectors to ALPs produced at the Fermilab Main Injector, the CERN SPS, and in the forward direction of the LHC. We show that SHiP will be sensitive to 𝑓 𝑎 values that are over an order of magnitude above the existing constraints.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Tensor network formulation of two-dimensional gravity

We show how to formulate a lattice gauge theory whose naive continuum limit corresponds to two-dimensional (Euclidean) quantum gravity including a positive cosmological constant. More precisely the resultant continuum theory corresponds to gravity in a first-order formalism in which the local frame and spin connection are treated as independent fields. Recasting this lattice theory as a tensor network allows us to study the theory at strong coupling without encountering a sign problem. In two dimensions this tensor network is exactly soluble and we show that the system has a series of critical points that occur for pure imaginary coupling and are associated with first order phase transitions. We then augment the action with a Yang-Mills term which allows us to control the lattice spacing and show how to apply the tensor renormalization group to compute the free energy and look for critical behavior. Finally we perform an analytic continuation in the gravity coupling in this extended model and show that its critical behavior in a certain scaling limit depends only on the topology of the underlying lattice. We also show how the lattice gauge theory can be naturally generalized to generate the Polyakov or Liouville action for two dimensional quantum gravity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Gifts from anomalies: Exact results for Landau phase transitions in metals

Non-Fermi liquid phenomena arise naturally near critical points of Landau ordering transitions in metallic systems, where strong fluctuations of a bosonic order parameter destroy coherent quasiparticles. Despite progress in developing controlled perturbative techniques, much of the low energy physics of such metallic quantum critical points remains poorly understood. We demonstrate that exact, non-perburbative results can be obtained for both optical transport and static susceptibilities in “Hertz-Millis” theories of Fermi surfaces coupled to critical bosons. Such models possess a large emergent symmetry and anomaly structure, which we leverage to fix these quantities. In particular, we show that in the infrared limit, the boson self energy at zero wave vector, \mathbf{q}=0 𝐪 = 0 , is a constant independent of frequency, and the real part of the optical conductivity, \sigma(\omega) σ ( ω ) , is purely a delta function Drude peak with no other corrections. Therefore, further frequency dependence in the boson self energy or optical conductivity can only come from irrelevant operators in a clean system. Exact relations between Fermi liquid parameters as the critical point is approached from the disordered phase are also obtained. The absence of a universal, power law frequency dependence in the boson self energy contrasts with previous perturbative calculations, and we explain the origin of this difference.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2024

Bioelectrochemical reactor to manage anthropogenic sulfate pollution for freshwater ecosystems: Mathematical modeling and experimental validation

Anthropogenic sulfate loading into otherwise low-sulfate freshwater systems can cause significant ecological consequences as a biogeochemical stressor. To address this challenge, in situ bioremediation technologies have been developed to leverage naturally occurring microorganisms that transform sulfate into sulfide rather than implementing resource-intensive physio-chemical processes. However, bioremediation technologies often require the supply of electron donors to facilitate biological sulfate reduction. Bioelectrochemical systems (BES) can be an alternative approach for supplying molecular hydrogen as an electron donor for sulfate-reducing bacteria through water electrolysis. Although the fundamental mechanisms behind BESs have been studied, limited research has evaluated the design and operational parameters of treatment systems when developing BESs on a scale relevant to environmental systems. Further, this study aimed to develop an application-based mathematical model to evaluate the performance of BESs across a range of reactor configurations and operational modes. The model was based on sulfate transformation by hydrogenotrophic sulfate-reducing bacteria coupled with the recovery of solid iron sulfide species formed by the oxidative dissolution of dissolved ferrous iron from a stainless steel anode. Sulfate removal closely corresponded to the rate of electrolytic hydrogen production and hydraulic residence time but was less sensitive to specific microbial rate constants. The mathematical model results were compared to experimental data from a pilot-scale BES tested with nonacidic mine drainage as a case study. The close agreement between the mathematical model and the pilot-scale BES experiment highlights the efficacy of using a mathematical model as a tool to develop a conceptual design of a scaled-up treatment system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Bimetallic Metal–Organic Framework Fe/Co-MIL-88(NH 2 ) Exhibiting High Peroxidase-like Activity and Its Application in Detection of Extracellular Vesicles

Metal–organic frameworks (MOFs) have many attractive features, including tunable composition, rigid structure, controllable pore size, and large specific surface area, and thus are highly applicable in molecular analysis. Depending on the MOF structure, a high number of un-saturated metal sites can be exposed to catalyze chemical reactions. In the present work, we report that by using both Co(II) and Fe(III) to prepare the MIL-88(NH 2 ) MOF, we can produce the bimetallic MOF that can catalyze the conversion of 3,3', 5,5"-tetramethylbenzidine (TMB) to a color product through reaction with H 2 O 2 at a higher reaction rate than the monometallic Fe-MIL-88(NH 2 ). The Michaelis constants (K m ) of the catalytic reaction for TMB and H 2 O 2 are 3-5 times smaller, and the catalytic constants (k cat ) are 5-10 times higher than those of the horse-radish peroxidase (HRP), supporting ultrahigh peroxidase-like activity. These values are also much more superior to those of the HRP-mimicking MOFs reported previously. Interestingly, the bimetallic MOF can be coupled with glucose oxidase (GOx) to trigger the cascade enzymatic reaction for highly sensitive detection of extracellular vesicles (EVs), a family of important biomarkers. Through conjugation to the aptamer that recognizes the marker protein on EV surface, the MOF can help isolate the EVs from biological matrices, which are subsequently labeled by GOx via antibody recognition. The cascade enzymatic reaction between MOF and GOx enables detection of EVs at a concentration as low as 7.8 × 10 4 particles/ml. Further, the assay can be applied to monitor EV secretion by cultured cells, and also can successfully detect the different EV quantities in the sera samples collected from cancer patients and healthy controls. Overall, we prove that the bimetallic Fe/Co-MIL-88(NH 2 ) MOF, with its high peroxidase-activity and high biocompatibility, is a valuable tool deployable in clinical assays to facility disease diagnosis and prognosis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗
DOE OSTI2020

Wireless Interrogation of High Temperature Surface Acoustic Wave Dynamic Strain Sensor

Dynamic strain sensing is necessary for high-temperature harsh-environment applications, including powerplants, oil wells, aerospace, and metal manufacturing. Monitoring dynamic strain is important for structural health monitoring and condition-based maintenance in order to guarantee safety, increase process efficiency, and reduce operation and maintenance costs. Sensing in high-temperature (HT), harsh-environments (HE) comes with challenges including mounting and packaging, sensor stability, and data acquisition and processing. Wireless sensor operation at HT is desirable because it reduces the complexity of the sensor connection, increases reliability, and reduces costs. Surface acoustic wave resonators (SAWRs) are compact, can operate wirelessly and battery-free, and have been shown to operate above 1000°C, making them a potential option for HT HE dynamic strain sensing. This paper presents wirelessly interrogated SAWR dynamic strain sensors operating around 288.8MHz at room temperature and tested up to 400°C. The SAWRs were calibrated with a high-temperature wired commercial strain gauge. The sensors were mounted onto a tapered-type Inconel constant stress beam and the assembly was tested inside a box furnace. The SAWR sensitivity to dynamic strain excitation at 25°C, 100°C, and 400°C was .439 μV/με, 0.363μV/με, and .136 μV/με, respectively. Furthermore, the experimental outcomes verified that inductive coupled wirelessly interrogated SAWRs can be successfully used for dynamic strain sensing up to 400°C.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Consequence Analysis of Residual Water in a Storage Canister - 20484

Recent observations from projects sponsored under the DOE-Nuclear Energy Spent Fuel and Waste Disposition campaign show that residual free water, well above the 0.4 gm-moles that correspond to a 3 torr rebound pressure after a 30-minute hold following prototypic drying, may remain within SNF canisters used for long-term dry storage. Specifically, experimental observations in an Integrated Research Project showed subzero temperatures and a local pocket of frozen water in a mockup fuel assembly subjected to an industry practiced drying process. In gas sampling and analysis from the High Burnup (HBU) Demonstration project, using a fully loaded and dried SNF canister, indicated presence of 100 g of residual water had remained in the canister. The effects of residual water in terms of potential adverse impacts on the condition of the SNF and the canister internals during long-term dry storage has been rigorously evaluated to address inadvertent residual water contents. The impact of residual water on the fuel and canister internals, with a focus on the oxidation of the cladding and of postulated exposed (breached-clad) fuel, were evaluated using an integrated set of time-dependent environmental conditions with semi-empirical materials oxidation models. An 'integration model,' with coupled sub-models, developed and previously used by the Center for Nuclear Waste Regulatory Analysis (CNWRA), was applied and expanded with improved sub-models for water radiolysis, and for cladding and fuel oxidations. Cladding oxidation models for various alloys were developed using the rate constants for various cladding alloys in the dry storage temperature range. A burnup-dependent fuel oxidation model was developed; it accounted for inhibition effect of fission products on fuel oxidation rates. The potential for extension of the initial breach flaw (1 mm hole) was evaluated by two separate mechanics criteria. In the first case, a critical cladding strain criterion was used. In the second case, a critical stress intensity was used to evaluate stability of a planar axial flaw postulate. A failure condition exists when a sufficient amount of fuel oxidation (at the postulated breached-clad location) causes fuel pellet swelling to load the cladding to meet a mechanics criterion for flaw extension. The integration model's results provide the following insights. The radiolysis-generated oxygen is consumed with a partitioning between the cladding and the fuel. The extent of cladding oxidation is no more than 2 μm additional consumption of the cladding metal even with 10 moles of residual water, indicating that changes in cladding conditions due to water are expected to be negligible. The conditions for cladding failure (breach opening size extension) directly correlates with the extent of fuel oxidation. Fuel oxidation significant to extend the breach could occur even with 5.5 moles of residual water. The results indicate that fuel oxidation vis-a-vis cladding oxidation is more likely under a higher radiation field, i.e., when radiolytic decomposition of the residual water occurs in few years compared to several decades, and that the conditions to cause failure are increased with fuel and cladding temperatures approaching the peak storage temperature limit of 400 deg. C. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Thermodynamically consistent algorithms for models of incompressible multiphase polymer solutions with a variable mobility

Here we present a general strategy for developing structure and property preserving numerical algorithms for thermodynamically consistent models of incompressible multiphase polymer solutions with a variable mobility. We first present a formalism to derive thermodynamically consistent, incompressible, multiphase polymer models. Then, we develop the general strategy, known as the supplementary variable method, to devise thermodynamically consistent numerical approximations to the models. We illustrate the numerical strategy using newly developed models of incompressible diblock copolymer solutions coupled with an electric and a magnetic field, respectively. Mesh refinement is conducted to verify convergence rates of the developed schemes. Some numerical examples are given to exhibit underlying dynamics absent from and driven by the external fields, respectively, highlighting differences between models with the variable and constant mobilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Pore connectivity influences mass transport in natural rocks: Pore structure, gas diffusion and batch sorption studies

For this work, six rocks (one granodiorite, one limestone, two chalks, one mudstone, and one dolostone) with different extents of heterogeneity at six different particle sizes (from 75 to 8000 μm) were studied to describe the effects of pore connectivity on mass transport. The methods applied were (i) porosity measurement of granular rocks, (ii) analyses of gas-phase diffusive transport in a bed of packed particles, along with a solid quartz method at these six particle sizes being developed to identify the contribution of intraparticle diffusion, and (iii) batch sorption tests of multiple ions (anions and cations) with subsequent analyses of inductively coupled plasma-mass spectrometry. Granular porosity measurement results reveal that with decreasing particle sizes, the effective porosities for the “heterogenous” group of rocks (Grimsel granodiorite and Edwards limestone) increase, whereas the porosities of another “homogeneous” group (two Israel chalk samples, Japan mudstone, and Wyoming dolostone) remain constant. Gas diffusion results show that the intraparticle gas diffusion coefficient among these two sample groups, varying in the magnitude of 10 -8 to 10 -6 m 2 /s, are not directly correlated to the porosity differences. Moreover, the batch sorption work displays a different affinity of rocks for various tracers. For Grimsel granodiorite, Japan mudstone, and Wyoming dolostone, the adsorption capacity of Sm 3+ and Eu 3+ increases as the particle size decreases. In general, this integrated research of grain size distribution, granular rock porosity, intraparticle diffusivity, and ionic sorption capacity gives insights into the pore connectivity effect on both physical and chemical transport behaviors for different lithologies and/or different particle sizes.

58 GEOSCIENCES↗