arxiv: 2605.12329 · v1 · submitted 2026-05-12 · ⚛️ physics.geo-ph

Recognition: 2 theorem links

· Lean Theorem

Calibration of stress-jump conditions for arbitrary flow directions in fluid-porous systems

Joscha Nickl (1) ((1) I2M), Philippe Angot (1)

Pith reviewed 2026-05-13 02:39 UTC · model grok-4.3

classification ⚛️ physics.geo-ph

keywords stress-jump conditionsStokes-Darcy flowporous mediainterface conditionsfriction tensorparameter calibrationnumerical validationfluid-porous systems

0 comments

The pith

The stress-jump conditions for arbitrary flow directions at fluid-porous interfaces reduce to a single-parameter calibration by exploiting structural properties of the porous medium.

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

The paper validates the stress-jump coupling conditions for Stokes-Darcy flow in two dimensions by calibrating the unknown coefficients of the friction tensor against reference data from processed pore-scale simulations. It shows that structural properties of the porous medium permit reducing the formally three-parameter problem to an effective one-dimensional calibration. This reduction incurs negligible loss in accuracy across tested configurations and flow regimes. A regional sensitivity analysis further establishes that even coarse parameter estimates produce reliable model performance.

Core claim

The stress-jump conditions, formulated for arbitrary flow directions at the interface between a porous medium and an adjacent free-flow region, involve a friction tensor whose coefficients are calibrated by matching the macroscopic model to processed pore-scale simulations. Exploiting structural properties of the porous medium enables an effective reduction to a one-dimensional calibration with negligible loss in accuracy. Regional sensitivity analysis indicates that even coarse parameter estimates can yield a well-performing model.

What carries the argument

The friction tensor within the stress-jump coupling conditions, whose three coefficients are determined by matching macroscopic predictions to pore-scale reference solutions.

Load-bearing premise

The processed pore-scale simulations supply sufficiently accurate reference solutions to calibrate the macroscopic stress-jump model without introducing systematic bias from discretization or processing choices.

What would settle it

Comparison of the one-parameter calibrated model's predicted interface stresses and velocity fields against independent pore-scale simulations for a porous geometry or flow direction excluded from the original calibration set.

read the original abstract

A numerical validation of the stress-jump coupling conditions for Stokes-Darcy flow in two dimensions is presented, addressing a gap that has remained since their introduction by Angot et al.. These conditions, formulated for arbitrary flow directions at the interface between a porous medium and an adjacent free-flow region, involve a friction tensor whose coefficients are not known a priori. We calibrate these parameters for a range of porous-medium configurations and flow regimes by matching the macroscopic model to reference solutions derived from processed pore-scale simulations. Several optimization strategies are assessed for this calibration task. The results show that, although three parameters are formally required, exploiting structural properties of the porous medium enables an effective reduction to a one-dimensional calibration with negligible loss in accuracy. A regional sensitivity analysis further indicates that even coarse parameter estimates can yield a well-performing model, highlighting the robustness and practical applicability of the stress-jump formulation.

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.

Desk Editor's Note private letter to a colleague

This paper calibrates the friction tensor in stress-jump conditions for 2D arbitrary flows and shows that medium structure lets you drop to one effective parameter with little accuracy loss.

read the letter

The core result is a set of calibrated values for the stress-jump friction tensor across several porous configurations and flow angles, plus the finding that the three-parameter tensor can be reduced to a single scalar fit without much penalty. They do this by matching the macroscopic Stokes-Darcy model to pore-scale reference data and by testing a few optimization routes. The regional sensitivity check is useful: it indicates that even coarse guesses for the remaining parameter still give acceptable interface behavior in many cases. That part is practical and directly addresses why the conditions have been hard to use when flow is not aligned with the medium axes. The work is an extension of the earlier Angot formulation rather than a new derivation, but the explicit numbers and the 1D-reduction evidence are new for arbitrary directions. The numerical evidence for retained accuracy after reduction is the strongest part of the contribution. The main soft spot is the dependence on processed pore-scale simulations as the target. Without reported grid-convergence checks at the interface, details on how the interface location was extracted, or quantitative error norms between the calibrated model and the references, it is difficult to judge whether the fitted coefficients reflect the physics or the discretization choices. The abstract does not give those metrics, so the claim of negligible loss rests on data whose quality is not fully visible here. This is a targeted methods paper for people who already model fluid-porous interfaces in geophysics or engineering and want concrete guidance on setting the jump parameters in 2D codes. A reader who needs to implement these conditions would find the calibration examples and the sensitivity results worth looking at. It is solid enough on its own terms to deserve a serious referee, mainly to check the reference-data processing and to see whether the 1D reduction holds under tighter error controls. I would send it to review rather than desk-reject.

