pith. sign in

arxiv: 2512.18787 · v1 · submitted 2025-12-21 · 🧮 math.AP

On the effects of surface roughness in non-isothermal porous medium flow

Mar\'ia Anguiano , Igor Pa\v{z}anin , Francisco J. Su\'arez-Grau This is my paper

Pith reviewed 2026-05-16 20:59 UTC · model grok-4.3

classification 🧮 math.AP
keywords surface roughnessnon-isothermal flowporous mediumDarcy-Brinkmanthin-film flowasymptotic analysisperiodic unfoldingviscous dissipation
0
0 comments X

The pith

In the critical roughness regime, non-isothermal thin-film porous flow converges to a strongly coupled effective model driven by the oscillatory boundary geometry.

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

The paper establishes convergence of velocity, pressure, and temperature fields in a Darcy-Brinkman thin-film flow that includes viscous dissipation as a heat source and a periodically oscillating boundary. As the combined small parameter for film thickness and roughness period tends to zero, the limit systems take the form of coupled elliptic equations whose structure depends on the relative scaling between roughness wavelength and film thickness. In the critical scaling, the oscillatory geometry forces a strong coupling between the flow and temperature equations that is absent when the boundary is smooth. This matters for applications in which rough porous surfaces interact with heat generation, because the derived effective models replace the need to resolve every small oscillation with simpler macroscopic equations that still capture the geometry-induced interaction.

Core claim

Using asymptotic analysis and the periodic unfolding method, the velocity, pressure, and temperature fields converge to solutions of limit problems that combine Reynolds-type equations with Darcy-Brinkman cell problems and a reduced energy equation. The degree of coupling between these equations is controlled by the relative scaling of the roughness wavelength and the film thickness; the critical regime produces a strong geometry-induced coupling that does not appear in the corresponding smooth-boundary case.

What carries the argument

Periodic unfolding method applied to the non-isothermal Darcy-Brinkman equations on a domain with an oscillating boundary, yielding scaling-dependent limit systems that encode roughness effects through cell problems.

Load-bearing premise

The boundary is periodically oscillating and the small parameter combining film thickness and roughness period tends to zero under the stated relative scalings that permit the asymptotic expansions and periodic unfolding.

What would settle it

A numerical simulation of the original non-isothermal Darcy-Brinkman equations on a domain with small but finite periodic roughness in the critical scaling regime that produces velocity, pressure, and temperature fields failing to approach the predicted coupled limit system would falsify the convergence result.

read the original abstract

We analyze a non-isothermal Darcy-Brinkman thin-film flow with a periodically oscillating boundary and viscous dissipation acting as a heat source. Using asymptotic analysis and the periodic unfolding method, we establish the convergence of velocity, pressure, and temperature fields as the small parameter (related to the film thickness and the period of the roughness) tends to zero. The limit problems depend on the relative scaling of the roughness wavelength and consist of coupled elliptic systems combining Reynolds-type equations with Darcy-Brinkman cell problems and reduced energy equation. In the critical roughness regime, the effective model exhibits a strong coupling induced by the oscillatory geometry, which does not occur in a smooth-boundary case.

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 / 3 minor

Summary. The manuscript analyzes non-isothermal Darcy-Brinkman thin-film flow over a periodically oscillating boundary with viscous dissipation as heat source. Asymptotic expansions combined with periodic unfolding are used to pass to the limit as the small parameter (film thickness and roughness period) tends to zero under different relative scalings. The resulting effective models are coupled elliptic systems consisting of a Reynolds-type pressure equation, Brinkman cell problems, and a reduced energy balance; in the critical roughness regime these systems contain additional geometry-induced coupling terms that are absent for a flat boundary.

Significance. If the convergence statements hold, the paper supplies rigorous effective equations that quantify how periodic roughness induces strong cross-coupling between flow and temperature fields in the critical scaling regime. This extends standard homogenization tools to a non-isothermal setting with viscous heating and provides explicit, geometry-dependent limit problems that can be used for reduced-order modeling in lubrication or porous-media heat transfer applications. The reliance on periodic unfolding and standard a-priori estimates is a methodological strength.

major comments (2)
  1. [§3.2, Theorem 3.1] §3.2, Theorem 3.1: the passage to the limit in the viscous-dissipation term inside the unfolded energy equation requires uniform integrability that is asserted but not verified in detail; the proof sketch relies on the a-priori bound (3.8) without showing that the dissipation term remains compact under the critical scaling.
  2. [§4.3, system (4.12)–(4.14)] §4.3, system (4.12)–(4.14): the cell problems for the Brinkman correctors are stated to be well-posed for temperature-dependent viscosity, yet the fixed-point argument used to close the coupled system is only sketched; an explicit contraction estimate or reference to a standard result for the Stokes-Brinkman operator with bounded coefficients would strengthen the claim.
