pith. sign in

arxiv: 2604.25722 · v1 · submitted 2026-04-28 · 🧮 math.NA · cs.NA

Computational homogenization of unsteady flows in a periodic porous medium

P.N. Vabishchevich This is my paper

Pith reviewed 2026-05-07 15:01 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords computational homogenizationporous mediaunsteady flowsmemory kernelDarcy lawfinite element methodnumerical approximation
0
0 comments X

The pith

Approximating the memory kernel as a sum of exponentials converts the nonlocal macroscale model of unsteady porous flows into a local differential system.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper presents a computational approach to homogenize unsteady flows of viscous incompressible fluids through periodic porous media. The macroscale description involves an integro-differential Darcy law whose memory kernel comes from solving unsteady cell problems on the periodic unit cell. By decomposing the conductivity tensor into steady and unsteady parts and approximating the kernel with exponentials, the nonlocal problem is reduced to a local system of differential equations. This allows the use of standard finite-element methods in space and stable two-level time integration schemes. The method is demonstrated on a two-dimensional test problem for filtration in porous media, showing how memory effects can be handled without direct microscale resolution at large scales.

Core claim

The work establishes that the unsteady homogenization of viscous flows in periodic porous media leads to a macroscopic integro-differential Darcy law with a tensor memory kernel computed from cell problems. This kernel is split into steady-state and unsteady components obtained from auxiliary boundary-value and spectral problems on the periodicity cell. Approximating the memory kernel as a sum of exponentials then transforms the nonlocal macro problem into a local system of differential equations, for which spatial finite-element approximations and stable two-level time schemes are constructed and tested on a two-dimensional example.

What carries the argument

The sum-of-exponentials approximation to the memory kernel determined from periodicity cell problems, which enables reduction of the integro-differential macroscale equations to a local differential system.

Load-bearing premise

The memory kernel from the cell problems admits a sufficiently accurate approximation by a finite sum of exponentials over the time scales of interest, without introducing significant errors in the macroscopic solution.

What would settle it

Numerical experiments on a periodic porous structure whose memory kernel decays slowly would reveal whether the exponential sum can capture the long-time behavior without large deviations from the full integro-differential solution.

Figures

Figures reproduced from arXiv: 2604.25722 by P.N. Vabishchevich.

Figure 1
Figure 1. Figure 1: Periodic porous medium and periodicity cell view at source ↗
Figure 2
Figure 2. Figure 2: Auxiliary functions w0 1 · e1 (left) and w0 2 · e1 (right) view at source ↗
Figure 3
Figure 3. Figure 3: Pressure for different pore geometries. The approximation of the memory kernel is performed by solving the spectral problem (3.1)–(3.3) in the domain Yf . The accuracy of the eigenvalue computation is controlled by calculations on a sequence of refined meshes. We use three computational meshes: Mesh 1 — characteristic mesh size h = 0.02, 2,365 nodes; Mesh 2 — h = 0.01, 8,973 nodes; Mesh 3 — h = 0.005, 35,0… view at source ↗
Figure 4
Figure 4. Figure 4: Eigenfunctions: top row — φ1 · e1 (left) and φ2 · e1 (right), bottom row — φ3 · e1 (left) and φ4 · e1 (right). 15 view at source ↗
Figure 5
Figure 5. Figure 5: Eigenvalues λk, k = 1, 2, . . . , m. 0 20 40 60 80 100 k 10 9 10 8 10 7 10 6 10 5 10 4 10 3 10 2 10 1 10 0 |a k 1 | 0 20 40 60 80 100 k 10 9 10 8 10 7 10 6 10 5 10 4 10 3 10 2 10 1 10 0 |a k 2 | view at source ↗
Figure 6
Figure 6. Figure 6: Coefficients |a k 1 | (top) and |a k 2 |, k = 1, 2, . . . , m (bottom). 16 view at source ↗
Figure 7
Figure 7. Figure 7: Elements of the steady-state conductivity tensor: diagonal (top) and off-diagonal (bottom). view at source ↗
Figure 8
Figure 8. Figure 8: Elements of the conductivity tensor Ke(ϵ) after filtering: diagonal (top) and off-diagonal (bottom). 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 x1 0.0 0.2 0.4 0.6 0.8 1.0 x2 t = 0.0 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 x1 0.0 0.2 0.4 0.6 0.8 1.0 x2 t = 0.00025 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00 x1 0.… view at source ↗
Figure 9
Figure 9. Figure 9: Solution of the macroscale problem taking into account memory effects at different time instants. view at source ↗
read the original abstract

