SEARCH · Search NASA

Results for “SATURATION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Collisionless conduction in a high-beta plasma: a collision operator for whistler turbulence

The regulation of electron heat transport in high-β, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as β$^{-1}_{e}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of β$_e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected β$^{-1}_{e}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known 90° scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗

DOE OSTI2025

Mitigating CaCO 3 crystal nucleation and growth through continuous ion displacement via alternating electric fields

Mineral crystal formation poses a challenge on surfaces (e.g., heat exchangers, pipes, membranes, etc.) in contact with super-saturated fluids. Applying alternating currents (AC) to such surfaces can prevent surface crystallization under certain conditions. Here, we demonstrate that ion displacement induced by periodic charging and discharging of the electrical double layer (EDL) inhibits both heterogeneous and homogeneous nucleation (and crystal growth) of CaCO 3 . Titanium sheets (meant to simulate metallic heat exchanger surfaces) are immersed in super-saturated CaCO 3 solutions with a saturation index >11. We show that at relatively high AC frequencies, incomplete EDL formation leads to an alternating electric field that propagates far into the bulk solution, inducing rapid ion migration that overwhelms the Brownian motion of ions. Electrochemical characterization reveals EDL charging/discharging under AC conditions that greatly inhibits precipitation. Operating at 4 V pp , 0.1–10 Hz reduces turbidity by over 96% and reduces CaCO 3 coverage on the metal plates by over 92%. Based on electrokinetic and crystallization models, the ion displacement velocity (exceeding the mean Brownian velocity) and displacement length disrupts ion collision and crystal nucleation. Overall, the technique has potential for preventing mineral crystal formation in heat exchangers and many other industrially relevant systems.

42 ENGINEERING↗

DOE OSTI2024

Belowground plant allocation regulates rice methane emissions from degraded peat soils

Carbon-rich peat soils have been drained and used extensively for agriculture throughout human history, leading to significant losses of their soil carbon. One solution for rewetting degraded peat is wet crop cultivation. Crops such as rice, which can grow in water-saturated conditions, could enable agricultural production to be maintained whilst reducing CO 2 and N 2 O emissions from peat. However, wet rice cultivation can release considerable methane (CH 4 ). Water table and soil management strategies may enhance rice yield and minimize CH 4 emissions, but they also influence plant biomass allocation strategies. It remains unclear how water and soil management influences rice allocation strategies and how changing plant allocation and associated traits, particularly belowground, influence CH 4 -related processes. We examined belowground biomass (BGB), aboveground biomass (AGB), belowground:aboveground ratio (BGB:ABG), and a range of root traits (root length, root diameter, root volume, root area, and specific root length) under different soil and water treatments; and evaluated plant trait linkages to CH 4 . Rice (Oryza sativa L.) was grown for six months in field mesocosms under high (saturated) or low water table treatments, and in either degraded peat soil or degraded peat covered with mineral soil. We found that BGB and BGB:AGB were lowest in water saturated conditions where mineral soil had been added to the peat, and highest in low-water table peat soils. Furthermore, CH 4 and BGB were positively related, with BGB explaining 60% of the variation in CH 4 but only under low water table conditions. Our results suggest that a mix of low water table and mineral soil addition could minimize belowground plant allocation in rice, which could further lower CH 4 likely because root-derived carbon is a key substrate for methanogenesis. Minimizing root allocation, in conjunction with water and soil management, could be explored as a strategy for lowering CH 4 emissions from wet rice cultivation in degraded peatlands.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

An octahedral Mach B-dot probe for 3D flows and magnetic fields in the edge of reversed field pinches

Measurements and simulations show that plasma relaxation processes in the reversed field pinch drive and redistribute both magnetic flux and momentum. To examine this relaxation process, a new 3D Mach B-dot probe has been constructed. This probe collects ion saturation currents through six molybdenum electrodes arranged on the flattened vertices of an octahedron made of boron nitride (BN). The ion saturation current flows through configurable voltage dividers for measurement and returns through one of six selectable return electrodes equally spaced along the 12 cm BN probe arm. In addition, the probe arm houses three B-dot magnetic pickup coils in the BN stalk immediately below to the octahedron, to measure the local magnetic field. Inserted in the Madison Symmetric Torus (MST) during deuterium discharges with 220 kA plasma current, density of 0.8 × 10 13 cm –3 , the probe collects ion saturation currents with sawtooth-like peaks correlated with relaxation events. This compact octahedral design fitting six Mach electrode surfaces within a 1 cm3 cube will enable future multi-point, multi-field probes compatible with the 1.5 in. ports of MST. Such probes will allow for flow circulation, current, and canonical vorticity to be calculated in the center of the finite difference stencil formed by the measurement locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

