Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heterogeneity”

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 127 records · Page 7

Heterogenous catalysis for oxygen tolerant photoredox atom transfer radical polymerization and small-molecule dehalogenation

Heterogeneous photocatalysts (PCs) have garnered attention for their sustainability and cost-effectiveness. Despite the existence of various types of these PCs, their synthesis often involves complex, multi-step procedures and laborious purification. Herein, we propose a simple method for attaching small-molecule photocatalytic species onto crosslinked 3-D polymer networks as insoluble scaffolds to create robust heterogeneous PCs. The highly swellable poly(ethylene glycol)-based ChemMatrix (CM) resin, known for its amphiphilic properties and high functional group loading, facilitated the covalent immobilization of the photoredox dye Eosin Y (EY), but also streamlined functionalization with Ir( III ) complexes. The resulting heterogeneous CM-EY demonstrated efficient photocatalytic performance in open-to-air dual photoredox catalysis of atom transfer radical polymerization (photo-ATRP) under green light. This was confirmed by the well-controlled synthesis of polymers with molecular masses ranging from 20 kDa to 300 kDa and low dispersities. Furthermore, CM-EY exhibited excellent photostability and recyclability over multiple cycles of ATRP. The heterogeneous catalysis of photo-ATRP provided high temporal control and enabled benign conditions for synthesizing protein-polymer hybrids (PPH). When combined with the initiator-modified CM (CM-BIB), CM-EY facilitated the solid-phase synthesis of homopolymers and block copolymers with recyclable performance. However, the coordinatively bound Ir@CM showed decreased catalytic activity and efficiency toward photoredox dehalogenation due to the leaching of active species during recycling. This study highlights the advantages of the covalent linking of catalysts to solid supports over non-covalent interactions, underscoring the potential of functionalized polymer resin as a promising scaffold. Such an approach offers customization and tunability, presenting opportunities for innovation in green chemistry.

Kapil, Kriti↗

IRIS-DMEM: Efficient Memory Management for Heterogeneous Computing

This paper proposes an efficient data memory management approach for the Intelligent RuntIme System (IRIS) heterogeneous computing framework along with new data transfer policies. IRIS provides a task-based programming model for extreme heterogeneous computing (e.g., CPU, GPU, DSP, FPGA) with support for today's most important programming languages (e.g., OpenMP, OpenCL, CUDA, HIP, OpenACC). However, the IRIS framework either forces the programmer to introduce data transfer commands for each task or relies on suboptimal memory management for automatic and transparent data transfers. The work described here extends IRIS with novel heterogeneous memory handling and introduces novel data transfer policies by employing the Distributed data MEMory handler (DMEM) for efficient and optimal movement of data among the various computing resources. The proposed approach achieves performance gains of up to 7× for tiled LU factorization and tiled DGEMM (i.e., matrix multiplication) benchmarks. Moreover, this approach also reduces data transfers by up to 71% when compared to previous IRIS heterogeneous memory management handlers. This work compares the performance results of the IRIS framework's novel DMEM with the StarPU runtime and MAGMA math library for GPUs. Experiments show a performance gain of up to 1.95× over StarPU and 2.1× over MAGMA.

Miniskar, Narasinga Rao↗

Pursuing sources of heterogeneity in modeling clustered population

Abstract Researchers often have to deal with heterogeneous population with mixed regression relationships, increasingly so in the era of data explosion. In such problems, when there are many candidate predictors, it is not only of interest to identify the predictors that are associated with the outcome, but also to distinguish the true sources of heterogeneity , that is, to identify the predictors that have different effects among the clusters and thus are the true contributors to the formation of the clusters. We clarify the concepts of the source of heterogeneity that account for potential scale differences of the clusters and propose a regularized finite mixture effects regression to achieve heterogeneity pursuit and feature selection simultaneously. We develop an efficient algorithm and show that our approach can achieve both estimation and selection consistency. Simulation studies further demonstrate the effectiveness of our method under various practical scenarios. Three applications are presented, namely, an imaging genetics study for linking genetic factors and brain neuroimaging traits in Alzheimer's disease, a public health study for exploring the association between suicide risk among adolescents and their school district characteristics, and a sport analytics study for understanding how the salary levels of baseball players are associated with their performance and contractual status.

Li, Yan↗

Interplay of Electrode Heterogeneity and Lithium Plating

