Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow regime”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Evaluation of idealized large-eddy simulations performed with the Weather Research and Forecasting model using turbulence measurements from a 250 m meteorological mast

Abstract. We investigate the ability of the Weather Research and Forecasting model to perform large-eddy simulation of canonical flows. This is achieved through comparison of the simulation outputs with measurements from sonic anemometers on a 250 m meteorological mast located at Østerild, in northern Denmark. Østerild is on a flat and rough area, and for the predominant wind directions, the atmospheric flow can be considered to be close to homogeneous. The idealized simulated flows aim at representing atmospheric boundary layer turbulence under unstable, neutral, and stable stability conditions at the surface, which are statistically significant conditions observed at Østerild. We found that the resolved fields from the simulations appear to have the characteristics of the three stability regimes. Vertical profiles of observed mean wind speeds and direction are well reproduced by the simulations, with the largest differences under near-neutral conditions, where the effect of the subgrid-scale model is evident on the vertical wind shear close to the surface. Vertical profiles of observed eddy fluxes are also well reproduced by the simulations, with the largest differences for the three velocity component variances under stable stability conditions, although nearly always within the observed variability. With regards to turbulent kinetic energy, we find good agreement between observations and simulations at all vertical levels. Simulated and observed velocity spectra match very well and show very similar behavior with height and with atmospheric stability within the low-frequency interval; at the effective resolution, the simulated spectra show the typical drop-off of finite differences. Our findings demonstrate that these idealized simulations reproduce the characteristics of atmospheric stability regimes often observed at a high turbulent and flat site within a direction sector, where the air flows over nearly homogeneous land.

17 WIND ENERGY↗
DOE OSTI2023

Influence of Temperature on Properties and Dynamics of Gas-Solid Flow in Fluidized-Particle Tubular Solar Receiver

A fluidized particle single-tube solar receiver has been tested for investigating the gas-particle characteristics that enable the best operating conditions in a commercial-scale plant. The principle of the solar receiver is to fluidize the particles in a vessel – the dispenser – in which the receiver tube is plunged. The particles are flowing upward in the tube, irradiated over 1-meter height, by applying an overpressure in the dispenser. Experiments with a concentrated solar flux varying between 188 and 358 kW/m² are carried out, and the particle mass flux varied from 0 to 72 kg/(m²s). The mean particles and external tube wall temperatures in the irradiated zone are heated from the ambient to respectively 700°C and 940°C. It is shown that the temperature rise leads to a decrease of the particle volume fraction. Furthermore, a self-regulation of the system is evidenced with a short transient regime. This characteristic is essential from the operational viewpoint. The thermal efficiency of the receiver increases with the particle flow rate, reaching between 60 and 75% above 30 kg/(m²s). Several fluidization regimes are identified thanks to pressure signal analyses, like slugging, turbulent and fast fluidization, showing that regimes transitions are strongly affected by the temperature.

Gueguen, Ronny (ORCID:0000000252368638)↗
DOE OSTI2021

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗
DOE OSTI2025

Automated nuclear cloud feature extraction from film

Chemical, biological, radiological, nuclear, and explosives incidents require rapid detection and characterization for appropriate response. For a nuclear detonation, visible-light cameras may be used to locate the cloud and characterize fallout deposition when coupled with numerical models. Films from the United States’ nuclear testing era compose the only sizeable collection of imagery depicting high-yield detonations. These films offer unique insights into characteristics of flows involving scales that are difficult to replicate experimentally, and they are a valuable source of data for the validation of models for nuclear fallout transport, either as part of emergency response or forensic activities. In this work, we implement modern computer vision and machine learning techniques to identify and track the cloud automatically and subsequently determine the time dependence of some of its features. We trained a ResNet-18 image classifier on hundreds of images to categorize nuclear cloud morphology. Each category or cloud regime is determined by early cloud evolution and is associated to constitutive properties of the flow, such as distribution of vorticity. Next, we identified keypoint features using the KAZE algorithm and tracked these keypoints in the images, allowing us to determine the dimensions and velocities of the cloud across film frames. These measurements converted to real-world units provide valuable experimental data that can be used in the development and validation of nuclear cloud models. We compared the results of this method against manual cloud rise measurements from two different films. In one, our automated method accelerated the feature extraction process without sacrificing measurement accuracy.