General kinetic ion-induced electron emission model for metallic walls applied to biased Z-pinch electrodes

A kinetic ion-induced electron emission (IIEE) model for general applications is developed to obtain the emitted electron energy spectrum for a distribution of ion impacts on a metallic surface. We assume an ionization cascade mechanism and use empirical models for the ion and electron stopping powers. The emission spectrum and the secondary electron yield (SEY) are validated for a variety of materials. The IIEE model is used to study the effect of IIEE on the plasma-material interactions of Z-pinch electrodes. Un-magnetized Boltzmann-Poisson simulations are performed for a Z-pinch plasma doubly bounded by two biased copper electrodes with and without IIEE at bias potentials from 0 to 9 kV. At the anode, the SEY decreases from 0 to 1 kV, but then increases at higher bias potentials. At the cathode, the SEY is much larger due to higher energy ion bombardment and grows with bias potential. As the bias potential increases, the emitted cathode electrons are accelerated to higher energies into the domain, collisionally heating the plasma. Above 1 kV, the heating is strong enough to increase the plasma potential. Despite SEY greater than 1, only a classical sheath forms as opposed to a space-charge limited or inverse sheath due to the emitted electron flux not reaching the space charge current saturation limits. Furthermore, the current in the emissionless cases saturates to a value lower than experiment. With IIEE, the current does not saturate and continues to increase with the 4 kV case, matching most closely with the experiment.

Carbon based materials↗

DOE OSTI2024

Examining transport and integrated modeling predictive capabilities for negative-triangularity scenarios

This paper investigates the predictive capabilities of TGYRO and TGLF models in assessing the performance of negative triangularity (NT) plasmas compared to positive triangularity (PT) plasmas in fusion devices. TGYRO predicts kinetic profiles, while TGLF analyzes turbulent transport. The study reveals that TGYRO reasonably predicts NT profiles similar to PT, although it overpredicts the high-power scenarios where there is increased experimental MHD activity. TGLF analysis finds reduced linear growth rates in NT and altered flux spectra relative to PT. Additionally, the TGLF SAT0 saturation model is observed to predict high-k transport and a reduction of particle transport with the electron temperature gradient. These findings are further corroborated by core-pedestal modeling using the Stability Transport Equilibrium Pedestal workflow, showing stronger confinement improvements in NT, particularly at higher power densities for the SAT0 saturation model. Furthermore, the study underscores the importance of accurately capturing turbulence saturation mechanisms for NT in order to project its performance accurately in fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗

DOE OSTI2025

Disentangling core and edge mechanisms of the density limit in DIII-D negative triangularity plasmas

The density limit is investigated in the DIII-D negative triangularity plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans show distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered by the onset of a non-disruptive, high-field side radiation front and the associated cooling, which clamps the edge density below . In contrast, the core density continues to rise until it saturates, a state characterized by enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis ≈ 6), increased Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative cooling and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

density limit↗

DOE OSTI2024

Out-of-plane magnetic phase diagram of the Kitaev quantum magnet Na 2 Co 2 TeO 6

We have investigated the magnetic properties and mapped out the phase diagram of the honeycomb magnet Na 2⁢ Co 2 ⁢TeO 6 with Co 3d 7 in out-of-plane magnetic fields. This material has previously been proposed to show nearest-neighbor Kitaev interactions between Co spins and maybe even Kitaev quantum spin liquid behavior in high fields. At 0.01 T, we observe a thermal phase transition at T N = 27K, transitioning from a paramagnetic state to a canonical ferrimagnetic state. Upon increasing the magnetic field, a spin-floplike phase transition occurred before saturation of J = 1/2 between 10 K and T N . Below 10 K, a peak-dip-peak structure emerges between 10 and 17 T in the magnetic susceptibility (d⁢M/⁢dH) before the magnetic saturation, reminiscent of magnetization plateau behavior. The measurement of the magnetocaloric effect also shows dip-peak-dip behavior in this field range. Our data can be explained by an XXZ model with a single ion anisotropy and possibly small Kitaev and Γ exchange interactions. Finally, we also determined the magnetization saturation field that helps constrain the energy scale of the exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Transverse energy-energy correlators in the color-glass condensate at the electron-ion collider

We investigate the transverse energy-energy correlators (TEEC) in the small- x regime at the upcoming Electron-Ion Collider (EIC). Focusing on the back-to-back production of electron-hadron pairs in both e p and e A collisions, we establish a factorization formula given in terms of the hard function, quark distributions, soft functions, and TEEC jet functions, where the gluon saturation effect is incorporated. Numerical results for TEEC in both e p and e A collisions are presented, together with the nuclear modification factor R A . Our analysis reveals that TEEC observables in deep inelastic scattering provide a valuable approach for probing gluon saturation phenomena. Our findings underscore the significance of measuring TEEC at the EIC, emphasizing its efficacy in advancing our understanding of gluon saturation and nuclear modifications in high-energy collisions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗

