Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heterogeneous media”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2023

Synthesis of PHI-type zeolite harmotome using barium as an inorganic structure-directing agent inspired from natural mineral compositions

Zeolites are nanoporous aluminosilicates widely used as commercial adsorbents, heterogeneous catalysts, and ion-exchange materials due to their unique porosity, acidity, and thermal stability. The vast majority of synthetically-realized zeolite structures are prepared under hydrothermal conditions in alkaline media, and often in the presence of an organic structure-directing agent that facilitates the crystallization of diverse porous networks. Due to economic and environmental disadvantages of organic-based syntheses, it is often desirable to produce zeolites in organic-free media using alkali metal ions as the most commonly employed inorganic structure-directing agents (ISDAs). Here, in this study, we were inspired by the elemental compositions of natural zeolite minerals to examine the use of an alkaline earth metal ion as the ISDA to prepare zeolites with the PHI framework type. Two natural minerals exhibiting this crystal structure are phillipsite and harmotome. Both structures contain sodium and potassium as extra-framework cations, but only harmotome contains barium cations. Here we examined how the combination of alkali and alkaline earth ISDAs promote the formation of PHI-type zeolites. Our findings reveal that each isostructure is prepared under different synthesis conditions where simply adding barium to a growth mixture does not always lead to harmotome formation. Harmotome is formed in more dilute and alkaline growth media; however, time-resolved analyses of zeolite crystallization reveal that harmotome can be synthesized at low temperatures (ca. 65 °C) with rates that are 5-times faster than syntheses of phillipsite. This seemingly indicates that barium reduces the energetic barrier for PHI-type zeolite nucleation in ways that are not fully understood. Moreover, harmotome crystals exhibit an average particle size of ca. 60 nm, which is nearly 40-fold smaller than phillipsite crystals prepared under similar conditions. Collectively, these findings identify routes to control crystallization of PHI-type zeolites and may provide a broader strategy for using the composition of natural zeolite minerals to optimize the preparation of synthetic analogues.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An aggregation-based nonlinear multigrid solver for two-phase flow and transport in porous media

A nonlinear multigrid solver for two-phase flow and transport in a mixed fractional-flow velocity-pressure-saturation formulation is proposed. The solver, which is under the framework of the full approximation scheme (FAS), extends our previous work on nonlinear multigrid for heterogeneous diffusion problems. The coarse spaces in the multigrid hierarchy are constructed by first aggregating degrees of freedom, and then solving some local flow problems. The mixed formulation and the choice of coarse spaces allow us to assemble the coarse problems without visiting finer levels during the solving phase, which is crucial for the scalability of multigrid methods. Specifically, a natural generalization of the upwind flux can be evaluated directly on coarse levels using the precomputed coarse flux basis vectors. We find that the resulting solver is applicable to problems discretized on general unstructured grids. The performance of the proposed nonlinear multigrid solver in comparison with the standard single level Newton's method is demonstrated through challenging numerical examples. It is observed that the proposed solver is robust for highly nonlinear problems and clearly outperforms Newton's method in the case of high Courant-Friedrichs-Lewy (CFL) numbers.

58 GEOSCIENCES↗
DOE OSTI2024

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗
DOE OSTI2024

4D-STEM Mapping of Nanocrystal Reaction Dynamics and Heterogeneity in a Graphene Liquid Cell

Chemical reaction kinetics at the nanoscale are intertwined with heterogeneity in structure and composition. However, mapping such heterogeneity in a liquid environment is extremely challenging. Here, in this work, we integrate graphene liquid cell (GLC) transmission electron microscopy and four-dimensional scanning transmission electron microscopy to image the etching dynamics of gold nanorods in the reaction media. Critical to our experiment is the small liquid thickness in a GLC that allows the collection of high-quality electron diffraction patterns at low dose conditions. Machine learning-based data-mining of the diffraction patterns maps the three-dimensional nanocrystal orientation, groups spatial domains of various species in the GLC, and identifies newly generated nanocrystallites during reaction, offering a comprehensive understanding on the reaction mechanism inside a nanoenvironment. This work opens opportunities in probing the interplay of structural properties such as phase and strain with solution-phase reaction dynamics, which is important for applications in catalysis, energy storage, and self-assembly.

four-dimensional scanning transmission electron mi↗
DOE OSTI2020

Thermodynamic and Kinetic Transitions of Liquids in Nanoconfinement