Khristy, Joel [ORNL] (ORCID:0000000209963060)↗
DOE OSTI2023

Optimization with Neural Network Feasibility Surrogates: Formulations and Application to Security-Constrained Optimal Power Flow

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Directing the far-from-equilibrium assembly of nanoparticles in confined liquid crystals by hydrodynamic fields

The assembly of nematic colloids relies on long-range elastic interactions that can be manipulated through external stimuli. Confinement and the presence of a hydrodynamic field alter the defect structures and the energetic interactions between the particles. Here, the assembly landscape of nanoparticles embedded in a nematic liquid crystal confined in a nanochannel under a pressure-driven flow is determined. The dynamics of the liquid crystal tensor alignment field is determined through a Poisson-Bracket framework, namely the Stark–Lubensky equations, coupled with the zero-Reynolds momentum equations and the liquid crystal Landau-de Gennes free energy functional. A second order semi-implicit time integration and a three-dimensional Galerkin finite element method are used to resolve flow and nematic fields under several conditions. In general, the zero Reynolds flow displaces the defects around the particles in the upstream direction and renders the surface anchoring ineffective when the flow strength dominates over the nematic elasticity. More importantly, the potential of mean force for particle assembly is non-monotonic independent of surface anchoring. Our results show that the confinement length scale determines the repulsion/attraction transition between colloids, while the flow strength modifies the static defect structure surrounding the particles and determines the magnitude of the energetic barrier for successful assembly. In the attractive regime, the particles move at different rates through the nematic until one particle eventually catches up with the other. This process occurs against or along the direction of flow depending on the flow strength. Ultimately, these results provide a template for engineering and controlling the transport and assembly of nanoparticles under far-from equilibrium conditions in anisotropic media.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular Dynamics Simulation and Theoretical Analysis of Structural Relaxation, Bond Exchange Dynamics, and Glass Transition in Vitrimers

Vitrimers are a class of polymer networks featuring dynamic covalent cross-links that can undergo associative bond exchange. There has been recent interest in these materials due to their promise as recyclable thermosets or self-healing polymers because of the ability of vitrimer networks to rearrange at the molecular level and undergo macroscopic flow. However, the practical use of these materials often occurs in the supercooled regime or glassy state, where the implications of dynamic bonds are complicated by the interplay between slow activated segmental dynamics, cross-link (i.e., bond-exchange) kinetics, and ultimately material properties. Here, in this paper, we combine coarse-grained molecular dynamics simulation and microscopic statistical mechanical theory to understand how cross-linking kinetics affect material dynamics and how this couples to segmental relaxation of the polymeric network strands across a spectrum of length and time scales, especially in the supercooled regime. We characterize the Kuhn segmental alpha relaxation time and bond exchange time for vitrimer systems across various cross-link densities, temperatures, and bond exchange rates. Simulation and theory both exhibit a bending-up behavior for bond exchange time upon cooling, suggesting a coupling between bond exchange dynamics and segmental relaxation that intensifies with faster bond exchange kinetics. We also found bond exchange dynamics have an impact on Kuhn segment alpha relaxation time, which is most significant at higher cross-link densities. Both these effects are most prominent when the bond exchange time is similar to the Kuhn segment alpha relaxation time, and the resulting coupling of these two relaxation processes is tied to both the probability of a free end to find a bonded pair and the time scale of the constraints imposed by the dynamic cross-links. This relationship is reflected by a cross-link dependence of a theoretical parameter which represents the quantitative degree of coupling between bond exchange and segmental dynamics. Overall, the combination of simulation and theory clarifies the intricate interaction between bond kinetics and segmental relaxation and demonstrates the ability to provide molecular-level insights into vitrimer dynamics over a wide temperature range.

dynamic relaxation↗
DOE OSTI2022

Birefringence tomography for axion cloud

An axion cloud surrounding a supermassive black hole can be naturally produced through the superradiance process. Its existence can be examined by the axion induced birefringence effect. It predicts an oscillation of the electric vector position angle of linearly polarized radiations. Stringent constraints of the existence of the axion in a particular mass window has been obtained based on the recent Event Horizon Telescope measurement on M87$^\star$. The future Very-Long-Baseline Interferometry (VLBI) observations will be able to measure the vicinity of many supermassive black holes, thus it opens the possibility to search for the existence of axions in a wide mass regime. In this paper, we study how different black hole properties and accretion flows influence the signatures of the axion induced birefringence. We include the impacts of black hole inclination angles, spins, magnetic fields, plasma velocity distributions, the thickness of the accretion flows. We pay special attention to characterize the washout effects induced by the finite thickness of the accretion flows and the lensed photons. Finally, based on this study, we give prospects on how to optimize the axion search using future VLBI observations, such as the next-generation Event Horizon Telescope, to further increase the sensitivity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Mean-field control for efficient mixing of energy loads