Lateral microstructure heterogeneity in anodes is known to induce nonuniform current density, state of charge, and lithium plating. This means that such electrode heterogeneity can limit the fast charging of lithium-ion batteries. In this work, a combination of experiments and simulation is employed to understand the effect of mm scale lateral heterogeneity on cell aging. A previously developed model was extended to efficiently simulate SEI formation and Li plating for independent regions of an electrode. The model consists of three parallel regions each described under a P2D framework and with a distinct ionic resistance and possibly active material loading. The results suggest that during fast charge when the active material is uniformly distributed across the three regions, the region with the highest resistance reaches the end of life sooner than the other regions. There is also positive feedback from Li metal filling the pores near the separator interface that further accelerates lithium plating. Lastly, when there is a non-uniform active material distribution associated with the ionic resistance heterogeneity, tight competition between regions can occur, leading to less overall lithium plating and plating that is more uniform between regions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microstructure Scale Lithium-Ion Battery Modeling: Part II. On In-Plane Heterogeneities and the Mechanisms that Regulate Them

Li-ion batteries performance and degradation are typically modeled at the macroscopic scale, that is neglecting in-plane heterogeneities that can arise from non-uniform electrode microstructures. Herein, a microstructure scale electrochemical model is used to quantify the impact of microstructure heterogeneity on cell performance during fast charging. The model predicts the electrolyte and solid concentration in-plane standard deviation can reach, respectively, ≈200 mol·m −3 and 6–7 kmol·m −3 locally. Further, the intercalation current density in-plane relative standard deviation can reach extremely high values, around 100% in the cathode and well above 100% in the anode graphite. These denote highly non-uniform lithiation rates and material utilization within each slice of the microstructure along the cell thickness. Non-uniform curvatures, at the particle scale (surface roughness) and between particles (size distribution), were found to initiate these in-plane heterogeneities, while an OCP-induced mechanism subsequently regulates them. The present model provides new insights into small length scale heterogeneity impact on battery performance not available with standard macro-scale/P2D modeling.

25 ENERGY STORAGE↗

Multi-scale Analysis of Surface Heterogeneity Induced Convection on Isentropic Coordinates

This research uses semi-idealized simulations to investigate multi-scale processes induced by the heterogeneity of soil moisture observed during the 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign. The semi-idealized simulations have realistic land heterogeneity, but large-scale winds are removed. Analysis on isentropic coordinates enables the tracking of circulation that transports energy vertically and facilitates the identification of the primary convective processes induced by realistic land heterogeneity. The isentropes associated with upward motion are found to connect the ground characterized by high latent heat flux to cloud bases directly over the ground with high sensible heat flux, while isentropes associated with downward motion connect precipitation to the ground characterized by high sensible heat fluxes. The mixing of dry, warm parcels ascending from the ground with high sensible heat fluxes and moist parcels from high latent heat regions leads to cloud formation. This new mechanism explains how soil moisture heterogeneity provides the key ingredients such as buoyancy and moisture for shallow cloud formation. We also found that the sub-mesoscale dominates upward energy transport in the boundary layer, while mesoscale circulations contribute to vertical energy transport above the boundary layer. Our novel method better illustrates and elucidates the nature of land atmospheric interactions under irregular and realistic soil moisture patterns.

54 ENVIRONMENTAL SCIENCES↗

Coreflood Tests to Evaluate Enhanced Oil Recovery Potential of Wettability-Altering Surfactants for Oil-Wet Heterogeneous Carbonate Reservoirs

Summary Oil-wetness and heterogeneity are two main factors that result in low oil recovery (OR) by waterflood in carbonate reservoirs. The injected water is likely to flow through high-permeability regions and bypass the oil in the low-permeability matrix. In this study, systematic coreflood tests were carried out in both “homogeneous” cores and “heterogeneous” cores with a wettability-altering surfactant. The homogeneous coreflood tests were conducted to evaluate surfactant retention, as well as to compare tertiary surfactant flooding with secondary surfactant flooding. The heterogeneous coreflood test was proposed to model bypassing in low-permeability matrix during waterfloods, and dynamic imbibition of surfactant into the low-permeability matrix. Surfactant retention results suggest that retention increases as initial oil saturation decreases. The retention of selected surfactant in the target reservoir cores was measured to be within a range of 0.07–0.12 mg/g-rock, which is economically acceptable. The results of homogeneous coreflood tests showed that both secondary waterflood and secondary surfactant flood can achieve high OR (>50%) from relatively homogeneous oil-wet cores. A shut-in phase after the surfactant injection resulted in a surge in oil production, which suggests that enough time should be given for wettability alteration by surfactants. The results of heterogeneous coreflood tests showed that more oil is bypassed in the tighter matrix by waterflood if the permeability is higher in the flooded layer and this bypassed oil is the target for the wettability-altering surfactant floods. Slow wettability-altering surfactant injection leads to imbibition into bypassed regions. When the oil-wet carbonate reservoirs have large unswept regions after waterflood, wettability-altering surfactants can significantly improve OR if enough time is given for imbibition.

