SEARCH · Search NASA

Results for “Flow characteristics”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Electrode biasing maintains the edge shear layer at high density in the J-TEXT tokamak

Collapse of the edge flow shear as the line-averaged density approaches the Greenwald density limit has been observed as a precursor to the enhanced edge particle flux characteristic of proximity to the density limit regime. Here, we report the use of a biased electrode to sustain the edge shear layer in high density discharges, in which the shear layer would otherwise collapse. A stable increase in line-averaged density is observed along with a strong increase in edge density. These experiments were carried out on the J-TEXT tokamak. The Reynolds stress at the edge is enhanced, and the zonal flow sustained, while density perturbation levels, the flux of turbulence internal energy (i.e., turbulence spreading), and particle and heat flux all decrease significantly. Electron adiabaticity increases, and bias voltage modulation experiments show that an increase in the edge shear leads the increase in adiabaticity. Furthermore, these results suggest that external edge E × B flow shear drive may be of interest for sustaining edge plasma states at high density, and support the hypothesis that collapse of the edge shear layer triggers the onset of the strong transport and turbulence characteristic of the density limit regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Alteration in micro-mechanical characteristics of wellbore cement fracture surfaces due to fluid exposure

The stiffness of wellbore cement fracture surfaces was measured after exposing to the advective flow of nitrogen, silicone oil, and medium sweet dead crude oil for different exposure periods. The test specimens were extracted from fractured cement cylinders, where the cement fracture surfaces were exposed to the different fluids up to 15 weeks. A nanoindenter with a Berkovich indenter tip was used to measure load-indentation depth data, which was used to extract the elastic modulus (E) and nano-hardness (H) of the cement fracture surfaces. A reduction in the elastic modulus compared with an unexposed specimen were observed in all the specimens. Both elastic modulus and nano-hardness for the specimens exposed to silicone oil were lower than specimens exposed to nitrogen gas and varied with the period of exposure. The elastic modulus and nano-hardness of the specimens exposed to crude oil were the lowest with a significant decrement with the exposure period. The frequency distribution of the nanoindentation measurements shows that the volume-fraction ratio of the two types of cement hydrated nanocomposites for both the unexposed and test specimens is about 70:30%. Phase transformation beneath the indenter is observed for all of the specimens, with more obvious plastic deformation in specimens exposed to crude oil. Analytical measurements (SEM, EDS, FT-IR, and XRD) on exposed cement fracture surfaces reveal different levels of physical and chemical alteration that are consistent with the reduction in stiffness measured by nanoindentation. Furthermore, the study suggests that cement stiffness will decrease due to crude oil exposure, and the fracture will be sensitive to stress and pore pressure with time.

42 ENGINEERING↗

DOE OSTI2025

A 3D helical filament surrogate model for 3D tokamak equilibria

A novel approach for efficient representation of three-dimensional (3D) tokamak equilibria is investigated, where a set of helical current filaments occupying the plasma region are employed to resolve deviations from the two-dimensional (2D) axi-symmetric state. A discrete set of 3D filaments, located at rational surfaces for a given toroidal mode number n and following the 2D equilibrium field lines (thus forming closed current loops), are found to provide a surrogate model of 3D equilibria with reasonable accuracy. Specifically, application of the filament model to 3D perturbed equilibria, due to the resonant magnetic perturbation (RMP) in DIII-D and MAST-U discharges, reveals that (1) a single helical filament per rational surface is sufficient; (2) 21 such helical filaments are capable of representing the n = 2 3D response field in MAST-U with less than 10% relative error as compared to that computed by a full magnetohydrodynamic code; (3) optimizing currents (both amplitude and phase) flowing in 3D filaments with fixed geometry, the highest accuracy fitting is found to depend on the characteristics of the 3D equilibria such as the coil current phasing of the RMP coils in our case studies. Here, whis filament approach is also applicable for generating surrogate models of other type of 3D tokamak equilibria, including those during the initial phase of the plasma disruption.

MARS-F↗

DOE OSTI2022

Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point

Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

A&L Annual Report: GRM-MAPS