minor comments (3)
  1. [§2.1] The definition of the unfolding operator U_ε in §2.1 is introduced without recalling its basic properties (e.g., the integral identity (2.4)); a short reminder would improve readability for readers unfamiliar with the method.
  2. [Figure 2] Figure 2 (critical-regime geometry) uses a shading that is difficult to distinguish in black-and-white print; a clearer line-style distinction between the oscillating boundary and the effective flat interface is recommended.
  3. [Theorem 5.1] The statement of the main convergence result (Theorem 5.1) lists the topologies but omits the precise statement of the boundary conditions satisfied by the limit temperature; this should be added for completeness.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comments. We address each major point below and will revise the manuscript accordingly to strengthen the proofs.

read point-by-point responses

  1. Referee: [§3.2, Theorem 3.1] §3.2, Theorem 3.1: the passage to the limit in the viscous-dissipation term inside the unfolded energy equation requires uniform integrability that is asserted but not verified in detail; the proof sketch relies on the a-priori bound (3.8) without showing that the dissipation term remains compact under the critical scaling.

    Authors: We agree that a more detailed verification of uniform integrability is needed for the viscous dissipation term under the critical scaling. In the revised manuscript we will add an explicit argument: starting from the a-priori bound (3.8) on the temperature and the uniform L^2 bound on the velocity gradients furnished by the Darcy-Brinkman cell problems, we apply Hölder's inequality together with the boundedness of the viscosity function to obtain a uniform L^{1+δ} bound on the dissipation for some δ>0. This implies uniform integrability via the de la Vallée-Poussin criterion, allowing passage to the limit inside the unfolded energy equation by Vitali's convergence theorem. revision: yes

  2. Referee: [§4.3, system (4.12)–(4.14)] §4.3, system (4.12)–(4.14): the cell problems for the Brinkman correctors are stated to be well-posed for temperature-dependent viscosity, yet the fixed-point argument used to close the coupled system is only sketched; an explicit contraction estimate or reference to a standard result for the Stokes-Brinkman operator with bounded coefficients would strengthen the claim.

    Authors: We thank the referee for this observation. In the revision we will replace the sketch with an explicit contraction estimate: let μ(θ) be Lipschitz continuous with constant L. For two temperatures θ1, θ2 we consider the difference of the corresponding Brinkman solutions and use the standard energy estimate for the linear Stokes-Brinkman operator with bounded coefficients (which yields a constant C independent of the temperature) to obtain ||u1−u2||_{H^1} ≤ C L ||θ1−θ2||_∞ ||u2||_{H^1}. Choosing the ball radius appropriately then shows that the fixed-point map is a contraction on a sufficiently small ball in C(Ω̄), thereby closing the argument rigorously. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivations rely on standard asymptotic expansions and periodic unfolding applied to the given non-isothermal Darcy-Brinkman system. Convergence is obtained via a-priori estimates, compactness arguments, and limit passage in the unfolded formulation. The geometry-induced coupling terms in the critical regime emerge directly from the oscillatory boundary under the stated scalings; no quantities are defined in terms of their own outputs, no fitted parameters are renamed as predictions, and no load-bearing steps reduce to self-citations or ansatzes imported from prior author work. The limit systems are independent of the target results and follow from the PDE structure and unfolding operator.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard PDE assumptions for the Darcy-Brinkman system with viscous dissipation and periodic boundary data; no free parameters, new entities or ad-hoc axioms are introduced beyond the usual homogenization framework.

axioms (1)
  • domain assumption Darcy-Brinkman equations govern the flow with viscous dissipation as the sole heat source and the upper boundary is periodically oscillating
    These form the starting microscopic model whose limit is derived.

pith-pipeline@v0.9.0 · 5416 in / 1109 out tokens · 28692 ms · 2026-05-16T20:59:18.783337+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