Engineering↗

Single molecule-based fliFISH validates radial and heterogeneous gene expression patterns in pancreatic islet ß cells

Single cell RNA sequencing (scRNA-Seq) technologies have greatly enhanced our understanding of islet cell transcriptomes and have revealed the existence of ß cell heterogeneity. However, comparison of scRNA-Seq datasets from different groups have highlighted inconsistencies in gene expression patterns, primarily due to variable detection of lower abundance transcripts. Furthermore, such analyses are unable to uncover the spatial organization of heterogeneous gene expression. Here we used fluctuation localization imaging-based fluorescence in situ hybridization (fliFISH) to quantify transcripts in single cells in mouse pancreatic islet sections. We compared the expression patterns of Insulin 2 (Ins2) with Mafa and Ucn3 – two genes expressed in ß cells as they mature, as well as Rgs4 – a factor with variably reported expression in the islet. This approach accurately quantified transcripts across a wide range of expression levels - from single copies to over hundred copies per cell in one islet. Importantly, fliFISH allowed evaluation of transcript heterogeneity in the spatial context of an intact islet. These studies confirm the existence of a high degree of heterogeneous gene expression levels within the islet and highlight relative and radial expression patterns that likely reflect distinct ß cell maturation states along the radial axis of the islet.

Li, Fangjia↗

Design and analysis of a wake model for spatially heterogeneous flow

Abstract. Methods of turbine wake modeling are being developed to more accurately account for spatially variant atmospheric conditions within wind farms. Most current wake modeling utilities are designed to apply a uniform flow field to the entire domain of a wind farm. When this method is used, the accuracy of power prediction and wind farm controls can be compromised depending on the flow-field characteristics of a particular area. In an effort to improve strategies of wind farm wake modeling and power prediction, FLOw Redirection and Induction in Steady State (FLORIS) was developed to implement sophisticated methods of atmospheric characterization and power output calculation. In this paper, we describe an adapted FLORIS model that features spatial heterogeneity in flow-field characterization. This model approximates an observed flow field by interpolating from a set of atmospheric measurements that represent local weather conditions. The objective of this method is to capture heterogeneous atmospheric effects caused by site-specific terrain features, without explicitly modeling the geometry of the wind farm terrain. The implemented adaptations were validated by comparing the simulated power predictions generated from FLORIS to the actual recorded wind farm output from the supervisory control and data acquisition (SCADA) recordings and large eddy simulations (LESs). When comparing the performance of the proposed heterogeneous model to homogeneous FLORIS simulations, the results show a 14.6 % decrease for mean absolute error (MAE) in wind farm power output predictions for cases using wind farm SCADA data and a 18.9 % decrease in LES case studies. The results of these studies also indicate that the efficacy of the proposed modeling techniques may vary with differing site-specific operational conditions. This work quantifies the accuracy of wind plant power predictions under heterogeneous flow conditions and establishes best practices for atmospheric surveying for wake modeling.

17 WIND ENERGY↗

MEPHESTO: Modeling Energy-Performance in Heterogeneous SoCs and Their Trade-Offs

Integrated shared memory heterogeneous architectures are pervasive because they satisfy the diverse needs of mobile, autonomous, and edge computing platforms. Although specialized processing units (PUs) that share a unified system memory improve performance and energy efficiency by reducing data movement, they also increase contention for this memory since the PUs interact with each other. Prior work has investigated performance degradation due to memory contention, but few have studied the relationship of power and energy to memory contention. Moreover, a comprehensive solution that models memory contention for kernel placement on contemporary heterogeneous systems on chip (SoCs) in response to energy and performance has been largely unaddressed.This paper presents MEPHESTO, a novel and holistic approach for managing this balance. The authors characterize applications and PUs in terms of two memory contention factors - time factors and power factors - to achieve the desired trade-off between energy and performance for collocated kernel execution on heterogeneous systems. The authors believe that this investigation is the first to combine all of these factors and present a simple knob-based approach that expresses the target trade-off. The approach is evaluated on a diverse integrated shared memory heterogeneous system with a CPU, GPU, and programmable vision accelerator. By using an empirical model for memory contention that provides up to 92% accuracy, the kernel collocation approach can provide a near-optimal ordering and placement based on the user-defined, energy-performance trade-off parameter. Moreover, the dynamic programming-based heuristics provide up to 30% better energy or 20% performance benefits when compared with the greedy approaches commonly employed by previous studies.

