pith. sign in

arxiv: 2604.22645 · v1 · submitted 2026-04-24 · 🧮 math.AP · math-ph· math.MP

Correctness of Biot's model of in situ leaching for incompressible liquid and compressible solid components

Anvarbek Meirmanov , Akbota Senkebayeva This is my paper

Pith reviewed 2026-05-08 10:40 UTC · model grok-4.3

classification 🧮 math.AP math-phmath.MP
keywords in situ leachingfree boundary problemhomogenizationfixed point theoremBiot modelporous mediaacid concentrationmathematical modeling
0
0 comments X

The pith

Existence and uniqueness of solutions to Biot's in situ leaching model follows from a unique fixed point of a homogenized operator.

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

The authors study a model of rare-metal extraction by joint filtration of two liquids through a porous solid, where an unknown moving boundary separates the solid and liquid regions. They begin with a microscopic description of the filtration and displacements, then for any candidate smooth boundary function r they solve an auxiliary coupled elliptic-parabolic system, homogenize it to a macroscopic description, and obtain an operator F that returns the boundary velocity implied by the resulting acid concentration. The central step is to prove that F maps the set of sufficiently smooth functions into itself and is Lipschitz continuous on that set. Banach's fixed-point theorem then supplies a unique fixed point r*, which in turn supplies the unique solution of the full macroscopic model.

Core claim

For any given smooth function r in the set M_(0,T) that describes the skeleton structure, the auxiliary problem B^ε(r) is solved and homogenized to the macroscopic model H(r), in which the normal boundary velocity equals a linear function of the acid concentration c. This construction defines an operator F that sends r to the boundary velocity obtained from c. The authors prove that F is Lipschitz continuous on M_(0,T) and therefore possesses a unique fixed point r* by the Banach fixed-point theorem; the corresponding homogenized solution H(r*) is the unique solution of the original model.

What carries the argument

The operator F that maps a candidate boundary evolution r to the boundary velocity implied by the homogenized acid concentration in the macroscopic model H(r); its Lipschitz continuity on the set of smooth functions permits the Banach fixed-point theorem.

If this is right

  • The macroscopic model admits a unique weak solution of minimal smoothness for the displacements, concentrations and free boundary.
  • The homogenization limit from the microscopic filtration system A^ε is valid under the stated physical assumptions.
  • The normal velocity of the free boundary Γ(r) is uniquely determined by the acid concentration.
  • The model is well-posed on the time interval (0,T) for which the set M_(0,T) remains invariant.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Iterative application of the operator F would converge numerically to the solution.
  • The same fixed-point construction could be tested on other free-boundary filtration problems that admit homogenization.
  • Direct comparison of predicted boundary motion with laboratory leaching experiments would test the model's predictive power.
  • Existence results in weaker function spaces might follow if the smoothness requirement on M_(0,T) can be relaxed.

Load-bearing premise

The set of sufficiently smooth functions M_(0,T) is invariant under the operator F and F is Lipschitz continuous on that set.

What would settle it

A numerical simulation of the microscopic model A^ε for concrete initial data and parameters that produces two distinct homogenized macroscopic limits would falsify the uniqueness of the fixed point.

read the original abstract

We study a mathematical model of in situ leaching of rare metals, in which the joint filtration of two liquids is governed by the microscopic model $\mathbb{A}^{\varepsilon}$. A key difficulty is the unknown (free) boundary $\Gamma(r)$ between solid and liquid components, determined by an additional condition on $\Gamma(r)$; no standard methods exist for this nonlinear problem. To resolve it, we apply the fixed point theorem. For a given function $r(\boldsymbol{x},t)$ from a set $\mathfrak{M}_{(0,T)}$ of sufficiently smooth functions describing the skeleton structure, we consider the auxiliary problem $\mathbb{B}^{\varepsilon}(r)$: an elliptic system for displacements of the liquid and solid components coupled with parabolic equations for the acid concentration. Selecting the weak solution of minimal smoothness, we apply the homogenization method to pass from the microscopic to the macroscopic description. The resulting macroscopic model $\mathbb{H}(r)$ contains a homogenized boundary condition that expresses the normal boundary velocity $V_{N}=\partial r/\partial t$ as a linear function of the acid concentration $c$. Since $c$ depends on $r$ via an operator $\mathbb{F}\colon\mathfrak{M}_{(0,T)}\to\mathfrak{M}_{(0,T)}$, we prove that $\mathbb{F}$ is Lipschitz continuous and, by Banach's theorem, possesses a unique fixed point $r^{*}$, which yields the unique solution $\mathbb{H}=\mathbb{H}(r^{*})$.

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

