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 163 records · Page 9

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↗

Longitudinal, Lateral, Vertical, and Temporal Thermal Heterogeneity in a Large Impounded River: Implications for Cold-Water Refuges

Dam operations can affect mixing of the water column, thereby influencing thermal heterogeneity spatially and temporally. This occurs by restricting or eliminating connectivity in longitudinal, lateral, vertical, and temporal dimensions. We examined thermal heterogeneity across space and time and identified potential cold-water refuges for salmonids in a large impounded river in inland northwestern USA. To describe these patterns, we used thermal infrared (TIR) imagery, in situ thermographs, and high-resolution, 3-D hydraulic mapping. We explained the median water temperature and probability of occurrence of cool-water areas using generalized additive models (GAMs) at reach and subcatchment scales, and we evaluated potential cold-water refuge occurrence in relation to these patterns. We demonstrated that (1) lateral contributions from tributaries dominated thermal heterogeneity, (2) thermal variability at confluences was approximately an order of magnitude greater than of the main stem, (3) potential cold-water refuges were mostly found at confluences, and (4) the probability of occurrence of cool areas and median water temperature were associated with channel geomorphology and distance from dam. These findings highlight the importance of using multiple approaches to describe thermal heterogeneity in large, impounded rivers and the need to incorporate these types of rivers in the understanding of thermal riverscapes because of their limited representation in the literature.

Mejia, Francine H.↗

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↗

Mapping of Heterogeneous Catalyst Degradation in Polymer Electrolyte Fuel Cells

Pt catalysts in polymer electrolyte fuel cells degrade heterogeneously as the catalyst particles are exposed to local variations throughout the catalyst layer during operation. State-of-the-art analytical techniques for studying degradation of Pt catalysts do not possess fine spatial resolution to elucidate such non-uniform degradation behavior at a large electrode level. A new methodology is developed to spatially resolve and quantify the heterogeneous Pt catalyst degradation over a large area (several cm ) of aged MEAs based on synchrotron X-ray microdiffraction. PEFC single cells are aged using voltage cycling as an accelerated stress test and the degradation heterogeneity at a micrometer length scale is visualized by mapping Pt catalyst particle size after voltage cycling. It is demonstrated in detail that the Pt catalyst particle size growth is non-uniform and follows the flow field geometry. The Pt particle size growth is greater in the area under the flow field land, while it is minimal in the area under the flow field channel. Additional non-uniformity is observed with the Pt particle size increasing more rapidly at the air outlet area than the Pt particle size at the inlet area.

36 MATERIALS SCIENCE↗

OpenABLext: An automatic code generation framework for agent-based simulations on CPU-GPU-FPGA heterogeneous platforms

The execution of agent-based simulations (ABSs) on hardware accelerator devices such as graphics processing units (GPUs) has been shown to offer great performance potentials. However, in heterogeneous hardware environments, it can become increasingly difficult to find viable partitions of the simulation and provide implementations for different hardware devices. To automate this process, we present OpenABLext, an extension to OpenABL, a model specification language for ABSs. By providing a device-aware OpenCL backend, OpenABLext enables the co-execution of ABS on heterogeneous hardware platforms consisting of central processing units, GPUs, and field programmable gate arrays (FPGAs).We present a novel online dispatching method that efficiently profiles partitions of the simulation during run-time to optimize the hardware assignment while using the profiling results to advance the simulation itself. In addition, OpenABLext features automated conflict resolution based on user-specified rules, supports graph-based simulation spaces, and utilizes an efficient neighbor search algorithm. We show the improved performance of OpenABLext and demonstrate the potential of FPGAs in the context of ABS. We illustrate how co-execution can be used to further lower execution times. OpenABLext can be seen as an enabler to tap the computing power of heterogeneous hardware platforms for ABS.

97 MATHEMATICS AND COMPUTING↗

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↗

Heterogeneous solute segregation suppresses strain localization in nanocrystalline Ag-Ni alloys

Solute segregation to individual grain boundaries is used by design to produce strong and stable nanocrystalline metallic alloys. Grain-boundary segregation, however, is known to cause adverse embrittlement effects from a strain-localization failure mechanism that imposes significant material limitations for structural applications. Here, using atomistic simulations, it is discovered that heterogeneous Ni segregation in nanocrystalline Ni-mixed Ag alloys dramatically shuts down localized shear bands during plastic deformation, while simultaneously increasing the tensile strength. Nanocrystalline Cu-mixed Ag metals are predicted to exhibit standard homogeneous Cu segregation and a tensile strength that saturates above a solute concentration of 8 at.% due to glass-like shear localization induced by grain boundaries. By contrast, it is found that heterogeneous Ni segregation in nanocrystalline Ag-Ni alloys forms solute-rich clusters along interfaces leading to strain delocalization at high strain and continuous strengthening at high solute concentrations up to 15 at.%. As a result, this study reveals the importance of heterogeneous versus homogeneous segregation behaviors on strain localization and points to a fundamentally new strategy to design failure-resistant nanostructured materials through grain boundary segregation engineering.

36 MATERIALS SCIENCE↗

Deformation-enhanced hierarchical multiscale structure heterogeneity in a Pd-Si bulk metallic glass

Here, the multiscale structures in a Pd 82 Si 18 binary bulk metallic glass before and after deformation were studied using electron microscopies, high-energy synchrotron X-ray diffraction, and small-angle scattering techniques. The experimental results revealed an enhancement of hierarchical structure heterogeneities on multiple length scales after deformation. Hierarchical multiple shear bands of high number density were observed after bending, introducing complex but periodically distributed residual strain. Pair distribution function analysis revealed that the connectivity of the short-range clusters on the medium-range scale determines the packing density difference between the tension side and the compression side in the sample after bending. In-situ synchrotron X-ray diffraction study also revealed a transformation of connection modes among short-range clusters under uniaxial tension and compression, which is consistent with those of triaxial tension/compression parts upon bending in Pd 82 Si 18 glassy alloys. The nanoscale heterogeneities for metallic glasses after deformation observed by small-angle scattering and transmission electron microscopy may be attributed to the nanoscale amorphous phase separation and interacting multiple shear bands enhanced by plastic deformation. Our findings suggested that the enhancement of hierarchical heterogeneous structure on multiple length scales may explain the excellent plasticity of Pd-Si glassy alloys, deepening the understanding of structure-property relation during plastic deformation in metallic glasses.

36 MATERIALS SCIENCE↗

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↗