Core principles of chemistry are ubiquitously invoked to shed light on the nature of molecular level interactions in nanoconfined fluids, which play a pivotal role in a wide range of processes in geochemistry, biology, and engineering. A detailed understanding of the physicochemical processes involved in the flow, structural transitions, and freezing or melting behavior of fluids confined within nanometer-sized pores of solid materials is thus of enormous importance for both basic research and technological applications. This Account provides a perspective on new insights into the thermodynamic and kinetic transitions of nanoconfined fluids in their stable and metastable forms. After briefly introducing the unique properties of mesoporous silicas from the SBA, MCM, and FDU families that serve as the confinement matrices, combining highly ordered single and bimodal mesopore architectures with tunable pore sizes in the ~2-15 nm range and narrow size distributions, recent studies on melting/freezing behavior of water confined in these host matrices are reviewed. While differential scanning calorimetry (DSC) reveals a linear relationship between melting point depression and pore size (independent of the pore shape), as predicted by the Gibbs-Thomson relation, variable temperature 2 H wide-line nuclear magnetic resonance (NMR) spectroscopy studies confirm the core-shell model of water and give evidence for a layer-by-layer freezing mechanism, which gives rise to an apparent fragile-to-strong transition in the solidification dynamics. In contrast to the freezing/melting behavior of water, the effect of nanoconfinement on the glass transition of supercooled liquids is nonuniversal and the glass transition temperature T g can either increase or decrease with the dimensionality and extent of confinement. This nonuniversal behavior is exemplified by the two glass-forming molecular liquids, glycerol and ortho-terphenyl (OTP). While glycerol shows an increase in T g and a pronounced slowdown of the rotational dynamics of the constituent molecules due to a change in the molecular packing between the bulk and the confined liquid, OTP displays a linear and confining-media-dependent depression of T g with increased confinement that is strongly influenced by the pore-liquid interface characteristics. This Account concludes with a focus on recent experimental evidence of extreme spatial and dynamical heterogeneity in both freezing and glass transition processes. This discovery was enabled by the unique mesoporous structures of SBA-16 and FDU-5, possessing bimodal architectures with two interconnected pore types of different size and shape (spherical and cylindrical). For the very first time, two melting points for water and two glass transitions for supercooled OTP, corresponding to a specific pore type, were observed. Collectively, these observations strongly suggest a close mechanistic connection between the local fluctuations in the structure and dynamics of nanoconfined liquids. Finally, while the findings reviewed in this Account provide new insights into thermodynamic and kinetic transitions of fluids, there remain many unanswered questions regarding the effects of nanoconfinement on the fundamental properties of fluids, which offer exciting future opportunities in chemical research.

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

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
DOE OSTI2022

Macrovoid resolved simulations of transport through HPRO relevant membrane geometries

Modeling the transport properties such as diffusivity and permeability of high pressure reverse osmosis (HPRO) membranes is critical for the selection and manufacture of membranes suitable for operation under high pressure. These properties can be significantly affected by the changes in heterogeneous pore structures due to compaction. Here, the modeling platform presented in this work resolves the two scale porosity in HPRO relevant membranes. A synthetic membrane geometry is constructed based on available experimental visualizations of pore structures. The simulations directly capture the flow channeling that results from a combination of material properties and geometric features of the macrovoids. A parametric study is presented to account for the transition of flow characteristics from a material governed regime to a macrovoid governed regime. Permeability of the membrane is evaluated using the simulation data and compared with an existing model that scales with the square of tortuosity over a range of material properties. The model is found to perform well under a narrow range of tortuosity while deviating from the calculated permeabilities at extreme conditions. The effective permeability of the membrane is found to vary by at least two orders of magnitude between the two flow regimes. It is also observed that tortuosity is a bounded property with its upper limit determined by the macrovoid geometry. Consequently, the tortuosity based correlations fail near a flow regime that is mainly governed by the macrovoids. The modeled permeability can be more than an order of magnitude smaller than the simulation result. A new model based on flux-weighted porosity of a membrane is introduced and its correlation with tortuosity is studied. The model agrees with the simulated data as, in addition to tortuosity, it also accounts for the flux partition within and outside the flow channels. Such correlations enable extending the existing understanding of flow characteristics to enhance predictability of porous media models.

42 ENGINEERING↗
DOE OSTI2022

Media for culturing epithelial cells

A method of growing primary human epithelial cells, in particular human epithelial cells using a basal formula containing individual (a) amino acids, (b) vitamins, (c) trace elements, and (d) other organics such as linoleic acid. The basal medium may be a mixture of amino acids, vitamins, and salts that constitute the basic media that is used to culture epithelial cells over a number of population doublings, e.g., over at least one week, while maintaining a normal phenotype and exerting low stress on the cultured cells, and maintaining lineage heterogeneity.

LaBarge, Mark A.↗
DOE OSTI2020