In this work, we pose an engineering challenge of controlling an ensemble of energy devices via coordinated, implementation-light, and randomized on/off switching as a problem in nonequilibrium statistical mechanics. We show that mean-field control with nonlinear feedback on the cumulative consumption, assumed available to the aggregator via direct physical measurements of the energy flow, allows the ensemble to recover from its use in the demand response regime, i.e., transition to a statistical steady state, significantly faster than in the case of the fixed feedback. Moreover when the nonlinearity is sufficiently strong, one observes the phenomenon of “super-relaxation,” where the total instantaneous energy consumption of the ensemble transitions to the steady state much faster than the underlying probability distribution of the devices over their state space, while also leaving almost no devices outside of the comfort zone.

25 ENERGY STORAGE↗
DOE OSTI2021

A High-Fidelity Modeling Tool to Support the Design of Oxy-Combustors for Direct-Fired Supercritical CO2 Cycles

The challenge in the design of oxy-combustors for direct-fired supercritical CO 2 (sCO 2 ) cycles is in addressing disparate performance metrics and objectives. Key design parameters to consider include, among others, injector design for mixing and flame stability, split of recycled CO 2 diluent between injectors and cooling films, target flame temperatures to control noncondensable products, and strategies to inject the diluent CO 2 for film cooling and thermal control. In order to support novel oxy-combustor designs, a high-fidelity yet numerically efficient modeling framework based on the CRUNCH CFD® flow solver is presented, featuring key physics-based submodels relevant in this regime. For computational efficiency in modeling large kinetic sets, a flamelet/progress variable (FPV) based tabulated-chemistry approach is utilized featuring a three-stream extension to allow for the simulation of the CO 2 film cooling stream in addition to the fuel and oxidizer streams. Finite rate chemistry effects are modeled in terms of multiple progress variables for the primary flame as well as for slower-evolving chemical species such as NO x and SO x contaminants. Real fluid effects are modeled using advanced equations of states. The predictive capabilities of this computationally tractable design support tool are demonstrated on a conceptual injector design for an oxy-combustor operating near 30 MPa. Simulations results provide quantitative feedback on the effectiveness of the film cooling as well as the level of contaminants (CO, NO, and N) in the exhaust due to impurities entering from the injectors. These results indicate that this framework would be a useful tool for refining and optimizing the oxy-combustor designs as well as risk mitigation analyses.

Engineering↗
DOE OSTI2020

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Visualization of fast ion phase-space flow in plasmas well-below, near and well-above Alfvén eigenmode stability threshold in tokamak

Abstract Transport of fast ions along certain local phase space paths, referred to as fast ion phase-space flow, has been systematically measured by an imaging neutral particle analyzer (INPA) in three plasma regimes, which are well-below, near, and well-above the Alfvén eigenmode stability threshold (Du et al 2021 Phys. Rev. Lett. 127 235002). (1) In plasmas well-below the Alfvén eignenmodes (AE) stability threshold, fast ions are well-confined on passing particle orbits without noticeable transport over the phase space. The observed INPA images agree well with the synthetic INPA images, using the fast ion distribution predicted by neoclassical theory. (2) In plasmas near the AE stability threshold, INPA images in the presence of AE activity moderately deviate from those without AE activity. The image difference can be well interpreted by AE-driven, phase-space fast ion flow. Paths of this flow over the velocity space (or streamlines) are reconstructed by the intersection lines of curved E ʹ and µ surfaces, referred to as E ʹ and µ line (where E ′ ≡ E − ( ω / n ) P ζ ; E , P ζ and µ are the energy, canonical toroidal momentum and magnetic moment of ions; ω and n are the angular frequency and toroidal mode number of AEs, respectively). Resonant fast ions move radially inward by gaining energy and move radially outward by losing energy and the trajectory well aligns with the E ʹ and µ lines that pass through the mode resonances near the injection energy of neutral beams. (3) In plasmas well-above the AE stability threshold, fast ion phase-space dynamics shows additional features. Fast ions are transported out of the birth positions so promptly along the streamlines that the slowing-down process from the injection energy is not observable, exhibiting strong critical gradient behavior at local phase space. As a result, the increase of electron temperature is very small, in spite of an increase of beam power by ∼ 45 % . It should be emphasized that the directions of phase-space transport, induced by AEs with different frequencies, structures and mode numbers, do not largely differ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A critical dislocation velocity for serration mechanism transition in a nickel-chromium solid solution alloy