DOE OSTI2025

Quantifications of growth rate and particle size of magnesium silicate hydrate cement using small‐angle X‐ray scattering

The focus of recent research on cement binders has shifted to alternatives such as magnesium silicate hydrate (M-S-H) to mitigate the environmental impact of conventional concrete production. However, a detailed understanding of the morphology and growth rate of M-S-H, which informs its applicability as a cementing agent, is lacking. These properties are important as they influence cement hydration and consequently the development of macroscopic properties. Here, in this study, the particle size distribution and growth rate of M-S-H synthesized by mixing aqueous stock solutions were analyzed using synchrotron-based small-angle X-ray scattering. The reaction solutions contained dissolved magnesium and silicon with concentrations ranging from 2 to 500 mM, corresponding to saturation indices with respect to M-S-H endmembers, M 0.75 -S-H and M 1.50 -S-H, between 5.8 and 10.3. The precipitation rates of M-S-H align with affinity-based kinetic models and can be described by the following equations: For M 0.75 -S-H: 𝑟 M 0.75 −S−H =10 −9.54±0.70 ⁢(Ω$^{0.59±0.08}_{M_{0.75}⁢\textrm{SH}}$−1), and for M 1.50 -S-H: 𝑟 M 1.50 −S−H =10 −7.56±0.37 ⁢(Ω$^{0.34±0.04}_{M_{1.50}⁢\textrm{SH}}$−1), where 𝑟 is the rate in mol/L/s, and Ω represents the saturation ratio with respect to M-S-H. Moreover, we observed that changing the [Mg]/[Si] ratio affects particle size, a trend that cannot be fully explained by changes in the saturation index alone. These insights are critical for improving the predictability and control of cementation in concrete that utilizes M-S-H-based binders.

M-S-H↗

DOE OSTI2025

Experimental Investigation of Buoyant Flow in Realistic Bedforms With Heterogeneous Wettability

Submeter-scale geologic heterogeneity greatly affects CO 2 plume migration and retention. In this work, we present meter-scale laboratory experiments that can capture the impact of realistic submeter-scale geologic heterogeneity on multiphase flow and trapping. We produce realistic sedimentary formations consisting of ripple deposits with varying grain size contrast and wettability in a meter-scale slab chamber. Then, we conduct multiphase flow experiments with analog fluids through these structures and measure the saturation patterns, capillary heterogeneity trapping (CHT), and overall trapping performance. When we alter the ripple bedform architecture, variations in trapped saturation and CHT (10–20%) increment are exhibited. Similar growth in trapping performance is also observed when grain size contrast increases. Finally, wettability changes (water- to oil-wet) can increase nonwetting saturation and CHT up to 5% and 10–20%, respectively. These results emphasize the importance of correctly characterizing the impact of small-scale heterogeneities and wettability changes. We believe this is the first time that multiphase flow experiments were conducted in meter-scale domains with realistic ripple bedforms and heterogeneous wettability to investigate plume migration and trapping.

58 GEOSCIENCES↗

DOE OSTI2024

UT-GOM2-2 Preliminary Report: Terrebonne Basin Northern Gulf of Mexico, 30 July-28 September 2023