An inelastic homogenization framework for layered materials with planes of weakness

Many geologic materials have a composite structure, in which macroscopic mechanical behavior is determined by the properties, shape, and heterogeneous distribution of individual constituents. In particular, sedimentary rocks commonly exhibit a layered microstructure, with distinct bedding planes that can also form planes of weakness. In this work, we present a homogenization framework for modeling inelastic layered media. The proposed constitutive model allows for distinct micro-constitutive laws for each layer, explicit representation of layer distributions, as well as incorporation of imperfect bonding at the interface between adjacent layers. No a priori assumptions are needed regarding the specific constitutive models used for the layers and interfaces, providing significant modeling flexibility. The overall framework provides a simple and physically-motivated way of defining anisotropic material behavior as an emergent property of the layered microstructure. The model is calibrated using triaxial and true-triaxial experimental data to demonstrate its ability to describe anisotropic deformation and multiple modes of failure.

58 GEOSCIENCES↗
DOE OSTI2023

Model-Based Reconstruction for Multi-Frequency Collimated Beam Ultrasound Systems

Collimated beam ultrasound systems are a technology for imaging inside multi-layered structures such as geothermal wells. These systems work by using a collimated narrow-band ultrasound transmitter that can penetrate through multiple layers of heterogeneous material. A series of measurements can then be made at multiple transmit frequencies. However, commonly used reconstruction algorithms such as Synthetic Aperture Focusing Technique (SAFT) tend to produce poor quality reconstructions for these systems both because they do not model collimated beam systems and they do not jointly reconstruct the multiple frequencies. Here, in this article, we propose a multi-frequency ultrasound model-based iterative reconstruction (UMBIR) algorithm designed for multi-frequency collimated beam ultrasound systems. The combined system targets reflective imaging of heterogeneous, multi-layered structures. For each transmitted frequency band, we introduce a physics-based forward model to accurately account for the propagation of the collimated narrow-band ultrasonic beam through the multi-layered media. We then show how the joint multi-frequency UMBIR reconstruction can be computed by modeling the direct arrival signals, detector noise, and incorporating a spatially varying image prior. Results using both simulated and experimental data indicate that multi-frequency UMBIR reconstruction yields much higher reconstruction quality than either single frequency UMBIR or SAFT.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

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

Effects of process intensification on homogeneity of an IgG1:κ monoclonal antibody during perfusion culture

The pharmaceutical industry employs various strategies to improve cell productivity. These strategies include process intensification, culture media improvement, clonal selection, media supplementation and genetic engineering of cells. However, improved cell productivity has inherent risk of impacting product quality attributes (PQA). PQAs may affect the products’ efficacy via stability, bioavailability, or in vivo bioactivity. Variations in manufacturing process may introduce heterogeneity in the products by altering the type and extent of N-glycosylation, which is a PQA of therapeutic proteins. We investigated the effect of different cell densities representing increasing process intensification in a perfusion cell culture on the production of an IgG1-κ monoclonal antibody from a CHO-K1 cell line. This antibody is glycosylated both on light chain and heavy chain. Our results showed that the contents of glycosylation of IgG1-κ mAb increased in G0F and fucosylated type glycans as a group, whereas sialylated type glycans decreased, for the mAb whole protein. Overall, significant differences were observed in amounts of G0F, G1F, G0, G2FS1, and G2FS2 type glycans across all process intensification levels. G2FS2 and G2 type N-glycans were predominantly quantifiable from light chain rather than heavy chain. It may be concluded that there is a potential impact to product quality attributes of therapeutic proteins during process intensification via perfusion cell culture that needs to be assessed. Since during perfusion cell culture the product is collected throughout the duration of the process, lot allocation needs careful attention to process parameters, as PQAs are affected by the critical process parameters (CPPs).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

2-D seismic wave propagation using the distributional finite-difference method: further developments and potential for global seismology

SUMMARY We present a time-domain distributional finite-difference scheme based on the Lebedev staggered grid for the numerical simulation of wave propagation in acoustic and elastic media. The central aspect of the proposed method is the representation of the stresses and displacements with different sets of B-splines functions organized according to the staggered grid. The distributional finite-difference approach allows domain-decomposition, heterogeneity of the medium, curvilinear mesh, anisotropy, non-conformal interfaces, discontinuous grid and fluid–solid interfaces. Numerical examples show that the proposed scheme is suitable to model wave propagation through the Earth, where sharp interfaces separate large, relatively homogeneous layers. A few domains or elements are sufficient to represent the Earth’s internal structure without relying on advanced meshing techniques. We compare seismograms obtained with the proposed scheme and the spectral element method, and we show that our approach offers superior accuracy, reduced memory usage, and comparable efficiency.

