A mathematical preprint has derived a closed analytical solution for water infiltrating a porous material while a hydration reaction consumes some of the water. Under the model’s stated conditions, the equations have a unique solution for all time. The formulas describe the position of the infiltration front, the hydration field and the pressure, while showing that the front moves strictly forward. The result is a description of what this idealized mathematical system does, rather than an experimental measurement.
The model follows water entering a homogeneous porous medium along a one-dimensional horizontal domain. A fixed pressure head is imposed at the origin, and the material is divided by a sharp boundary between the wet and dry regions. As water advances, a first-order hydration reaction consumes infiltrating water. The central question is how that moving boundary, written as s(t), changes over time when the reaction is included.
A solvable moving boundary
To obtain an explicit result, the authors formulate the problem in dimensionless form as a quasi-steady free-boundary problem. In ordinary language, the boundary is “free” because its location is not prescribed in advance; it is part of the solution. The formulation combines pressure boundary conditions with a Stefan-type mass-conservation condition at the infiltration front, linking the front’s motion to the water balance.
The derivation introduces an auxiliary function called J and uses it to reduce the moving-boundary calculation to an ordinary differential equation. That reduction supplies the route to the exact solution and to the existence argument. The resulting theorem states that, when porosity is positive and Y is real-valued, the model has one global-in-time solution rather than multiple competing mathematical solutions.
The solution is expressed in closed formulas for the front position s, the hydration field Γ and the pressure p. Because the free boundary is strictly increasing, the model does not produce a retreating infiltration front under the stated conditions. The explicit form also allows the analysis to examine the front at both early and late times, including the points where its time dependence changes.
The front keeps changing its pace
The first result is a square-root law at early times. The front position grows in proportion to the square root of time, and the coefficient in that early-time behavior depends only on porosity. Within the model, this makes the initial scaling universal across the hydration settings considered: the reaction does not alter the leading square-root form at the start of infiltration.
The no-hydration case provides the simplest benchmark. When Y = 0, the infiltration front follows a single square-root law for all positive times, not only during the early-time approximation. In that parameter case, the same basic time dependence continues throughout the modeled process.
The late-time behavior is organized by two model quantities. The symbol ϕ represents porosity, while Y controls the reaction contribution. In the physical interpretation used by the analysis, positive Y represents water consumption and negative Y represents water release. The mathematical treatment allows any real-valued Y as long as porosity remains positive, so its regime map is broader than the positive-consumption hydration case.
When Y is greater than −ϕ, the front eventually returns to square-root-in-time behavior, but the coefficient now depends on both porosity and water consumption. The late-time law therefore has the same broad shape as the early-time law while carrying information about the reaction through its prefactor.
At the boundary Y = −ϕ, the late-time front instead moves linearly with time. Below that boundary, when Y < −ϕ, the later motion is exponential. These are conditional mathematical results from the model. Because negative Y is interpreted as water release, the latter regimes belong to the analysis’s wider mathematical parameter range rather than to the positive-Y interpretation of hydration consumption.
What the formulas can—and cannot—say
For physical hydration cases with positive Y, the authors report that the analytical solution validates a numerical solution presented earlier. That links the new formulas to the numerical behavior the paper discusses, including the transition between early- and late-time regimes. The supplied analysis gives no quantitative discrepancy measure, so the reported consistency is not accompanied by a numerical error estimate.
The assumptions define the limits of the result. The domain is one-dimensional and horizontal, the porous medium is homogeneous, the pressure conditions are fixed, the wet–dry boundary is sharp, and the reaction is first-order. Those choices make a closed treatment possible, but they also mean the formulas describe this particular idealized configuration rather than every porous material or infiltration geometry.
This is a deterministic modeling analysis, not a study of participants, physical specimens or an empirical dataset. It reports no sample size, parameter calibration or statistical uncertainty. The uniqueness result is therefore conditional on the equations and their assumptions; it does not establish that the same front behavior has been measured in natural subsurface materials.
The authors also propose using the closed solution as a basis for studying reaction-fracture mechanics linked to hydration-swelling eigenstrains in porous media. That is a possible downstream application, not a demonstrated fracture result. The supplied analysis does not present a fracture-mechanics calculation, quantify fracture risk or validate the proposed connection experimentally.
Several questions consequently remain open: whether the solution matches measured infiltration and hydration in heterogeneous rocks, whether it can support quantitative predictions of reaction-driven fracture, and whether the same regimes survive in multidimensional settings with variable properties or non-first-order kinetics. Independent comparison with the earlier numerical studies would also be needed to assess the reported analytical–numerical consistency.
Publication status
The record identifies the work as an arXiv version 1 preprint dated 20 Aug 2026. Its conclusions should therefore be read as the results of the stated theoretical analysis and its idealized infiltration–hydration model.
Paper data and sources
Original title: On the closed solution of a problem coupling fluid infiltration with a hydration reaction
Authors: Diego Guevara, Sabrina Roscani, Piotr Rybka, Vaughan Voller
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text