Referee Report

1 major / 2 minor

Summary. The manuscript presents a numerical validation of the stress-jump coupling conditions for Stokes-Darcy flow in two dimensions, addressing a gap since their introduction by Angot et al. The conditions involve a friction tensor with coefficients not known a priori. The authors calibrate these parameters for a range of porous-medium configurations and flow regimes by matching the macroscopic model to reference solutions from processed pore-scale simulations, assess several optimization strategies, and conclude that structural properties of the porous medium enable an effective reduction to a one-dimensional calibration with negligible loss in accuracy. A regional sensitivity analysis indicates that even coarse parameter estimates yield a well-performing model.

Significance. If the pore-scale reference solutions are free of systematic bias, the work would be significant for fluid-porous media modeling by providing practical calibration guidance and demonstrating that the stress-jump formulation is robust and can be simplified for arbitrary flow directions. The assessment of multiple optimization strategies and the regional sensitivity analysis are useful strengths that support applicability. The paper supplies concrete numerical evidence for the dimensionality reduction claim rather than purely theoretical arguments.

major comments (1)

[Section describing pore-scale simulation processing and reference extraction] The central claim—that structural properties permit reduction from three parameters to a one-dimensional calibration with negligible accuracy loss, and that the model is robust per the sensitivity analysis—rests on the processed pore-scale simulations supplying unbiased reference data. The manuscript must provide explicit details on interface extraction, averaging procedures, grid resolution at the fluid-porous interface, and any convergence checks to confirm that discretization or processing artifacts are not absorbed into the calibrated friction tensor coefficients. Without this, the reported negligible loss and robustness conclusions cannot be fully evaluated.

minor comments (2)

[Abstract] The abstract states that 'several optimization strategies are assessed' but does not name them; a brief enumeration would improve clarity for readers.
[Introduction or model formulation section] Notation for the friction tensor and its coefficients should be introduced with a clear equation reference early in the text to aid readers unfamiliar with the Angot et al. formulation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive and detailed review. We agree that the central claims regarding the dimensionality reduction and model robustness depend on the reliability of the processed pore-scale reference data, and that the manuscript requires additional explicit details on the processing procedures to allow full evaluation of potential biases or artifacts.

read point-by-point responses

Referee: [Section describing pore-scale simulation processing and reference extraction] The central claim—that structural properties permit reduction from three parameters to a one-dimensional calibration with negligible accuracy loss, and that the model is robust per the sensitivity analysis—rests on the processed pore-scale simulations supplying unbiased reference data. The manuscript must provide explicit details on interface extraction, averaging procedures, grid resolution at the fluid-porous interface, and any convergence checks to confirm that discretization or processing artifacts are not absorbed into the calibrated friction tensor coefficients. Without this, the reported negligible loss and robustness conclusions cannot be fully evaluated.

Authors: We acknowledge that the current description of the pore-scale simulation processing and reference extraction is at a summary level and does not include the requested specifics. In the revised manuscript, we will add a dedicated subsection that explicitly details: the interface extraction algorithm (including how the sharp interface is identified from the pore geometry), the spatial and temporal averaging procedures used to compute the macroscopic velocity and stress fields for comparison, the grid resolution at the fluid-porous interface (with quantitative values such as minimum cell size relative to pore diameter), and the convergence checks performed by successive grid refinement to quantify discretization error in the extracted reference data. These additions will confirm that processing artifacts are not inadvertently absorbed into the calibrated friction tensor and will thereby strengthen the support for the one-dimensional calibration claim and the regional sensitivity analysis. revision: yes