Geochemistry & Geophysics↗
DOE OSTI2020

The role of the spatial heterogeneity and correlation length of surface wettability on two-phase flow in a CO 2 -water-rock system

This study characterized and modeled heterogeneous surface wettability in sandstone and investigated the role of spatial heterogeneity and correlation length of surface wettability on relative permeability in a supercritical CO 2 (scCO 2 )-brine-rock system. Understanding the role of wettability heterogeneity on relative permeability is essential to geological CO 2 sequestration, oil and gas recovery, and contaminated groundwater remediation. Although numerous studies have attempted to understand the influences of surface wettability, capillary number (Ca), and viscosity ratio, the role of the spatial variation and correlation length of surface wettability on two-phase flow in three-dimensional (3D) porous media has not been unraveled due to the challenges in the measurement and representation of realistic rock surface wettability. In this work, we conducted in-situ measurements of surface contact angle (CA) in a Bentheimer sandstone after CO 2 flooding using micro-computed tomography (micro-CT), and found that the pore-scale CA distribution on rock surfaces followed a log-normal distribution associated with a spatial correlation length. Based on the statistical information from CT scanning, a Gaussian random field was used to model CA distributions that had desired standard deviations and spatial correlation lengths, which were then adjusted within a certain range of values for sensitivity analyses to study their combined effects on the two-phase flow in the porous medium using the lattice Boltzmann (LB) method. The LB two-phase flow simulation was accelerated using hybrid, multicore parallel computing to overcome the challenges in simulating multiphase flow in a large 3D domain having 800 × 800 × 600 nodes. The simulation results showed that the surface wettability heterogeneity (i.e., standard deviation of CA) had a lesser effect on the relative permeability of the wetting fluid (water) but a more significant impact on the relative permeability of the non-wetting fluid (scCO 2 ). The Corey model was used to fit the LB-simulated relative permeability curves of water and scCO 2 and showed that the variations in the relative permeability curves for both water and scCO 2 increased as the standard deviation and spatial correlation length of CA increased. This study illustrated that the assumption of homogeneous surface wettability may cause errors in multiphase flow simulations. Furthermore, the impacts of both the standard deviation and spatial correlation length of CAs should be accounted for. This is the first study that explored the spatial correlation lengths associated with CA distributions on sandstone surfaces and comprehensively investigated the roles of both spatial variation and correlation length of CA on two-phase flow properties in 3D porous media. The optimized LB multiphase flow model was proved a powerful tool to study the interplays and combined effects of these statistical parameters, which had critical applications in numerous natural and engineering processes that involved multiphase flow in porous media.

58 GEOSCIENCES↗
DOE OSTI2021

Advances and challenges in CO2 foam technologies for enhanced oil recovery in carbonate reservoirs

Utilization of surfactants for generating foam has the potential for fluid-mobility control during CO2 flooding processes, leading to improved oil recovery. This manuscript is a review of the most important aspects for the design of CO2-foams with mobility control in carbonate reservoirs, including recent advances in novel surfactants and analytical techniques for analysis of surfactant adsorption on carbonate minerals and their thermal stability. Several key parameters and properties regarding foam transport in porous media are reviewed such as, the minimum pressure gradient required for foam generation, the effects of the partition coefficient of a CO2-soluble surfactant on foam transport in rock cores, visualization of foam flow by either microfluidic or PET/CT imaging processes in heterogeneous or fractured systems and the effect of oil on foam transport under reservoir conditions. Also are discussed the most recent advances of two major foam modeling methods, including the semi-empirical STARS foam model and the population balance model, as well as CO2 foam pilot tests and the factors governing their successes.

Jian, Guoqing↗
DOE OSTI2021

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

Probability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude-based model for describing advective-dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub-continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a manyto-one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub-continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Additionally, differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude-based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characteristics of flow through randomly packed impermeable and permeable particles using pore resolved simulations

Pore resolved simulations are performed to study the mean flow characteristics in porous media formed by random distribution of impermeable and permeable particles. Here, permeability of the medium and reactive surface areas are analyzed for a range of solid fractions, packing patterns and flow rates. The results are compared with models presented in the literature and with simulations of flow through ordered packings. Permeability of a randomly packed porous medium is found to depend on the heterogeneous distribution of porosity which results from particle agglomeration and is typically ignored by simpler models. The reactive surface area is dependent on the flow field at pore scale in addition to geometry. Moreover, the wake formed downstream of the particles at moderately high flow rates also affects the reactive surface area. Similar analysis is also performed by assuming particles to be permeable to account for dual scale porosity relevant for various engineering applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