[Paper Review] Homogenization of a locally-periodic medium with areas of low and high diffusivity

T.L. van Noorden, Adrian Muntean|arXiv (Cornell University)|Mar 21, 2010

Advanced Mathematical Modeling in Engineering3 citations

TL;DR

This paper develops a formal homogenization framework for locally-periodic heterogeneous media with spatially varying circular inclusions of low diffusivity, deriving a two-scale reaction-diffusion model with x-dependent effective coefficients. The key contribution is a rigorous weak solvability proof for the resulting model using x-dependent Bochner spaces and a Schauder fixed-point argument, enabling analysis of micro-macro transport in porous materials with structured diffusivity.

ABSTRACT

We aim at understanding transport in porous materials including regions with both high and low diffusivities. For such scenarios, the transport becomes structured (here: {\em micro-macro}). The geometry we have in mind includes regions of low diffusivity arranged in a locally-periodic fashion. We choose a prototypical advection-diffusion system (of minimal size), discuss its formal homogenization (the heterogenous medium being now assumed to be made of zones with circular areas of low diffusivity of $x$-varying sizes), and prove the weak solvability of the limit two-scale reaction-diffusion model. A special feature of our analysis is that most of the basic estimates (positivity, $L^\infty$-bounds, uniqueness, energy inequality) are obtained in $x$-dependent Bochner spaces.

Motivation & Objective

• To model transport in porous materials with spatially varying diffusivity, particularly when regions of low diffusivity are arranged in a locally-periodic fashion.
• To extend classical homogenization by accounting for microstructure that varies with macroscopic position x, moving beyond purely periodic assumptions.
• To derive a two-scale reaction-diffusion model that preserves microstructural geometry while capturing effective macroscopic behavior.
• To establish weak solvability of the two-scale model using x-dependent Bochner spaces, ensuring positivity, L∞-bounds, and energy inequalities.
• To lay the groundwork for future rigorous convergence rate analysis as ε → 0, where ε is the microscale parameter.

Proposed method

• Formal asymptotic expansion using a level-set-based description of the microstructure boundary, resembling the boundary unfolding operator.
• Modeling the medium as a collection of circular low-diffusivity inclusions whose size varies with macroscopic position x.
• Deriving a two-scale model where the solution depends on both macroscopic variable x and microscopic variable y, with x-dependent effective coefficients (porosity, tortuosity).
• Employing x-dependent Bochner spaces to formulate the solution space, enabling estimation of key properties like positivity and L∞-bounds.
• Proving weak solvability via a Schauder fixed-point argument, offering an alternative to semigroup or Banach fixed-point approaches in transformed domains.
• Using energy estimates, trace inequalities, and Gronwall’s inequality to control the growth of solutions and ensure boundedness in the limit.

Experimental results

Research questions

• RQ1 How can homogenization be formalized in a deterministic framework when the microstructure is locally periodic but not globally periodic?
• RQ2 What is the structure of the effective macroscopic model when the size of low-diffusivity inclusions varies with position x?
• RQ3 How can weak solvability be established for a two-scale reaction-diffusion system with x-dependent coefficients and microstructure?
• RQ4 What mathematical tools are necessary to maintain estimates (positivity, L∞-bounds, energy inequality) in x-dependent function spaces?
• RQ5 Can the formal homogenization result be made rigorous, and what techniques are needed to prove convergence rates as ε → 0?

Key findings

• A two-scale reaction-diffusion model is formally derived, with effective transport coefficients (porosity, permeability, tortuosity) that depend explicitly on the macroscopic position x.
• The model preserves the detailed geometry of the microstructure, including x-dependent circular inclusions of low diffusivity, avoiding oversimplification via classical homogenization.
• Weak solvability of the two-scale model is proven using a Schauder fixed-point argument in x-dependent Bochner spaces, ensuring existence of global-in-time solutions.
• Key estimates—positivity, L∞-bounds, uniqueness, and energy inequality—are rigorously established in x-dependent Bochner spaces, a novel feature of the analysis.
• The analysis provides a framework for handling free-boundary problems and time-dependent microstructures, as demonstrated by the applicability to dissolution-precipitation and fast-reaction slow-diffusion scenarios.
• The results lay the foundation for future rigorous convergence rate analysis as ε → 0, with the authors suggesting a fusion of corrector estimates and two-scale structure as a path forward.