The work is devoted to the development and computational implementation of the homogenization method for modeling unsteady flows of a viscous incompressible fluid in periodic porous media taking into account memory effects. At the macrolevel, the flow is described by an integro-differential Darcy law with a tensor memory kernel determined by solving unsteady problems on the periodicity cell. The developed approach to computational homogenization is based on finding the steady-state and unsteady components of the conductivity tensor from solving auxiliary boundary value and spectral problems on the periodicity cell. The nonlocal macroscopic problem is transformed into a local system of differential equations by approximating the memory kernel as a sum of exponentials. Issues of spatial finite element approximation are discussed, and stable two-level schemes in time are constructed. The results of applying the developed computational homogenization technology for unsteady filtration problems in porous media to a two-dimensional test problem are presented.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper develops a computational homogenization procedure for unsteady flows of a viscous incompressible fluid through periodic porous media. It derives an integro-differential macroscale model featuring a memory kernel computed via auxiliary steady and spectral problems on the periodicity cell. The key step is approximating this kernel by a finite sum of exponentials to obtain an equivalent local system of differential equations. Finite element discretization in space and stable two-level time schemes are presented, along with results from a two-dimensional numerical test.

Significance. Should the exponential approximation be shown to introduce controllable errors and the overall method validated with convergence studies, this approach would offer an efficient alternative to direct simulation of nonlocal macroscale problems in porous media flow. The decomposition into steady and unsteady conductivity components via cell problems is a solid technical foundation. The emphasis on stable time discretization addresses an important practical aspect of the unsteady problem. This could be significant for applications in filtration and groundwater modeling where memory effects are relevant.

major comments (2)
  1. [Section 3] The localization of the integro-differential macroscale problem via exponential approximation of the memory kernel (described after the cell problem solutions) lacks any a priori error estimate or numerical study of the approximation error's impact on the macroscale velocity or pressure fields. This is central to validating the method's accuracy.
  2. [Section 5 (Numerical results)] The two-dimensional test problem is presented without quantitative error measures, reference solutions, or convergence rates with respect to spatial mesh size, time step, or the number of exponential terms used in the kernel approximation. This makes it difficult to assess the practical performance and reliability of the proposed schemes.
minor comments (2)
  1. [Notation] The distinction between the full conductivity tensor and its steady/unsteady parts could be made clearer, perhaps with a summary table of definitions.
  2. [Time schemes] The stability proof for the two-level schemes would benefit from explicit statement of the conditions on the time step size.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive review of our manuscript. The comments identify important aspects for strengthening the validation of the proposed computational homogenization approach. We address each major comment below and indicate the revisions we will incorporate.

read point-by-point responses

  1. Referee: [Section 3] The localization of the integro-differential macroscale problem via exponential approximation of the memory kernel (described after the cell problem solutions) lacks any a priori error estimate or numerical study of the approximation error's impact on the macroscale velocity or pressure fields. This is central to validating the method's accuracy.

    Authors: We agree that a dedicated analysis of the exponential approximation error is essential. The approximation arises from truncating the spectral expansion of the cell-problem operator, which guarantees that the kernel error decays exponentially with the number of retained terms. While the current manuscript does not contain an a priori bound or a systematic numerical study of its propagation to the macroscale fields, we will add both: (i) a brief discussion of the spectral truncation error together with a reference to existing approximation theory for exponential sums, and (ii) a new numerical subsection that varies the number of exponential terms and reports the resulting changes in macroscale velocity and pressure for the test geometry. These additions will be placed in Section 3 and cross-referenced in the numerical results. revision: yes

  2. Referee: [Section 5 (Numerical results)] The two-dimensional test problem is presented without quantitative error measures, reference solutions, or convergence rates with respect to spatial mesh size, time step, or the number of exponential terms used in the kernel approximation. This makes it difficult to assess the practical performance and reliability of the proposed schemes.

    Authors: The numerical section was intended to illustrate the overall workflow rather than to serve as a comprehensive convergence study. We acknowledge that the absence of quantitative error tables and convergence rates limits the assessment of the method. In the revised manuscript we will augment Section 5 with: (i) a reference solution obtained on a sequence of successively refined meshes using a direct (non-homogenized) fine-scale simulation, (ii) L2-norm errors for velocity and pressure as functions of spatial mesh size h, time step τ, and the number N of exponential terms, and (iii) observed convergence rates extracted from these data. The two-level time schemes will also be tested for stability on the refined meshes. These quantitative results will be presented in tables and figures. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper derives the integro-differential macroscale model directly from independent auxiliary boundary-value and spectral problems solved on the periodicity cell; these cell problems furnish the conductivity tensor and memory kernel without reference to the final macroscale solution. The subsequent replacement of the kernel by a finite sum of exponentials is presented explicitly as a numerical approximation to obtain a local differential system, not as an identity or a fitted parameter whose value is taken from the target macroscale output. No load-bearing self-citations, uniqueness theorems imported from prior author work, or re-labeling of known results appear in the derivation. The overall procedure remains self-contained: cell-problem data are computed first, the approximation is applied afterward, and the resulting scheme is discretized and tested on an independent two-dimensional example.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The approach rests on standard domain assumptions of periodicity and fluid incompressibility plus the practical choice of how many exponential terms to retain; no new physical entities are introduced.