Alaul haque monil, Mohammad↗

Impact of sub-core scale heterogeneity on CO2/brine multiphase flow for geological carbon storage in the Minnelusa sandstone

CO2 geological storage in deep saline aquifers is a mitigation option for CO2 emissions due to its large storage capacity and immediate accessibility. Accurately determining the CO2-brine relative permeability curves is key to the evaluation of CO2 injectivity and sweep efficiency in reservoir simulation as well as the CO2 injection in the field. This study highlights the remarkable effects of sub-core scale heterogeneity on the CO2-brine multiphase flow properties of the Minnelusa Sandstone in Wyoming. Two unsteady state CO2-brine drainage experiments were performed on the two samples. The first sample exhibits slanted laminated structure, while the second one represents a more homogeneous sandstone system. The CO2 saturation distributions during drainage reveals that the main variation in multiphase flow properties of two core samples were attributed to the porosity distribution that leads to the capillary pressure heterogeneity. Assisted history matching was used to obtain the respective relative permeability curves, which suggests heterogeneity-dependent behavior. In addition, sensitive and uncertainty analyses indicate that physical and petro-physical properties of low-porosity and low-permeability bedding layers exert marked effects on CO2 front breakthrough time and brine production. The results presented in this study help to gain insight into the CO2-brine multiphase flow properties in heterogeneous sandstones and can pave the way for the upscaling of CO2 migration and field-scale simulation accurately. This work is funded under the Department of Energy CarbonSAFE program (awards DE-FE0031624 DE0031891).

Kou, Zuhao↗

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗

Workshop on Addressing Rigor and Reproducibility in Thermal, Heterogeneous Catalysis

Heterogeneous catalysis has long served as the bedrock of the manufacturing of energy carriers, fuels and chemicals, and various technologies for pollution abatement. The significant complexity and variability spanning the entire breadth of catalyst material properties, synthesis methods, characterization techniques, and evaluation procedures, has focused attention on the need to establish community-accepted best practices for ensuring high-quality, benchmarked, and reproducible data. In addition, increased societal urgency to transition to clean energy and reduce greenhouse gas concentrations has incentivized interdisciplinary, convergent, and translational approaches to catalysis research in recent years. Research engineers and scientists with expertise cutting broadly across materials science, chemical synthesis, interfacial science, spectroscopy, and methods of data science and computational simulation, all bring diverse and important perspectives to catalysis research, but often with little awareness of the complexity of catalytic systems, especially in their working environment. As has already occurred in other scientific fields, there has been growing recognition and consensus in the heterogeneous catalysis research community that mechanisms are needed to improve the rigor and reproducibility (R&R) of experimental measurements, to ensure alignment of the broader research community with a common core of best practices specific to the realization of high-quality catalysis research. Similarly, the field is moving rapidly toward computationally informed and data science-driven catalyst design, but the success of implementing such predictive tools hinges on model training and validation rooted in rigorously obtained and reproducible experimental data that are benchmarked to common specifications. As such, this workshop was convened to prepare a report summarizing best practices for reporting data and performing experiments that researchers can use to benchmark, validate, and reproduce data in specific sub-fields of thermal, heterogeneous catalysis. Additionally, we discussed recommendations for future actions that may improve R&R in this field. The workshop organizers and participants include a diverse range of catalysis researchers from various employment sectors (e.g., academia, industry, national laboratory), institutional mission and resources (e.g., PhD-granting research universities, non-PhD-granting teaching universities), career stage (e.g., early, mid and late-career), technical expertise, and demographic background. This diverse group was involved in the discussion of workshop agenda items, writing this report, and discussing possible future action items for the community to consider, which helped ensure that a broad range of perspectives were captured in the description of the problems at hand and the creation of actionable solutions that may be effectively adopted by the diverse practitioners in catalysis research. Importantly, this group of workshop participants also included very early career researchers (e.g., senior PhD students, postdoctoral scholars) who will become the next generation of scientific leaders in various sectors, thus capturing emerging perspectives of newcomers to the field to shape its future while positively impacting the development of its future workforce. We envision that this effort will help advance the field of catalysis science by improving the rigor and reproducibility of experimental data collected by current researchers and future newcomers to the field, which is of broad importance to health and vitality of any scientific discipline. Therefore, best practices identified in this endeavor for thermal heterogeneous catalysis can be translated to such efforts in other areas of catalysis and other scientific fields involving the study of materials, and vice versa. We also envision this to be an ongoing effort, with future workshops that are convened to discuss issues of rigor and reproducibility on technical topics that were unable to be covered in this workshop due to its scope limitations, and as emerging methods and materials become more prevalent in the research community.