Here, the influence of strain rate across three orders of magnitude (1.70 × 10 –5 /s to 1.43 × 10 –2 /s) along with the effect of the plastic strain accumulation (up to 10%) on the serrated plastic flow were investigated in the nickel-chromium (Ni-Cr) solid solution alloy Nimonic 75 by performing constant-strain-rate tension testing at 600 °C. As the strain rate decreased, the critical strain for the onset of serrations transitioned from normal behavior to inverse behavior. The serrated flow was characterized as Type A+B serration at high strain rate (1.43 × 10 –2 /s). In the intermediate strain-rate regime (1.43 × 10 –3 /s and 1.45 × 10 –4 /s), Type B serrations were observed and followed by a transformation to Type C+B serrations. At the low strain rate (1.70 × 10 –5 /s), the plastic flow immediately displayed Type C serrations, which later evolved into Type C+B serrations. Regardless of the strain rate, plastic strain, or dislocation density, a critical dislocation velocity falling in the range of 1.2 × 10 –6 – 2.2 × 10 –6 m/s was identified to signify the onset of Type C serration, whereby the mobile dislocations break free from the solute cloud for short bursts of deformation. Finally, a novel model by solute rearrangement across dislocation cores was used to understand how the critical dislocation velocity is quantitatively determined by the rate at which solute atoms are able to hop across the glide plane as a partial dislocation core moves through the lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Experimental study of the edge radial electric field in different drift configurations and its role in the access to H-mode at ASDEX Upgrade

The formation of the equilibrium radial electric field (Er) has been studied experimentally at ASDEX Upgrade (AUG) in L-modes of “favorable” (ion ∇ B-drift toward primary X-point) and “unfavorable” (ion ∇ B-drift away from primary X-point) drift configurations, in view of its impact on H-mode access, which changes with drift configurations. Edge electron and ion kinetic profiles and impurity velocity and mean-field Er profiles across the separatrix are investigated, employing new and improved measurement techniques. The experimental results are compared to local neoclassical theory as well as to a simple 1D scrape-off layer (SOL) model. It is found that in L-modes of matched heating power and plasma density, the upstream SOL Er and the main ion pressure gradient in the plasma edge are the same for either drift configurations, whereas the Er well in the confined plasma is shallower in unfavorable compared to the favorable drift configuration. The contributions of toroidal and poloidal main ion flows to Er, which are inferred from local neoclassical theory and the experiment, cannot account for these observed differences. Furthermore, it is found that in the L-mode, the intrinsic toroidal edge rotation decreases with increasing collisionality and it is co-current in the banana-plateau regime for all different drift configurations at AUG. This gives rise to a possible interaction of parallel Pfirsch–Schlüter flows in the SOL with the confined plasma. Thus, the different H-mode power threshold for the two drift configurations cannot be explained in the same way at AUG as suggested by LaBombard et al. [Phys. Plasmas 12, 056111 (2005)] for Alcator C-Mod. Finally, comparisons of Er profiles in favorable and unfavorable drift configurations at the respective confinement transitions show that also the Er gradients are all different, which indirectly indicates a different type or strength of the characteristic edge turbulence in the two drift configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effects of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer

The influences of ion trapping and fluctuations of electron temperature and plasma flow on cross-beam energy transfer (CBET) are examined using two- and three-dimensional particle-in-cell simulations in parameter regimes relevant to recent CBET experiments at the OMEGA laser facility. In mid-Z plasma irradiated by an intense pump beam and weaker probe beam, ion trapping, collisional de-trapping, and plasma flow induced by thermal effects are shown to affect the CBET gain. Ion trapping can enhance or detune the CBET resonance [Nguyen et al., Phys. Plasmas 28, 082705 (2021)]. Collisional de-trapping can affect the CBET gain at low seed beam intensity near the onset threshold for ion trapping. Thermal-effects-induced flow can also detune the CBET resonance at a level comparable to that from trapping at low seed beam intensity. As a consequence, the CBET gain is sensitive to collisions and dimensionality at low seed beam intensity where ion trapping is weak but is insensitive to collisions and dimensionality at high seed beam intensity where ion trapping is strong.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Universal Power-law Prescription for Variability from Synthetic Images of Black Hole Accretion Flows

We present a framework for characterizing the spatiotemporal power spectrum of the variability expected from the horizon-scale emission structure around supermassive black holes, and we apply this framework to a library of general relativistic magnetohydrodynamic (GRMHD) simulations and associated general relativistic ray-traced images relevant for Event Horizon Telescope (EHT) observations of Sgr A*. We find that the variability power spectrum is generically a red-noise process in both the temporal and spatial dimensions, with the peak in power occurring on the longest timescales and largest spatial scales. When both the time-averaged source structure and the spatially integrated light-curve variability are removed, the residual power spectrum exhibits a universal broken power-law behavior. On small spatial frequencies, the residual power spectrum rises as the square of the spatial frequency and is proportional to the variance in the centroid of emission. Beyond some peak in variability power, the residual power spectrum falls as that of the time-averaged source structure, which is similar across simulations; this behavior can be naturally explained if the variability arises from a multiplicative random field that has a steeper high-frequency power-law index than that of the time-averaged source structure. We briefly explore the ability of power spectral variability studies to constrain physical parameters relevant for the GRMHD simulations, which can be scaled to provide predictions for black holes in a range of systems in the optically thin regime. We present specific expectations for the behavior of the M87* and Sgr A* accretion flows as observed by the EHT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of particle size and AQDS on the flow of electron equivalents between magnetite and aqueous Fe2+

Magnetite can occur naturally in nano- to micro-size regimes and widely coexists with aqueous Fe2+ (Fe2+ (aq)) in natural environments. However, the effects of magnetite particle size on its interaction with Fe2+ (aq) in anoxic subsurface environments, particularly with redox-active organics, remain unclear. In this study, the interactions of Fe2+ (aq) with magnetite particles of 12 nm versus 109 nm (Mag-12 vs. Mag-109), with/without anthraquinone- 2,6-disulfonate (AQDS), were studied based on equilibrium Fe2+ (aq) concentrations, kinetics of AQDS reduction, and structural versus surface-localized Fe(II)/Fe(III) ratios (xstru and xsurf) of magnetite. In the absence of AQDS, Mag-12 tends to release Fe2+ (aq) at pH 7 but sorb Fe2+ (aq) at pH 8, while Fe2+ (aq) uptake by Mag-109 is observed at both pH 7 and 8. The amounts of Fe2+ (aq) adsorbed per unit area of Mag-109 is higher than that of Mag-12, due to the higher electron-accepting capacity of Mag-109 that facilitates interfacial electron transfer (IET) from surfaceassociated Fe(II) to structural Fe(III). The increases of xstru and xsurf in Mag-109 after reaction with Fe2+ (aq) at pH 7 and 8 suggest Fe2+ (aq) incorporation or electron injection into the structure of Mag-109. The presence of AQDS promotes Fe2+ (aq) uptake by both Mag-12 and Mag-109. However, AQDS reduction by Fe2+-amended Mag-12 results in the decrease of xstru and inhibits Fe2+ (aq) incorporation or electron injection into the structure. On the contrary, the increase of xstru observed in Fe2+-amended Mag-109 after reaction with AQDS suggests that Fe2+ (aq) incorporation or electron injection into the surface structure and then consequently into the interiors is more favorable for magnetite with larger particle sizes. The different flow directions of electron equivalents across the solid-solution interfaces can be attributed to the relatively higher electron-accepting capacity, i.e. redox potential, of Mag-109 than Mag-12; larger particle sizes facilitate IET from surface-associated Fe(II) to structural Fe(III) and promotes further Fe2+ (aq) uptake, culminating in the pronounced changes of redox potentials in magnetitebearing solutions. The results demonstrate that particle size and redox-active organics are important factors to affect reductive activity of Fe2+-magnetite system in redox-oscillating environments.

Peng, Huan↗