Circularity Check

0 steps flagged

Minor self-citation to introduction of conditions; calibration independent of inputs

full rationale

The paper's core results consist of fitting the friction tensor coefficients in the stress-jump conditions to reference data obtained from separate pore-scale simulations, followed by empirical demonstration that structural properties of the medium permit reduction to a one-parameter calibration with little accuracy loss. This process is not self-referential: the target quantities are external simulation outputs, not quantities defined from the macroscopic model itself. The single self-citation to Angot et al. supplies only the starting form of the interface conditions (whose coefficients are stated to be unknown a priori) and does not carry the load of the calibration outcomes, optimization comparisons, or regional sensitivity conclusions. No equation or claim reduces to its own inputs by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the validity of the Stokes-Darcy macroscopic model and the accuracy of pore-scale references used for fitting; the three friction-tensor coefficients are the fitted quantities.

free parameters (1)

friction tensor coefficients
Three coefficients calibrated by optimization against pore-scale data; structural symmetry reduces effective number to one.

axioms (1)

domain assumption Stokes-Darcy equations and interface jump conditions hold for the chosen flow regimes
Invoked throughout as the macroscopic model being calibrated.

pith-pipeline@v0.9.0 · 5457 in / 1211 out tokens · 49038 ms · 2026-05-13T02:39:18.381915+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
although three parameters are formally required, exploiting structural properties of the porous medium enables an effective reduction to a one-dimensional calibration with negligible loss in accuracy
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
The friction tensor β is an unknown quantity. Thus, the stress-jump model needs to be calibrated.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages

[1]

Angot, P., Goyeau, B., Ochoa-Tapia, J.A.: Asymptotic modeling of transport phenom- ena at the interface between a fluid and a porous layer: Jump conditions. Phys. Rev. E 95, 063302 (2017) https://doi.org/10.1103/ PhysRevE.95.063302

work page 2017
[2]

Das, D.B.: Hydrodynamic modelling for groundwater flow through permeable reactive barriers. Hydrol. Process.16(17), 3393–3418 (2002) https://doi.org/10.1002/hyp.1107

work page doi:10.1002/hyp.1107 2002
[3]

ESAIM: M2AN 45(3), 447–476 (2011) https://doi

D’Angelo, C., Zunino, P.: Robust numeri- cal approximation of coupled Stokes’ and Darcy’s flows applied to vascular hemody- namics and biochemical transport. ESAIM: M2AN 45(3), 447–476 (2011) https://doi. org/10.1051/m2an/2010062

work page doi:10.1051/m2an/2010062 2011
[4]

Finite Elem

Dereims, A., Drapier, S., Bergheau, J.-M., De Luca, P.: 3D robust iterative coupling of Stokes, Darcy and solid mechanics for low permeability media undergoing finite strains. Finite Elem. Anal. Des. 94, 1–15 (2015) https://doi.org/10.1016/j.finel.2014.09.003

work page doi:10.1016/j.finel.2014.09.003 2015
[5]

Narsilio, G.A., Buzzi, O., Fityus, S., Yun, T.S., Smith, D.W.: Upscaling of Navier–Stokes equations in porous media: Theoretical, numerical and experimen- tal approach. Comput. Geotech. 36(7), 1200–1206 (2009) https://doi.org/10.1016/j. compgeo.2009.05.006

work page doi:10.1016/j 2009
[6]

Fredrich, J.T., DiGiovanni, A.A., Noble, D.R.: Predicting macroscopic transport prop- erties using microscopic image data. J. Geo- phys. Res. 111(B3), 2005–003774 (2006) https://doi.org/10.1029/2005JB003774

work page doi:10.1029/2005jb003774 2005
[7]

Springer, New York (1997)

Hornung, U.: Homogenization and Porous Media. Springer, New York (1997). https:// doi.org/10.1007/978-1-4612-1920-0

work page doi:10.1007/978-1-4612-1920-0 1997
[8]