36 MATERIALS SCIENCE↗

Results from single-shock Marble experiments studying thermonuclear burn in the presence of heterogeneous mix on the National Ignition Facility

The Marble campaign on the National Ignition Facility investigates the effect of heterogeneous mix on thermonuclear burn for comparison to a probability distribution function (PDF) burn model. Marble utilizes plastic capsules filled with deuterated plastic foam and a fill gas containing tritium. As the capsules implode, the deuterium in the foam mixes with the tritium gas, and DT neutrons are produced as the shocks compress and heat the mixture. The yield of DT neutrons is dependent on the uniformity of the mix, with more heterogeneous mix producing fewer neutrons. In Marble, the heterogeneity of the mix is controlled by varying the diameter of voids introduced into the foam. The first NIF Marble campaign has been executed in which the Marble capsules were indirectly driven with a single strong shock using NIF hohlraums. The experiments produce a low-convergence, high-ion-temperature implosion. The ratio of DT to DD neutron yield is largely consistent with uniform atomic mix for fine-pore foam, and increases slightly with void diameter, contrary to 1D simulations using the PDF burn model. Recent 3D high-resolution simulations of similar experiments performed on the Omega Laser Facility suggest an explanation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Tailoring the Reaction Heterogeneity for Robust Li‐Rich Cathodes

The practical application of Li-rich Mn-based layered oxides (LLO) cathode is hindered by severe capacity and voltage degradation resulting from severe oxygen release and irreversible phase transition. In this work, the reaction heterogeneity, which describes the spatially resolved electrochemical divergence within individual cathode particles, is engineered through compositional gradient design to couple Li + transport kinetics and the anion redox activity between particle interiors and surfaces. It is revealed that Co/Mn concentration gradient within particles creates heterogeneous phase content distribution and structural ordering, inducing surface-bulk reaction heterogeneity that significantly impacts the overall electrochemical performance. Specifically, Li 2 MnO 3 -poor and Co-enriched surface effectively mitigates the oxygen loss and enhances electrochemical reaction kinetics, benefited from the reduced surface redox reactivity and induced highly ordered intra-layered cationic arrangement. Meanwhile, the Li 2 MnO 3 -enriched core with slight Li/Ni intermixing provides high reversible capacity and strong mechanical stability. Consequently, the greatly enhanced anion redox reversibility, Li + diffusion dynamics, and structure stability endow LLO with exceptional electrochemical properties, showing a capacity retention of 86.0% and a reduced voltage decay of 0.518 mV per cycle after 500 cycles at 1 C. This work provides a valuable strategy to tailor the redox chemistry and achieve robust LLO.

25 ENERGY STORAGE↗

Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media

Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.

97 MATHEMATICS AND COMPUTING↗

The Method of Finite Averages: A rigorous upscaling methodology for heterogeneous porous media

Rigorous upscaling techniques offer accurate and computationally-efficient strategies for modeling the average behaviors of multi-physical, multiscale phenomena in geological porous media. However, such techniques often rely on a variety of methodological assumptions that prohibit their rigorous application to practical systems (e.g., systems involving heterogeneous porous media, system-scale boundary conditions, and fine-scale dynamics that are not diffusion-dominant). In this work, we aim to formulate an upscaling methodology with few methodological assumptions to provide high levels of model generality and foster the utilization of rigorously-derived upscaled models in practice. In particular, we introduce the Method of Finite Averages (MoFA), a novel upscaling methodology for rigorously modeling heterogeneous porous media and system-scale boundary conditions. We then detail MoFA’s implementation for the advective–diffusive transport of a single species and compare the methodology with classic numerical techniques, as well as other rigorous upscaling techniques, to highlight MoFA’s unique combination of rigor and generality. We then validate the derived model while demonstrating its benefits in three numerical experiments. The results suggest that (1.) the applicability and a priori error guarantees of MoFA models do not directly depend on system geometry, (2.) a model’s applicability and error guarantees can be can arbitrarily expanded and reduced, respectively, with further computational expense, and (3.) downscaling with MoFA provides an efficient strategy for generating accurate pore-scale solutions from upscaled results. Ultimately, the results evidence that upscaled models can be rigorously derived for heterogeneous porous media systems and resolved in a fraction of the time it takes to perform the equivalent pore-scale simulations.

58 GEOSCIENCES↗

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