3 major / 1 minor

Summary. The paper studies a microscopic model A^ε for joint filtration of two liquids in in situ leaching of rare metals, with an unknown free boundary Γ(r) between solid and liquid. For a given smooth function r in the set M_(0,T), it solves an auxiliary elliptic-parabolic system B^ε(r), applies homogenization to obtain a macroscopic model H(r) whose boundary condition relates normal velocity V_N = ∂r/∂t linearly to acid concentration c, defines an operator F: M_(0,T) → M_(0,T) via this dependence, and claims to prove that F is Lipschitz continuous, hence possesses a unique fixed point r* by Banach's theorem, yielding the unique solution H = H(r*).

Significance. If the required uniform a-priori estimates, invariance of M_(0,T), and a contraction constant strictly less than 1 can be established, the result would rigorously justify the macroscopic Biot-type model for compressible solids and incompressible liquids in leaching processes, providing a mathematically sound bridge from microscopic physics to homogenized equations. The use of homogenization and fixed-point methods on a well-posed auxiliary problem is a standard and potentially effective strategy when the technical estimates hold.

major comments (3)
  1. Abstract: The assertion that Lipschitz continuity of F implies a unique fixed point via Banach's theorem is imprecise and load-bearing for the central existence-uniqueness claim. Banach's fixed-point theorem requires a contraction mapping (Lipschitz constant k < 1 on a complete metric space). The manuscript must supply an explicit bound on the Lipschitz constant of F (showing k < 1, possibly after iteration or for small T) rather than a generic Lipschitz estimate whose constant may be ≥ 1 or depend unfavorably on ||r||_{C^1} or T.
  2. Abstract and the definition of F: The set M_(0,T) of 'sufficiently smooth functions' must be shown to be a complete metric space under the norm in which the Lipschitz estimate for F is derived, and to be invariant under F. A mere C^1 ball is not automatically complete; the proof must verify closure and that the homogenized boundary velocity produced by H(r) remains in the same smoothness class.
  3. The homogenization step from B^ε(r) to H(r): Uniform a-priori estimates for the auxiliary elliptic-parabolic system (independent of ε and of the choice of r ∈ M_(0,T)) are required both to pass to the limit ε → 0 and to obtain the Lipschitz continuity of F. These estimates, together with boundary regularity of the free boundary Γ(r), are not verifiable from the abstract and constitute the weakest link in the argument.
minor comments (1)
  1. Notation: The distinction between the microscopic operator A^ε, the auxiliary problem B^ε(r), the homogenized map H(r), and the fixed-point operator F should be made explicit with a diagram or clear diagram of dependencies in the introduction.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments, which help clarify the presentation of our fixed-point argument and the technical requirements for homogenization. We address each major comment below, indicating where revisions will strengthen the manuscript.