Whitaker, S.: Flow in porous media I: A the- oretical derivation of Darcy’s law. Transp. Porous Media 1(1), 3–25 (1986) https://doi. org/10.1007/BF01036523

work page doi:10.1007/bf01036523 1986
[9]

Elhakim, A.F.: Estimation of soil permeabil- ity. Alex. Eng. J. 55(3), 2631–2638 (2016) https://doi.org/10.1016/j.aej.2016.07.034

work page doi:10.1016/j.aej.2016.07.034 2016
[10]

Vervoort, R.W., Cattle, S.R.: Linking hydraulic conductivity and tortuosity parameters to pore space geometry and pore-size distribution. J. Hydrol. 272(1- 4), 36–49 (2003) https://doi.org/10.1016/ S0022-1694(02)00253-6

work page 2003
[11]

Valdés-Parada, F.J., Lasseux, D.: Flow near porous media boundaries including inertia and slip: A one-domain approach. Phys. Flu- ids 33(7), 073612 (2021) https://doi.org/10. 1063/5.0056345

work page 2021
[12]

Hernandez-Rodriguez,R.,Angot,P.,Goyeau, B., Ochoa-Tapia, J.A.: Momentum transport in the free fluid-porous medium transition layer: One-domain approach. Chem. Eng. Sci. 248, 117111 (2022) https://doi.org/10.1016/ j.ces.2021.117111

work page arXiv 2022
[13]

Lyu, W., Wang, X.: Stokes–Darcy system, small-Darcy-number behaviour and related 27 interfacial conditions. J. Fluid Mech.922, 4 https://doi.org/10.1017/jfm.2021.509

work page doi:10.1017/jfm.2021.509 2021
[14]

Carraro, T., Goll, C., Marciniak-Czochra, A., Mikelić, A.: Pressure jump interface law for the Stokes-Darcy coupling: Confirmation by direct numerical simulations. J. Fluid Mech. 732, 510–536 (2013) https://doi.org/ 10.1017/jfm.2013.416 1301.6580

work page doi:10.1017/jfm.2013.416 2013
[15]

Beavers, G.S., Joseph, D.D.: Boundary condi- tions at a naturally permeable wall. J. Fluid Mech. 30(1),197–207(1967)https://doi.org/ 10.1017/S0022112067001375

work page doi:10.1017/s0022112067001375 1967
[16]

Jones, I.P.: Low Reynolds number flow past a porous spherical shell. Math. Proc. Camb. Phil. Soc. 73(1), 231–238 (1973) https://doi. org/10.1017/S0305004100047642

work page doi:10.1017/s0305004100047642 1973
[17]

Nield, D.A.: The Beavers–Joseph boundary condition and related matters: A histori- cal and critical note. Transp. Porous Media 78(3), 537–540 (2009) https://doi.org/10. 1007/s11242-009-9344-y

work page 2009
[18]

Yang, G., Coltman, E., Weishaupt, K., Terzis, A., Helmig, R., Weigand, B.: On the Beavers–Joseph interface condition for non- parallel coupled channel flow over a porous structure at high Reynolds numbers. Transp. Porous Media 128(2), 431–457 (2019) https: //doi.org/10.1007/s11242-019-01255-5

work page doi:10.1007/s11242-019-01255-5 2019
[19]

Strohbeck, P., Eggenweiler, E., Rybak, I.: A modification of the Beavers–Joseph condition for arbitrary flows to the fluid–porous interface. Transp. Porous Media 147(3), 605–628 (2023) https: //doi.org/10.1007/s11242-023-01919-3

work page doi:10.1007/s11242-023-01919-3 2023
[20]

Mierzwiczak, M., Fraska, A., Grabski, J.K.: Determination of the slip constant in the Beavers-Joseph experiment for lami- nar fluid flow through porous media using a meshless method. Math. Probl. Eng. 2019(1), 1494215 (2019) https://doi.org/10. 1155/2019/1494215

work page 2019
[21]

Eggenweiler, E., Rybak, I.: Unsuitability of the Beavers–Joseph interface condition for filtration problems. J. Fluid Mech. 892, 10 (2020) https://doi.org/10.1017/jfm.2020.194