34 extracted references · 34 canonical work pages

  1. [1]

    Allaire,Homogenization of the Navier-Stokes equations in open sets perforated with thiny holes I

    A. Allaire,Homogenization of the Navier-Stokes equations in open sets perforated with thiny holes I. Abstract framework, a volume distribution of holes, Arch. Ration. Mech. Anal. 113 (1991), 209–259

    work page 1991

  2. [2]

    Al-Hadhrami, L

    A.K. Al-Hadhrami, L. Elliott, D.B. Ingham,A new model for viscous dissipation in porous media across a range of permeability values, Transp. Porous Med. 53 (2003), 117–122

    work page 2003

  3. [3]

    Anguiano, F.J

    M. Anguiano, F.J. Su´ arez-Grau,Nonlinear Reynolds equations for non-Newtonian thin-film fluid flows over a rough boundary, IMA J. Appl. Math. 84 (2019), 63–95

    work page 2019

  4. [4]

    Anguiano, F.J

    M. Anguiano, F.J. Su´ arez-Grau,Modeling of a micropolar thin film flow with rapidly varying thickness and non-standard boundary conditions, Acta Math. Sci. 46B(1) (2026), 209–242

    work page 2026

  5. [5]

    Bayada, M

    G. Bayada, M. Chambat,Homogenization of the Stokes system in a thin film flow with rapidly varying thickness, RAIRO Mod´ el. Math. Anal. Num´ er. 23 (1989) 205–234

    work page 1989

  6. [6]

    Beneˇ s, I

    M. Beneˇ s, I. Paˇ zanin,Rigorous derivation of the effective model describing a non-isothermal fluid flow through a vertical pipe filled with porous medium, Contin. Mech. Thermodyn. 30 (2018), 301–317

    work page 2018

  7. [7]

    Boukrouche, R

    M. Boukrouche, R. El Mir,On a non-isothermal, non-Newtonian lubrication problem with Tresca law: Existence and the behavior of weak solutions, Nonlinear Anal-Real World Appl. 9 (2008), 674–692

    work page 2008

  8. [8]

    Brinkman,A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles, Appl

    H. Brinkman,A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles, Appl. Sci. Res. A1 (1947), 27–34

    work page 1947

  9. [9]

    Casado-D´ ıaz, M

    J. Casado-D´ ıaz, M. Luna-Laynez, F.J. Su´ arez-Grau,A decomposition result for the pressure of a fluid in a thin domain and extensions to elasticity problems, SIAM J. Math. Anal. 52 (2020), 2201–2236

    work page 2020

  10. [10]

    Cioranescu, A

    D. Cioranescu, A. Damlamian, G. Griso,Periodic unfolding and homogenization, C.R. Acad. Sci. Paris Ser. I 335 (2002) 99–104

    work page 2002

  11. [11]

    Cioranescu, A

    D. Cioranescu, A. Damlamian, G. Griso,The periodic unfolding method in homogenization, SIAM J. Math. Anal. 40 (2008) 1585–1620

    work page 2008

  12. [12]

    Cioranescu, A

    D. Cioranescu, A. Damlamian, G. Griso,The periodic unfolding method: theory and applications to partial differential problems, Series in Contemporary Mathematics, 3, Springer, Singapore, 2018

    work page 2018

  13. [13]

    Ciuperca, E

    I.S. Ciuperca, E. Feireisl, M. Jai, A. Petrov,Transition to thermohydrodynamic lubrication problem, Q. Appl. Math. 75 (2017), 391–414

    work page 2017

  14. [14]

    Darcy,Les Fontaines Publiques de la ville de Dijon, Dalmont, Paris, 1856

    H. Darcy,Les Fontaines Publiques de la ville de Dijon, Dalmont, Paris, 1856

    work page

  15. [15]

    Gallou¨ et, R

    T. Gallou¨ et, R. Herbin, Existence of a solution to a coupled elliptic system, Appl. Math. Lett. 7 (2) (1994) 49–55

    work page 1994

  16. [16]

    Mikeli´ c,Remark on the result on homogenization in hydrodynamical lubrication by G

    A. Mikeli´ c,Remark on the result on homogenization in hydrodynamical lubrication by G. Bayada and M. Chambat, RAIRO Mod´ el. Math. Anal. Num´ er. 25 (1991) 363–370. 24 Mar´ ıa Anguiano, Igor Paˇ zanin and Francisco J. Su´ arez-Grau

    work page 1991

  17. [17]

    R. El Mir,Etude math´ ematique et analyse asymptotique de quelques probl` emes de lubrification par des fluides incompressibles essentiellement non-Newtoniens avec des conditions de non adh´ erence aux bords, Hal Archive Ouverte, tel-00011638, (2005)

    work page 2005

  18. [18]

    Levy,Fluid flow through an array of fixed particles, Int

    T. Levy,Fluid flow through an array of fixed particles, Int. J. Eng. Sci. 21 (1983), 11–23

    work page 1983

  19. [19]

    Magyari, D.A.S

    E. Magyari, D.A.S. Rees, B. Keller,Effect of viscous dissipation on the flow in fluid saturated porous media, In Handbook of Porous Media, 2nd Edition, (Vafai, K., ed.), pp. 373–406, Taylor & Francis, 2005

    work page 2005

  20. [20]

    Maruˇ si´ c-Paloka, I

    E. Maruˇ si´ c-Paloka, I. Paˇ zanin,Non-isothermal fluid flow through a thin pipe with cooling, Appl. Anal. 88 (4) (2009), 495–515

    work page 2009

  21. [21]

    Maruˇ si´ c-Paloka, I

    E. Maruˇ si´ c-Paloka, I. Paˇ zanin,Higher-order effective model desrcibing a non-isothermal thin-film flow, Int. J. Multiscale Comput. Eng. 16 (2018), 121–130

    work page 2018

  22. [22]

    Maruˇ si´ c-Paloka, I

    E. Maruˇ si´ c-Paloka, I. Paˇ zanin,On the thermal flow through a porous annular region, J. Eng. Math. 147, 9 (2014), 1–19

    work page 2014

  23. [23]

    Maruˇ si´ c-Paloka, I

    E. Maruˇ si´ c-Paloka, I. Paˇ zanin, S. Maruˇ si´ c,Comparison between Darcy and Brinkman laws in a fracture, Appl. Math. Comput. 218 (2012), 7538–7545

    work page 2012

  24. [24]

    Panasenko, K

    G. Panasenko, K. Pileckas,Multiscale analysis of viscous flows in thin tube structures, Birkh¨auser, 2024

    work page 2024

  25. [25]

    Paˇ zanin, M

    I. Paˇ zanin, M. Radulovi´ c,Effects of the viscous dissipation on the Darcy-Brinkman flow: Rigorous derivation of the higher-order asymptotic model, Appl. Math. Comput. 386 (2020) 125479, pp. 1–12

    work page 2020

  26. [26]

    Paˇ zanin, F.J

    I. Paˇ zanin, F.J. Su´ arez-Grau,Homogenization of the Darcy-Lapwood-Brinkman flow through a thin domain with highly oscillating boundaries, Bull. Malays. Math. Sci. Soc. 42 (2019), 3073–3109

    work page 2019

  27. [27]

    Paˇ zanin, F.J

    I. Paˇ zanin, F.J. Su´ arez-Grau,Roughness-induced effects on the thermomicropolar fluid flow through a thin domain, Stud. Appl. Math. 151 (2023), 716–751

    work page 2023

  28. [28]

    Maruˇ si´ c-Paloka, I

    E. Maruˇ si´ c-Paloka, I. Paˇ zanin,Asymptotic modelling of the non-isothermal fluid flow through a porous annulus in the presence of strong magnetic field, Z. Angew. Math. Mech. 105 (11) (2025), pp. 1–17

    work page 2025

  29. [29]

    Paˇ zanin, F.J

    I. Paˇ zanin, F.J. Su´ arez-Grau,Analysis of the Darcy-Brinkman flow with viscous dissipation and non- homogeneous thermal boundary condition, to appear in Math. Model. Anal (2025), pp. 1–23

    work page 2025

  30. [30]

    Sanchez-Palencia,On the asymptotics of the fluid flow past an array of fixed obstacles, Int

    E. Sanchez-Palencia,On the asymptotics of the fluid flow past an array of fixed obstacles, Int. J. Eng. Sci. 20 (1982), 1291–1301

    work page 1982

  31. [31]

    Su´ arez-Grau,Asymptotic behavior of a non-Newtonian flow in a thin domain with Navier law on a rough boundary, Nonlin

    F.J. Su´ arez-Grau,Asymptotic behavior of a non-Newtonian flow in a thin domain with Navier law on a rough boundary, Nonlin. Anal., 117 (2015) 99–123

    work page 2015

  32. [32]

    Su´ arez-Grau,Analysis of the roughness regime for micropolar fluids via homogenization, Bull

    F.J. Su´ arez-Grau,Analysis of the roughness regime for micropolar fluids via homogenization, Bull. Malays. Math. Sci. Soc., 44 (2021) 1613–1652

    work page 2021

  33. [33]

    Szeri,Fluid film lubrication: theory and design, Cambridge University Press, 1998

    A.Z. Szeri,Fluid film lubrication: theory and design, Cambridge University Press, 1998

    work page 1998

  34. [34]

    Tabeling,Introduction to microfluidics, 2nd Edition, Oxford University Press, 2023

    P. Tabeling,Introduction to microfluidics, 2nd Edition, Oxford University Press, 2023. 25

    work page 2023