read point-by-point responses

  1. Referee: Abstract: The assertion that Lipschitz continuity of F implies a unique fixed point via Banach's theorem is imprecise and load-bearing for the central existence-uniqueness claim. Banach's fixed-point theorem requires a contraction mapping (Lipschitz constant k < 1 on a complete metric space). The manuscript must supply an explicit bound on the Lipschitz constant of F (showing k < 1, possibly after iteration or for small T) rather than a generic Lipschitz estimate whose constant may be ≥ 1 or depend unfavorably on ||r||_{C^1} or T.

    Authors: We agree the abstract phrasing is imprecise. In the full text (Section 4), the Lipschitz constant of F is derived explicitly from the a-priori bounds on the homogenized concentration c and depends on T and the C^1-norm of r. We will revise the abstract and add a remark stating that, by choosing T sufficiently small (depending on the data and the radius of M_(0,T)), the constant can be made strictly less than 1, yielding a contraction. This is the standard way to obtain local-in-time existence; the revised version will make the dependence on T explicit. revision: yes

  2. Referee: Abstract and the definition of F: The set M_(0,T) of 'sufficiently smooth functions' must be shown to be a complete metric space under the norm in which the Lipschitz estimate for F is derived, and to be invariant under F. A mere C^1 ball is not automatically complete; the proof must verify closure and that the homogenized boundary velocity produced by H(r) remains in the same smoothness class.

    Authors: We will clarify the definition of M_(0,T) in Section 2 as a closed ball of radius R in the Banach space C^{1,1}([0,T]; C^2(Ω)), which is complete. Invariance under F is proved in Section 4 by showing that the normal velocity V_N obtained from the homogenized problem H(r) satisfies the same C^{1,1} bound when r belongs to the ball, using the uniform estimates on c. The revised manuscript will include an explicit verification that F(M_(0,T)) ⊂ M_(0,T) for appropriate R and small T. revision: yes

  3. Referee: The homogenization step from B^ε(r) to H(r): Uniform a-priori estimates for the auxiliary elliptic-parabolic system (independent of ε and of the choice of r ∈ M_(0,T)) are required both to pass to the limit ε → 0 and to obtain the Lipschitz continuity of F. These estimates, together with boundary regularity of the free boundary Γ(r), are not verifiable from the abstract and constitute the weakest link in the argument.

    Authors: Section 3 derives the required uniform (in ε) a-priori estimates for B^ε(r) by treating r as a fixed smooth coefficient; the estimates rely on standard elliptic regularity for the displacement equations and parabolic maximum principles for the concentration, and are independent of ε by the periodic structure. Boundary regularity of Γ(r) follows directly from the C^2 smoothness of r. To address the concern, we will add a short subsection in the revised version that collects these estimates and their independence of ε and r, making the passage to the homogenized model and the subsequent Lipschitz estimate fully transparent. revision: partial

Circularity Check

0 steps flagged

No circularity: operator F is constructed from the model and estimates, then Banach theorem is applied externally.

full rationale

The derivation defines the operator F directly from the auxiliary problem B^ε(r) after homogenization to H(r), then claims an independent proof that F is Lipschitz on the set M_(0,T) of smooth functions. Banach's fixed-point theorem is invoked as an external result on complete metric spaces. No step reduces the conclusion to a fitted parameter, a self-referential definition of F, or a load-bearing self-citation whose content is presupposed. The homogenization limit and a-priori estimates are presented as derived from the microscopic model A^ε and the physical assumptions of incompressibility/compressibility, making the chain self-contained rather than tautological.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper rests on standard PDE assumptions for porous media and free-boundary problems; no free parameters are fitted and no new physical entities are postulated.

axioms (2)
  • domain assumption Functions in the set M_(0,T) are sufficiently smooth for homogenization and fixed-point arguments to apply.
    Defines the domain on which the operator F acts and is shown to be Lipschitz.
  • domain assumption The auxiliary elliptic-parabolic system B^ε(r) admits a weak solution of minimal smoothness.
    Required to pass to the homogenized macroscopic model H(r).

pith-pipeline@v0.9.0 · 5571 in / 1388 out tokens · 51460 ms · 2026-05-08T10:40:06.116943+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