Exploding Bridge Wire (EBW) is an important design for initiating detonators relevant to several LLNL systems. In such designs, a metal bridge wire is in direct contact with a porous secondary high explosive (HE), typically pressed to around 50% of its crystalline density. Initiation is believed to occur through shock waves generated when the bridge wire is made to explode via the rapid passing of a strong current pulse. Permeametry is an important characterization tool for the porous powdered HE (PETN, HMX, etc.) component in such devices. Such powders are well known to coarsen with age, with an accompanying deterioration in performance, which includes increased function time and sensitivity. Among various material characteristics, one that has been shown to correlate the most with such performance parameters is the flow-permeable surface area (FSSA), measured using air-permeametry. Commercially available permeametry apparatuses, such as that from Fisher Scientific or a more modern version from Micromeritics suffer from a few drawbacks, most notably as related to monitoring age-related changes in FSSA. This is because it requires re-compaction of the aged powder into the permeametry tube, which in most cases leads to an unwanted increase in FSSA. To address this and a few other shortcomings, we have been developing a modern version of the permeameter, Gee-Reinstein-Maiti Modern Air-Permeametry System (GRM-MAPS). In this report, we summarize progress made in FY25 toward perfecting and calibrating the instrument, which we aim at deploying within surveillance in FY26.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Distinguishing Tidal Disruption Events from Impostors

Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that “TDE signatures” have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter both reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, “micro-TDEs,” and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Testing of Microchannels and Lab-Grown Stress Corrosion Cracks for Quantification of Aerosol Transmission

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using higher backfill pressures in the canister, up to approximately 800 kPa, compared to their horizontal counterparts. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. While the engineered microchannel has dimensions similar to actual SCCs, it does not reproduce the tortuous path the aerosol laden flow would have to traverse for eventual transmission. SCCs can be rapidly grown in a laboratory setting given the right conditions, and initial characterization and clean-flow testing has begun on lab grown crack samples provided to Sandia National Laboratories (SNL). Many such samples are required to produce statistically relevant transmission results, and SNL is developing a procedure to produce samples in welded steel plates. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Experimental and Simulation Analysis of Binary Mixtures of Biomass and Inert Material

A literature review was conducted to identify experiments on various scales of biomass combustion specific to circulating fluidized bed (CFB) combustors. Typical feedstock and inert bed material characteristics for biomass combustion in a CFB were gathered. This information was used to design laboratory-scale, non-reacting experiments on biomass-inert bed material mixtures in a 1-in. cylindrical test rig in the Multiphase Flow Analysis Laboratory (MFAL) in Morgantown, WV. The biomass chosen was a hardwood. Pellets from the hardwood were milled and sieved to obtain three size ranges—small, medium, and large. Three mass loadings of each of the biomass samples were fluidized with glass beads as the inert bed material to study the mixing characteristics under varying fluidizing gas velocities. Hydrodynamic data including pressure drop measurements at select vertical locations and video recordings of the mixture behavior were generated in this effort.

09 BIOMASS FUELS↗

DOE OSTI2022

Mineral scaling on brass and aluminum surfaces with a range of wettability

Crystallization fouling, or mineral scaling, is fundamental to various surfaces in energy and environmental applications. This paper considers, for the first time, a systematic investigation of scaling of nonwetting aluminum and brass surfaces relative to their bare counterparts. Metallic surfaces with a range of nonwetting characteristics from hydrophilic to superhydrophobic and lubricant infused surfaces were fabricated using a facile etching method. Systematic dynamic flow scaling experiments are conducted on the six different surface types at three flow levels and two temperature levels. The interactive effects of the different parameters are analyzed relative to their correspondence to the nucleation theory. The study reveals that nonwetting metallic superhydrophobic or hydrophobic surfaces reduce fouling by over 50% compared to bare counterparts. Further, lubricant infused metallic surfaces demonstrate superior anti-scaling performance by over 60% and up to 90% fouling reduction. The study is the first to present fundamental insights into the engineering and anti-scaling performance of nonwetting brass and aluminum surfaces.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The Impacts of Rotational Mixing on the Precipitation Simulated by a Convection Permitting Model

With increased availability of computational resources, regional and global scale convection-permitting model (CPM, Δx ~ 1–10 km) simulations are becoming more common. CPMs have improved accuracy in their representation of deep convection and mesoscale convective systems (MCSs) compared to coarser resolution models. However, CPMs still exhibit convective cloud and precipitation biases relative to observations, notably a lesser frequency of light precipitation rates and greater frequency of heavy precipitation rates. In this work we hypothesize that these CPM biases are related to under-resolved mixing between convective updrafts and their surrounding environment. To test this hypothesis, we introduce a parameterization to the Weather Research and Forecasting model (WRF) that adds a small angular rotation of the grid-scale flow about the axis perpendicular to the plane of convective drafts. This rotated flow is then allowed to alter advection of moisture and hydrometeors. The effects of such mixing on precipitation characteristics are evaluated in month-long 4-km grid spacing simulations over the Amazon. The enhanced mixing transports moisture and condensate from convective cores to other areas including downdrafts. This increases the frequency of low-precipitable water and light precipitation. It also decreases the frequency of intense precipitation from isolated deep convection and MCSs, increases cloud top temperatures, reduces radar echo-top heights, and increases overall precipitation by altering the relationship of precipitation with precipitable water, in better agreement with observations. The results suggest when optimized using multiple observations, such an approach may provide a path toward more accurate representation of convection and precipitation statistics in convection-permitting simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Supporting data for "Frozen flume experiments indicate rapid permafrost riverbank erosion depends on bank roughness"