work page doi:10.1017/jfm.2020.194 2020
[22]

Theo- reticaldevelopment.Int.J.HeatMassTransf

Ochoa-Tapia, J.A., Whitaker, S.: Momentum transfer at the boundary between a porous medium and a homogeneous fluid—I. Theo- reticaldevelopment.Int.J.HeatMassTransf. 38(14), 2635–2646 (1995) https://doi.org/10. 1016/0017-9310(94)00346-W

work page 1995
[23]

International Journal of Heat and Mass Transfer , volume =

Ochoa-Tapia, J.A., Whitaker, S.: Momentum transfer at the boundary between a porous medium and a homogeneous fluid—II. Com- parison with experiment. Int. J. Heat Mass Transf. 38(14), 2647–2655 (1995) https:// doi.org/10.1016/0017-9310(94)00347-X

work page doi:10.1016/0017-9310(94)00347-x 1995
[24]

Valdés-Parada, F.J., Aguilar-Madera, C.G., Ochoa-Tapia, J.A., Goyeau, B.: Velocity and stress jump conditions between a porous medium and a fluid. Adv. Water Resour. 62, 327–339 (2013) https://doi.org/10.1016/ j.advwatres.2013.08.008

work page 2013
[25]

https://hal.science/hal-04711298/ (in revi- sion) (2021)

Angot, P., Goyeau, B., Ochoa-Tapia, J.A.: An optimal stress jump interface condition for the fluid-porous multi-dimensional flow. https://hal.science/hal-04711298/ (in revi- sion) (2021)

work page 2021
[26]

working paper or preprint (2021)

Angot, P.: Solvability of the variable-viscosity fluid-porous flows coupled with an opti- mal stress jump interface condition. working paper or preprint (2021). https://hal.science/ hal-03172378

work page 2021
[27]

In: Proceedings of the 39th IAHR World Congress, pp

Hassan, A., Tügel, F., Özgen, I., Molkenthin, F., Hinkelmann, R.: Automatic calibration of a shallow water model for rainfall-runoff simulations using fast global optimization. In: Proceedings of the 39th IAHR World Congress, pp. 4647–4655 (2022). https://doi. org/10.3850/IAHR-39WC25217119202297

work page doi:10.3850/iahr-39wc25217119202297 2022
[28]

Com- put

Matott, L.S., Leung, K., Sim, J.: Applica- tion of MATLAB and Python optimizers to two case studies involving groundwater flow and contaminant transport modeling. Com- put. Geosci. 37(11), 1894–1899 (2011) https: //doi.org/10.1016/j.cageo.2011.03.017

work page doi:10.1016/j.cageo.2011.03.017 2011
[29]

Com- mun

Cao, Y., Gunzburger, M., Hua, F., Wang, X.: 28 Coupled Stokes–Darcy model with Beavers– Joseph interface boundary condition. Com- mun. Math. Sci. 8, 1–25 (2010) https://doi. org/10.4310/CMS.2010.v8.n1.a2

work page doi:10.4310/cms.2010.v8.n1.a2 2010
[30]

Angot, P.: On the well-posed coupling between free fluid and porous viscous flows. Appl. Math. Lett. 24(6), 803–810 (2011) https://doi.org/10.1016/j.aml.2010.07.008

work page doi:10.1016/j.aml.2010.07.008 2011
[31]

Hecht, F.: New development in freefem++. J. Numer. Math.20(3-4), 251–265 (2012) https: //doi.org/10.1515/jnum-2012-0013

work page doi:10.1515/jnum-2012-0013 2012
[32]

Correa, M.R., Loula, A.F.D.: Uncondition- ally stable mixed finite element methods for Darcy flow. Comput. Methods Appl. Mech. Engrg. 197(17-18), 1525–1540 (2008) https: //doi.org/10.1016/j.cma.2007.11.025

work page doi:10.1016/j.cma.2007.11.025 2008
[33]

Correa, M.R., Loula, A.F.D.: A unified mixed formulation naturally coupling Stokes and Darcy flows. Comput. Methods Appl. Mech. Engrg. 198(33-36), 2710–2722 (2009) https: //doi.org/10.1016/j.cma.2009.03.016