In the summer and fall of 2023, the Gulf of Mexico Deepwater Hydrate Coring Expedition (UT-GOM2-2) drilled, cored, made downhole measurements, and analyzed samples from the seafloor to the base of the gas hydrate stability zone in one location (Site H, WR313) in the Terrebonne basin, deepwater Gulf of Mexico. Analyses of data and samples from the expedition will inform biological, geochemical, and geomechanical models to constrain the role of gas hydrates in the carbon cycle and the potential for gas hydrates as an energy resource. Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were obtained. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of ~200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The biogenic ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years) with a pronounced increase in sedimentation rate with depth. Beneath 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C 1 /C 2 ) and the methane to ethane plus propane (C 1 /(C 2 +C 3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ 13 C isotopic signature of methane ranges between -69.9 and -78.5 ‰ Vienna Pee Dee Belemnite (VPDB). Pressure core recovery of all sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

DOE OSTI2025

Using Separation-Enhanced Isotope Ratio Mass Spectrometry to Enable Increased Renewable Carbon Content in Transportation Fuels (CRADA 525)

Stable isotope ratio measurements of carbon atoms using isotope ratio mass spectrometry (IRMS) can be an effective tool for quantifying biogenic carbon in co-processed fuels, with results approaching the precision and accuracy of accelerator mass spectrometry (AMS). The lower cost of an IRMS may enable deployment to refineries, improving access and analysis turnaround times (≤2 hours), and, by extension, provide data that can allow process optimization to maximize renewable carbon in desired refinery products. This project explored the integration of chemical separation with IRMS analyses to enable highly detailed tracking of biogenic carbon into fuel product streams separated by boiling point range, chemical class, or specific compound. Forty-nine fuels and fuel components of fossil and biogenic origin, spanning gasoline and diesel boiling point ranges, were received from three refiners and were analyzed for their δ 13 C values via IRMS. Results spanned a 13 C range from ca. 10‰ to 44‰ and reflect materials derived from sustainable sources (e.g., C4 or C3 plants, animal-based pathways, syngas) or from fossil-derived fuels. Common ranges are approximately 18‰ to 9‰ and approximately 30‰ to 20‰ for C4 and C3 plants, respectively, and approximately 34‰ to 24‰ and approximately 70‰ to 33‰ for petroleum-derived fuels and methane, respectively. Fuel-like standards were developed and tested using direct-injection elemental analyzer (EA) IRMS for liquid fuels. This method was compared with the published methods, yielding statistically similar results. Four blend curve sets were produced ranging from 0% to 100% of a fuel containing biogenic carbon, focusing on 0% to 10% biogenic carbon. Linear fits were the most applicable for two of the four blend curve sets; however, two sets were found to exhibit slightly quadratic behavior, which was more pronounced in low biogenic blend samples, necessitating second-order fits. The origin of the slight quadratic behavior remains unclear; however, the discussion points to possible interpretations. CanmetENERGY thoroughly characterized a majority of the samples using one- and two-dimensional gas chromatography (GC and GC×GC, respectively) and other analyses. Selected samples were subjected to solid phase extraction (SPE) for saturate, olefin, aromatic, and polar (SOAP) analysis, and the resulting solvent-diluted fractions containing saturates and aromatics were returned to Pacific Northwest National Laboratory (PNNL), where the solvent was removed via evaporation or physical separation using GC techniques. Characterization and separations provided an understanding of saturate and aromatic content, as well as boiling point ranges for each sample and sample fraction. Samples resulting from SPE were examined using EA-IRMS and gas chromatography combustion IRMS (GC-C-IRMS) analyses. Both approaches suggest that the range in values between end-members can be increased by selecting the paraffinic or aromatic fraction of the end-member or by selecting among individual compounds resulting from GC separation of the paraffinic fractions. Considerable work remains to put these approaches into practice and statistically validate the benefit for using a fraction or individual compound over bulk analysis of a sample. However, initial results suggest that separations provide advantages for samples having blend ratios of less than 10% biogenic blendstocks. 13 C results showed statistically similar biofuel blend results to those obtained at PNNL, although additional work is needed to obtain better reproducibility. Select samples were sent to Los Alamos National Laboratory (LANL) for IRMS measurements and Beta Analytics for AMS measurements. This work suggests that IRMS and AMS yield closely comparable results and in some circumstances, IRMS could serve as a surrogate for AMS. While additional work is needed to better resolve statistical advantages for separations and better show the comparable nature of IRMS and AMS in both the biogenic carbon analysis of bulk chemical classes, initial results from this study suggest that these should be pursued in order to proliferate this approach for quantifying biogenic carbon in transportation fuels to the refinery level, thereby potentially enabling process optimization in co-processing scenarios.

09 BIOMASS FUELS↗

DOE OSTI2025

Fluid/Fluid Interfacial Areas Measured for Different Non-Wetting/Wetting Fluid Pairs in Natural Porous Media

This study examined the impact of fluid type and grain diameter on the interfacial area between different pairs of non-wetting and wetting fluids in natural porous media. Synchrotron X-ray microtomography was used to obtain high-resolution, three-dimensional images of multi-phase porous media systems. Multiple porous media, comprising a range of physical and geochemical properties, were used in this study. The four pairs of non-wetting/wetting fluids used were dense OIL/water, light OIL/water, air/dense OIL, and air/water. Images were obtained over a broad range of wetting phase saturation and for both wetting phase drainage and imbibition conditions. The results showed that for each fluid pair, the total (capillary + film) interfacial area increased with decreasing wetting fluid saturation. Interfacial areas were similar among all fluid pairs for a given porous medium. They were also similar for drainage and imbibition conditions. The maximum specific interfacial area (A m ) was shown to correlate well with inverse median grain diameter. The physical properties of the porous medium appear to have a greater influence on the magnitude of specific total interfacial area for a given saturation than fluid properties or wetting phase history.

NAPL↗