free parameters (1)
  • Number of exponential terms
    Chosen to approximate the memory kernel; balances accuracy against the size of the resulting ODE system.
axioms (2)
  • domain assumption The porous medium possesses a periodic microstructure
    Invoked to justify reduction of the microscale Stokes problem to a single periodicity cell.
  • domain assumption The fluid is viscous and incompressible
    Standard assumption underlying the microscale flow equations whose solutions define the memory kernel.

pith-pipeline@v0.9.0 · 5435 in / 1295 out tokens · 59038 ms · 2026-05-07T15:01:57.921769+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Bear, Dynamics of Fluids in Porous Media, Courier Corporation, 2013

    J. Bear, Dynamics of Fluids in Porous Media, Courier Corporation, 2013

    work page 2013

  2. [2]

    M. P. Anderson, W. W. Woessner, R. J. Hunt, Applied Groundwater Modeling: Simulation of Flow and Advective Transport, Academic Press, 2015

    work page 2015

  3. [3]

    J. Bear, A. H.-D. Cheng, Modeling Groundwater Flow and Contaminant Transport, Springer, 2010

    work page 2010

  4. [4]

    M. B. Allen, The Mathematics of Fluid Flow Through Porous Media, John Wiley & Sons, 2021

    work page 2021

  5. [5]

    Sánchez-Palencia, Non-Homogeneous Media and Vibration Theory, Springer, 1980

    E. Sánchez-Palencia, Non-Homogeneous Media and Vibration Theory, Springer, 1980

    work page 1980

  6. [6]

    Auriault, C

    J.-L. Auriault, C. Boutin, C. Geindreau, Homogenization of Coupled Phenomena in Heterogenous Media, John Wiley & Sons, 2010

    work page 2010

  7. [7]

    Hornung, Homogenization and Porous Media, Springer Science & Business Media, 1997

    U. Hornung, Homogenization and Porous Media, Springer Science & Business Media, 1997

    work page 1997

  8. [8]

    Mikelić, Homogenization theory and applications to filtration through porous media, in: Filtration in Porous Media and Industrial Application, Springer, 2007, pp

    A. Mikelić, Homogenization theory and applications to filtration through porous media, in: Filtration in Porous Media and Industrial Application, Springer, 2007, pp. 127–214

    work page 2007

  9. [9]

    Lions, Some Methods in the Mathematical Analysis of Systems and Their Control, Gordon And Breach, 1981

    J. Lions, Some Methods in the Mathematical Analysis of Systems and Their Control, Gordon And Breach, 1981

    work page 1981

  10. [10]

    G. Allaire, Homogenization of the unsteady stokes equations in porous media, in: Progress in Partial Differential Equations: Calculus of Variations, Applications, Longman Scientific & Technical, 1992, pp. 109–123

    work page 1992

  11. [11]

    Mikelic, Mathematical derivation of the Darcy-type law with memory effects, governing transient flow through porous medium, Glasnik Matematicki 29 (49) (1994) 57–77

    A. Mikelic, Mathematical derivation of the Darcy-type law with memory effects, governing transient flow through porous medium, Glasnik Matematicki 29 (49) (1994) 57–77

    work page 1994

  12. [12]

    G. V. Sandrakov, Homogenization of non-stationary Stokes equations with viscosity in a perforated domain, Izvestiya: Mathematics 61 (1) (1997) 113–141

    work page 1997

  13. [13]

    Knabner, L

    P. Knabner, L. Angermann, Numerical Methods for Elliptic and Parabolic Partial Differential Equa- tions, Springer, New York, 2003

    work page 2003

  14. [14]

    Quarteroni, A

    A. Quarteroni, A. Valli, Numerical Approximation of Partial Differential Equations, Springer-Verlag, Berlin, 1994. 19

    work page 1994

  15. [15]

    C.Chen, T.Shih, FiniteElementMethodsforIntegrodifferentialEquations, WorldScientific, Singapore, 1998

    work page 1998

  16. [16]

    Lubich, Convolution quadrature and discretized operational calculus

    C. Lubich, Convolution quadrature and discretized operational calculus. I, Numerische Mathematik 52 (2) (1988) 129–145

    work page 1988

  17. [17]

    B. Jin, Z. Zhou, Numerical Treatment and Analysis of Time-Fractional Evolution Equations, Springer, 2023

    work page 2023

  18. [18]

    Linz, Analytical and Numerical Methods for Volterra Equations, SIAM, 1985

    P. Linz, Analytical and Numerical Methods for Volterra Equations, SIAM, 1985

    work page 1985

  19. [19]

    Vabishchevich, Numerical solution of theCauchy problemfor Volterraintegrodifferentialequations with difference kernels, Applied Numerical Mathematics 174 (2022) 177–190

    P.N. Vabishchevich, Numerical solution of theCauchy problemfor Volterraintegrodifferentialequations with difference kernels, Applied Numerical Mathematics 174 (2022) 177–190

    work page 2022

  20. [20]

    P. N. Vabishchevich, Approximate solution of the Cauchy problem for a first-order integrodifferential equation with solution derivative memory, Journal of Computational and Applied Mathematics 422 (2023) 114887

    work page 2023

  21. [21]

    D. G. Roychowdhury, Computational Fluid Dynamics for Incompressible Flows, CRC Press, 2020

    work page 2020

  22. [22]

    Glowinski, T.-W

    R. Glowinski, T.-W. Pan, Numerical Simulation of Incompressible Viscous Flow: Methods and Appli- cations, Walter de Gruyter GmbH, 2022

    work page 2022

  23. [23]

    Meinardus, Approximation of Functions: Theory and Numerical Methods, Springer, Berlin, Heidel- berg, 2012

    G. Meinardus, Approximation of Functions: Theory and Numerical Methods, Springer, Berlin, Heidel- berg, 2012

    work page 2012

  24. [24]

    L. N. Trefethen, Approximation Theory and Approximation Practice, Extended Edition, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2019

    work page 2019

  25. [25]

    Pereyra, G

    V. Pereyra, G. Scherer (Eds.), Exponential Data Fitting and Its Applications, Bentham Science Pub- lishers, Sharjah, UAE, 2010

    work page 2010

  26. [26]

    Holmström, J

    K. Holmström, J. Petersson, A review of the parameter estimation problem of fitting positive exponen- tial sums to empirical data, Applied Mathematics and Computation 126 (1) (2002) 31–61

    work page 2002

  27. [27]

    P.N.Vabishchevich, A.N.Semenov, Computationalalgorithmfornonlinearapproximationoffunctions, Moscow University Computational Mathematics and Cybernetics 50 (1) (2026) 25–31

    work page 2026

  28. [28]

    Braess, Nonlinear Approximation Theory, Springer, Berlin, Heidelberg, 1986

    D. Braess, Nonlinear Approximation Theory, Springer, Berlin, Heidelberg, 1986

    work page 1986

  29. [29]

    Nakatsukasa, L

    Y. Nakatsukasa, L. N. Trefethen, Applications of AAA rational approximation, arXiv preprint (2025). arXiv:2510.16237

    work page arXiv 2025

  30. [30]

    P. N. Vabishchevich, Computational homogenization of parabolic equations with memory effects for a periodic heterogeneous medium, Journal of Computational and Applied Mathematics 477 (2026) 1–16

    work page 2026

  31. [31]

    L. Wang, Q. Xu, Z. Zhang, Corrector estimates and homogenization errors of unsteady flow ruled by Darcy’s law, Mathematische Annalen (2025) 1–65

    work page 2025

  32. [32]

    Galdi, An Introduction to the Mathematical Theory of the Navier-Stokes Equations: Steady-State Problems, Springer Science & Business Media, 2011

    G. Galdi, An Introduction to the Mathematical Theory of the Navier-Stokes Equations: Steady-State Problems, Springer Science & Business Media, 2011

    work page 2011

  33. [33]

    Babuška, J

    I. Babuška, J. E. Osborn, Eigenvalue problems, Handbook of Numerical Analysis II (1991) 641–787

    work page 1991

  34. [34]

    Boffi, F

    D. Boffi, F. Brezzi, M. Fortin, Mixed Finite Element Methods and Applications, Springer, 2013

    work page 2013

  35. [35]

    Saad, Numerical Methods for Large Eigenvalue Problems, SIAM, 2011

    Y. Saad, Numerical Methods for Large Eigenvalue Problems, SIAM, 2011

    work page 2011

  36. [36]

    A. A. Samarskii, The Theory of Difference Schemes, Marcel Dekker, New York, 2001

    work page 2001

  37. [37]

    A. A. Samarskii, P. P. Matus, P. N. Vabishchevich, Difference Schemes with Operator Factors, Kluwer Academic, Dordrecht, 2002. 20

    work page 2002