work page doi:10.1016/j.cma.2009.03.016 2009
[34]

PhD thesis, Thèse de doctorat de l’Université Paris-Saclay pré- parée à l’Ecole polytechnique (2019)

Feppon, F.: Shape and topology optimization of multiphysics systems. PhD thesis, Thèse de doctorat de l’Université Paris-Saclay pré- parée à l’Ecole polytechnique (2019)

work page 2019
[35]

arXiv (2021)

Rinehart, A., L¯ acis, U., Bagheri, S.: The Brinkman viscosity for porous media exposed to a free flow. arXiv (2021). https://doi.org/ 10.48550/ARXIV.2106.01879 . https://arxiv. org/abs/2106.01879

work page doi:10.48550/arxiv.2106.01879 2021
[36]

Xiang, Y., Sun, D.Y., Fan, W., Gong, X.G.: Generalized simulated annealing algorithm and its application to the Thomson model. Phys. Lett. A233(3), 216–220 (1997) https: //doi.org/10.1016/S0375-9601(97)00474-X

work page doi:10.1016/s0375-9601(97)00474-x 1997
[37]

Soritz, S., Moser, D., Gruber-Wölfler, H.: Comparison of derivative-free algorithms for their applicability in self-optimization of chemical processes. Chem. Methods 2(5), 202100091 (2022) https://doi.org/10.1002/ cmtd.202100091

work page 2022
[38]

Virtanen, P., Gommers, R., Oliphant, T.E., Haberland, M., Reddy, T., Cournapeau, D., Burovski, E., Peterson, P., Weckesser, W., Bright, J., van der Walt, S.J., Brett, M., Wilson, J., Millman, K.J., Mayorov, N., Nel- son, A.R.J., Jones, E., Kern, R., Larson, E., Carey, C.J., Polat, I., Feng, Y., Moore, E.W., VanderPlas, J., Laxalde, D., Perk- told, J., Cim...

work page doi:10.1038/s41592-019-0686-2 2020
[39]

Pianosi, F., Beven, K., Freer, J., Hall, J.W., Rougier, J., Stephenson, D.B., Wagener, T.: Sensitivity analysis of environmental mod- els: A systematic review with practical workflow. Environ. Model. Softw. 79, 214– 232 (2016) https://doi.org/10.1016/j.envsoft. 2016.02.008

work page doi:10.1016/j.envsoft 2016
[40]

Wagener, T., Boyle, D.P., Lees, M.J., Wheater, H.S., Gupta, H.V., Sorooshian, S.: A framework for development and applica- tion of hydrological models. Hydrol. Earth Syst. Sci.5(1), 13–26 (2001) https://doi.org/ 10.5194/hess-5-13-2001

work page doi:10.5194/hess-5-13-2001 2001
[41]

Socio-Env

Iwanaga, T., Usher, W., Herman, J.: Toward SALib 2.0: Advancing the accessibility and interpretability of global sensitivity analy- ses. Socio-Env. Syst. Model.4, 18155 (2022) https://doi.org/10.18174/sesmo.18155

work page doi:10.18174/sesmo.18155 2022
[42]

Herman, J., Usher, W.: SALib: An open- source python library for sensitivity analysis. J. Open Source Softw. 2(9) (2017) https:// doi.org/10.21105/joss.00097

work page doi:10.21105/joss.00097 2017
[43]

Valdés-Parada, F.J., Lasseux, D.: A novel one-domain approach for modeling flow in a fluid-porous system including inertia and slip effects. Phys. Fluids33(2), 022106 (2021) https://doi.org/10.1063/5.0036812

work page doi:10.1063/5.0036812 2021
[44]

Sahraoui, M., Kaviany, M.: Slip and no-slip velocity boundary conditions at interface of porous, plain media 35(4), 927–943 (1991) https://doi.org/10.1016/0017-9310(92) 90258-T 29 Appendix A Experimental convergence analysis We implemented a Stokes–Darcy problem with right hand sides in the flow conditions and in the interface conditions. Because of these...

work page doi:10.1016/0017-9310(92 1991