This dataset contains supporting files detailing five frozen flume experiments conducted at the Caltech Earth Surface Dynamics Laboratory to investigate rates of ablation-limited permafrost riverbank erosion under controlled conditions. Water flowed past a bank of saturated, frozen sand and ice and gradually eroded the bank by thawing pore ice and immediately entraining sand and washing it downstream. Experiments were scaled for flow hydraulics and heat transfer allowing comparisons between our results and natural permafrost riverbanks. For each experiment, we measured the initial and final sand bank topography using a Keyence laser scanner, water surface slope at 3-min intervals throughout the experiment using a Massa sonar scanner, bank erosion using 10-sec overhead timelapse imagery taken by an overhead camera, water and bank temperature using thermistors frozen into the sand bank and sampling at 2 Hz, and water discharge using an in-line flow meter. We include calibration data for the carriage (engineered by the Saint Anthony Falls Laboratory) used to make sonar and laser topography measurements. We also include calibration data for temperature sensors, water discharge measurements, and images of a regular grid placed in the flume to align overhead camera images with the carriage datum. Grain size analysis for the channel bed (gravel) was produced using a pebble count and bank sand was measured using a Camsizer X2. In addition to the five frozen experiments, we include sonar scans of water surface slope and Keyence scans of bed and bank topography for calibration experiments ran with an immobile gravel bank and bed.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗

DOE OSTI2020

Towards understanding the development and characteristics of under-expanded flash boiling jets

Here in this research, the n-hexane jets discharged from a straight channel were studied numerically to reveal the development and characteristics of under-expanded flashing jets. Firstly, the development of flash boiling jets was investigated and compared to that of two-phase jets without phase transition. It was shown that the two-phase jets without phase transition can be under-expanded only if the flow is choked. In contrast, flash boiling jets can be under-expanded whether or not the choked conditions are verified. Secondly, the factors affecting the flashing jet behaviors were discussed. The relaxation time in phase transition, the expansion ratio, and the injection pressure have been identified as three key factors determining the jet characteristics. With high enough expansion ratio and low enough relaxation time, the expansion was so strong that the flow became transonic and shock structures similar to those of under-expanded gaseous jets appeared near the orifice exit. Finally, two types of choking have been identified for flash boiling jets, i.e., vaporization-induced choking and expansion-induced choking. The former one is mainly ascribed to the vapor formation in the orifice which reduces the local speed of sound, while the latter one is due to the expansion of under-expanded two-phase flow and the resultant flow acceleration.

42 ENGINEERING↗

DOE OSTI2023

Two-dimensional glass transition–like behavior of Janus particle–laden interface

Abstract Understanding the interactive behavior of Janus particles (JPs) is a growing field of research. The enhancement in binding energy, in comparison to homogenous particles, and the dual characteristic of JPs open up new possibilities for novel applications. In many such applications, interfacial materials become subjected to flows that produce dilational and shear stresses. Therefore, it is important to understand the impact that the Janus character brings to interfaces. In this work, we study the microstructure of two-dimensional (2D) JP monolayers formed at the air–water interface and examine the shear viscoelasticity with an interface rheometer that was adapted for in situ surface pressure control via a Langmuir trough. We extend concepts from bulk rheology to data obtained from interfacial rheology as a tool to understand and predict the monolayer’s viscoelastic behavior. Finally, by calculating the time relaxation spectrum from the measured 2D dynamic moduli, we conclude that a phenomenon similar to glass transition is taking place by analogy.

42 ENGINEERING↗

DOE OSTI2026

Boundary-layer receptivity to oblique freestream vorticity waves for a high-enthalpy hypersonic flow