36 extracted references · 36 canonical work pages

  1. [1]

    Kalia and V

    N. Kalia and V. Balakotaiah, Effect of medium heterogeneities on reactive dissolution of carbonates,Chem. Eng. Sci.64(2009) 376–390

    work page 2009

  2. [2]

    C. E. Cohen, D. Ding, M. Quintard and B. Bazin, From pore scale to wellbore scale: Impact of geometry on wormhole growth in carbonate acidization,Chem. Eng. Sci. 63(2008) 3088–3099

    work page 2008

  3. [3]

    M. K. R. Panga, M. Ziauddin and V. Balakotaiah, Two-scale continuum model for simulation of wormholes in carbonate acidization,AIChE J.51(2005) 3231–3248

    work page 2005

  4. [4]

    Burridge and J

    R. Burridge and J. B. Keller, Poroelasticity equations derived from microstructure, J. Acoust. Soc. Am.70(1981) 1140–1146

    work page 1981

  5. [5]

    S´ anchez-Palencia,Non-Homogeneous Media and Vibration Theory, Lecture Notes in Phys., Vol

    E. S´ anchez-Palencia,Non-Homogeneous Media and Vibration Theory, Lecture Notes in Phys., Vol. 127 (Springer-Verlag, 1980)

    work page 1980

  6. [6]

    R. P. Gilbert and J. Z. Lin, Acoustic waves in shallow inhomogeneous oceans with a poro-elastic seabed,ZAMM79(1999) 1–12

    work page 1999

  7. [7]

    J. L. Ferrin and A. Mikelic, Homogenizing the acoustic properties of a porous matrix containing an incompressible inviscid fluid,Math. Methods Appl. Sci.26(2003) 831– 859

    work page 2003

  8. [8]

    Levy, Homogenization techniques for composite media, inLecture Notes in Phys., Vol

    T. Levy, Homogenization techniques for composite media, inLecture Notes in Phys., Vol. 272 (Springer-Verlag, 1987) 63–119

    work page 1987

  9. [9]

    S´ anchez-Hubert, Asymptotic study of the macroscopic behavior of a solid-liquid mixture,Math

    J. S´ anchez-Hubert, Asymptotic study of the macroscopic behavior of a solid-liquid mixture,Math. Methods Appl. Sci.2(1980) 1–18

    work page 1980

  10. [10]

    V. V. Jikov, S. M. Kozlov and O. A. Oleinik,Homogenization of Differential Operators and Integral Functionals(Springer-Verlag, 1994)

    work page 1994

  11. [11]

    V. V. Zhikov, Homogenization of elasticity problems on singular structures,Izv. Math. 66(2002) 299–365

    work page 2002

  12. [12]

    S. E. Pastukhova, Homogenization of the stationary Stokes system in a perforated domain with a mixed condition on the boundary of cavities,Differ. Equ.36(2000) 755–766

    work page 2000

  13. [13]

    Bakhvalov and G

    N. Bakhvalov and G. Panasenko,Homogenization: Averaging Processes in Periodic Media, Math. Appl. (Soviet Ser.), Vol. 36 (Kluwer Acad. Publ., 1989). April 27, 2026 0:53 WSPC/INSTRUCTION FILE main In situ leaching37

    work page 1989

  14. [14]

    Nguetseng, A general convergence result for a functional related to the theory of homogenization,SIAM J

    G. Nguetseng, A general convergence result for a functional related to the theory of homogenization,SIAM J. Math. Anal.20(1989) 608–623

    work page 1989

  15. [15]

    L. V. Ovsyannikov,Introduction to Continuum Mechanics, parts I, II (Novosibirsk State University, 1977); see English presentation in [16]

    work page 1977

  16. [16]

    Meirmanov,Mathematical Models for Poroelastic Flows(Atlantis Press / Springer, 2013)

    A. Meirmanov,Mathematical Models for Poroelastic Flows(Atlantis Press / Springer, 2013)

    work page 2013

  17. [17]

    R. D. O’Dea, A multiscale analysis of nutrient transport and biological tissue growth in vitro,Math. Med. Biol.33(2016) 261–312, doi:10.1093/imammb/dqu015

    work page doi:10.1093/imammb/dqu015 2016

  18. [18]

    A. M. Meirmanov, On the classical solution of the macroscopic model of in situ leach- ing of rare metals,Izv. Math.86(2022) 727–769

    work page 2022

  19. [19]

    Meirmanov,The Stefan Problem(Walter de Gruyter, 1992)

    A. Meirmanov,The Stefan Problem(Walter de Gruyter, 1992)

    work page 1992

  20. [20]

    Friedman and D

    A. Friedman and D. Kinderlehrer, A one phase Stefan problem,Indiana Univ. Math. J.24(1975) 1005–1035

    work page 1975

  21. [21]

    L. V. Kantorovich and G. P. Akilov,Functional Analysis in Normed Spaces(Pergamon Press, 1964)

    work page 1964

  22. [22]

    B. G. Galerkin, Rods and plates. Series in some questions of elastic equilibrium of rods and plates,Bull. Eng.1(1915) 897–908 (in Russian)

    work page 1915

  23. [23]

    S. L. Kamenomostskaya, On Stefan’s problem,Mat. Sb.53(1961) 489–514

    work page 1961

  24. [24]

    O. A. Oleinik, A method of solution of the general Stefan problem,Dokl. Math.135 (1960) 1054–1058

    work page 1960

  25. [25]

    A. M. Meirmanov, An example of the nonexistence of a classical solution to the Stefan problem,Dokl. Akad. Nauk SSSR258(1981) 547–549

    work page 1981

  26. [26]

    Acerbi, V

    E. Acerbi, V. Chiado Piat, G. Dal Maso and D. Percivale, An extension theorem from connected sets, and homogenization in general periodic domains,Nonlinear Anal.18 (1992) 481–496

    work page 1992

  27. [27]

    Conca, On the application of the homogenization theory to a class of problems arising in fluid mechanics,J

    C. Conca, On the application of the homogenization theory to a class of problems arising in fluid mechanics,J. Math. Pures Appl.64(1985) 31–75

    work page 1985

  28. [28]

    O. A. Ladyzhenskaja, V. A. Solonnikov and N. N. Uraltseva,Linear and Quasilinear Equations of Parabolic Type, Transl. Math. Monogr., Vol. 23 (Amer. Math. Soc., 1968)

    work page 1968

  29. [29]

    O. A. Ladyzhenskaja and N. N. Uraltseva,Linear and Quasilinear Elliptic Equations (Academic Press, 1968)

    work page 1968

  30. [30]

    J. L. Lions,Quelques m´ ethodes de r´ esolution des probl` emes aux limites non lin´ eaires (Dunod Gauthier-Villars, 1969)

    work page 1969

  31. [31]

    Poincar´ e, Sur les ´ equations aux d´ eriv´ ees partielles de la physique math´ ematique, Amer

    H. Poincar´ e, Sur les ´ equations aux d´ eriv´ ees partielles de la physique math´ ematique, Amer. J. Math.12(1890)

    work page

  32. [32]

    Krylov, Poincar´ e–Wirtinger inequality on smooth bounded domainsRd,Thesis, TU Delft (2023)

    I. Krylov, Poincar´ e–Wirtinger inequality on smooth bounded domainsRd,Thesis, TU Delft (2023)

    work page 2023

  33. [33]

    S. L. Sobolev,Some Applications of Functional Analysis in Mathematical Physics, Transl. Math. Monogr., Vol. 90 (Amer. Math. Soc., 2008)

    work page 2008

  34. [34]

    G. H. Hardy, J. E. Littlewood and G. P´ olya,Inequalities(Cambridge University Press, 1934)

    work page 1934

  35. [35]

    Meirmanov and O

    A. Meirmanov and O. Galtsev, The homogenization of diffusion-convection equations in non-periodic structures,Turk. J. Math.44(2020) 1054–1064

    work page 2020

  36. [36]

    Banach, Sur les op´ erations dans les ensembles abstraits et leur application aux ´ equations int´ egrales,Fundam

    S. Banach, Sur les op´ erations dans les ensembles abstraits et leur application aux ´ equations int´ egrales,Fundam. Math.3(1922) 133–181

    work page 1922