Pith. sign in

REVIEW 2 major objections 2 minor 37 references

Computational homogenization of unsteady flows in a periodic porous medium

T0 review · 2 major / 2 minor · reviewed 2026-05-07 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗

Pith's one-line read Approximating the memory kernel as a sum of exponentials converts the nonlocal macroscale model of unsteady porous flows into a local differential system.

desk verdict The paper gives a concrete numerical recipe for unsteady porous-media homogenization by pulling a memory kernel from cell problems and approximating it with exponentials to recover a local system. read the letter →

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

P.N. Vabishchevich This is my paper · ORCID
classification math.NAcs.NA
keywords computationalhomogenizationporousmediaunsteadyflowsmemorykernelDarcylawfiniteelementmethodnumericalapproximation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

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.

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.

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.

Watch

Extended reading notes

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.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

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 · score 0.0 of 10

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.

Assumptions & free parameters 1 free parameters · 2 assumptions · 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.
assumptions (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.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Computational homogenization of unsteady flows in a periodic porous medium." pith.science (2026). https://pith.science/paper/2604.25722

@misc{pith2026260425722,
  author       = {Pith},
  title        = {Pith review of: Computational homogenization of unsteady flows in a periodic porous medium},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2604.25722}},
  note         = {Machine review of arXiv:2604.25722}
}
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.

Figures

Figures reproduced from arXiv: 2604.25722 by the authors.

Figure 1
Figure 1. Periodic porous medium and periodicity cell view at source ↗
Figure 2
Figure 2. Auxiliary functions w0 1 · e1 (left) and w0 2 · e1 (right) view at source ↗
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,015 nodes. … view at source ↗
Figures from the paper (6 more)
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
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 |
Figure 6
Figure 6. Figure 6: Coefficients |a k 1 | (top) and |a k 2 |, k = 1, 2, . . . , m (bottom). 16
Figure 7
Figure 7. Figure 7: Elements of the steady-state conductivity tensor: diagonal (top) and off-diagonal (bottom).
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…
Figure 9
Figure 9. Figure 9: Solution of the macroscale problem taking into account memory effects at different time instants.

Discussion (0). Continue with ORCID to comment.

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

Show all 37 references
  1. [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

    1981

  2. [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

    1992

  3. [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

    1994

  4. [12]

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

    1997

  5. [13]

    Knabner, L

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

    2003

  6. [14]

    Quarteroni, A

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

    1994

  7. [15]

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

    1998

  8. [16]

    Lubich, Convolution quadrature and discretized operational calculus

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

    1988

  9. [17]

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

    2023

  10. [18]

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

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

    1985

  11. [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

    2022

  12. [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

    2023

  13. [21]

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

    2020

  14. [22]

    Glowinski, T.-W

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

    2022

  15. [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

    2012

  16. [24]

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

    2019

  17. [25]

    Pereyra, G

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

    2010

  18. [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

    2002

  19. [27]

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

    2026

  20. [28]

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

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

    1986

  21. [29]

    Nakatsukasa, L

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

    2025

  22. [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

    2026

  23. [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

    2025

  24. [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

    2011

  25. [33]

    Babuška, J

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

    1991

  26. [34]

    Boffi, F

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

    2013

  27. [35]

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

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

    2011

  28. [36]

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

    2001

  29. [37]

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

    2002

Pith tools

Reviewed May 7, 2026 · model on record in the stance chip above.