The receptivity of a Mach 15 straight-cone boundary layer to oblique freestream vorticity waves is investigated using direct numerical simulation (DNS) alongside linear stability theory and the linear parabolized stability equations. A thermochemical nonequilibrium gas model is used. Oblique freestream vorticity waves at frequencies of 400, 800, and 1200 kHz are considered, with incident angles ranging from 0° to 29.4° at 400 kHz, 0° to 15.7° at 800 kHz, and 0° to 20.6° at 1200 kHz. The 400 kHz case is of primary interest due to the strong second-mode amplification at this frequency. Although the underlying base flow is axisymmetric, the oblique vorticity waves lead to a fully three-dimensional boundary-layer disturbance whose characteristics vary depending on the azimuthal ray relative to the freestream wave. Moving downstream within the second-mode instability region, some clear trends emerge in terms of the boundary-layer disturbance amplitudes; that is, disturbance amplitudes are highest at the leeward ray (relative to the freestream wave), but weakest about halfway between the windward and leeward rays. Increasing the incident angle causes the amplitudes to increase on the leeward ray and decrease on the windward ray. Moreover, the boundary-layer disturbance throughout contains a wide spectrum of azimuthal wavenumbers in which the disturbance energy falls off at higher wavenumbers. Increasing the incident angle causes the azimuthal spectrum of the boundary-layer disturbance to broaden overall. Qualitatively similar results are found for the two higher frequencies leading up to the peak-amplitude locations corresponding to the second-mode instability.

42 ENGINEERING↗

DOE OSTI2024

Solvent-Dependent Dynamics of Cellulose Nanocrystals in Process-Relevant Flow Fields

Flow-assisted alignment of anisotropic nanoparticles is a promising route for the bottom-up assembly of advanced materials with tunable properties. While aligning processes could be optimized by controlling factors such as solvent viscosity, flow deformation, and the structure of the particles themselves, it is necessary to understand the relationship between these factors and their effect on the final orientation. In this study, we investigated the flow of surface-charged cellulose nanocrystals (CNCs) with the shape of a rigid rod dispersed in water and propylene glycol (PG) in an isotropic tactoid state. In situ scanning small-angle X-ray scattering (SAXS) and rheo-optical flow-stop experiments were used to quantify the dynamics, orientation, and structure of the assigned system at the nanometer scale. The effects of both shear and extensional flow fields were revealed in a single experiment by using a flow-focusing channel geometry, which was used as a model flow for nanomaterial assembly. Due to the higher solvent viscosity, CNCs in PG showed much slower Brownian dynamics than CNCs in water and thus could be aligned at lower deformation rates. Moreover, CNCs in PG also formed a characteristic tactoid structure but with less ordering than CNCs in water owing to weaker electrostatic interactions. The results indicate that CNCs in water stay assembled in the mesoscale structure at moderate deformation rates but are broken up at higher flow rates, enhancing rotary diffusion and leading to lower overall alignment. Albeit being a study of cellulose nanoparticles, the fundamental interplay between imposed flow fields, Brownian motion, and electrostatic interactions likely apply to many other anisotropic colloidal systems.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multiple experimental studies of pore structure and mineral grain sizes of the Woodford shale in southern Oklahoma, USA

Pore structure study is an important part of unconventional shale reservoir characterization, since the pore system provides the primary petroleum storage space and fluid flow pathways. Previous studies have suggested that the pore structure is related to the total organic carbon (TOC) content, mineral compositions, and the maturity of the organic matter (OM). However, few studies have focused on the mineral grains, the primary grains being deposited but before cementation, which are the building blocks of shale. Eight Woodford Shale outcrop samples from southern Oklahoma were chosen to study the effects of mineral grain size on the pore structure characterization, using multiple and complementary experimental approaches, including laser diffraction, mineralogy, TOC, pyrolysis, liquid immersion porosimetry, mercury intrusion porosimetry, gas physisorption, (ultra) small angle X-ray scattering, scanning electron microscopy, and spontaneous imbibition. The results from different experiments of eight samples show that the Woodford Shale has the mean mineral grain diameters at 3–6 μm, a wide range of porosity at 3–40% and pore diameters at 50–1,000 nm, and various pore connectivity. Grain size variation was probably caused by the sea-level fluctuation during its deposition, which affect the porosity, pore size distribution, and pore connectivity. With decreasing mineral grain sizes, the porosity tends to increase while the pore connectivity worsens. The results also indicate that OM and carbonates in this low-maturity Woodford Shale could block the pores and decrease the porosity. Coupling with the grain size analyses, the control of depositional environment on grain sizes and subsequent effects on pore structure is identified. The pore structure characteristics over a wide pore-diameter range provided by multiple experiments could improve the understanding of storage space and fluid flow in the Woodford Shale to further increase its petroleum production.

54 ENVIRONMENTAL SCIENCES↗