{"id":"f4550ae0-b6eb-43ab-a4ea-be9e56162aff","arxiv_id":"1908.01771","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A quasi-3D Stokes free-flow model coupled to a pore-network model with a locally fitted slip coefficient reproduces integrated flow rates of a full 3D micromodel simulation within about 6% at a fraction of the runtime.","lead":"This paper speeds up simulations of fluid flowing over a porous material by replacing the intricate pore structure with a simplified pipe network, coupled to a fast quasi-3D flow model. The method cuts runtime from hours to minutes while keeping integrated flow rates within about 10% of full simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (3.16), (3.18), and (4.4) are dimensionally inconsistent: β is reported in m^-1 but appears in a stress balance requiring Pa·s/m, so the fitted slip coefficient's convention is unspecified.","rationale":"The paper has several independently supported elements: the 3D OpenFOAM model is checked against micro-PIV data, the quasi-3D model is compared to the 3D solution, and integrated throat fluxes in Fig. 20 are compared across models. These give real support to the model-reduction chain. However, the novel contribution is the slip closure, and that closure has an internal dimensional inconsistency. The reader focused on transferability of the fitted β and on the calibration/validation loop; those are legitimate, and I agree they weaken the claim. But the more direct problem is that the equations defining β cannot be simultaneously true. Eq. (3.16) and Eq. (3.18) require β in Pa·s/m, while Eq. (4.4) and the reported values give β in m^-1. Without knowing which convention the code uses, the fitted coefficient and all error metrics that rely on it are not reproducible. This is a correctness risk in the central model formulation, not a disagreement with the community consensus. It does not necessarily invalidate the numerical results — the code may use a consistent convention that is simply misreported — but it makes the paper's core equations ambiguous. Therefore the appropriate outcome remains CONDITIONAL: the authors must clarify or correct the units and confirm the implementation, in addition to the reader's request for independent validation. I keep the reader's verdict unchanged but for a more fundamental reason.","tokens_in":19919,"tokens_out":11855,"duration_ms":131013,"concrete_test":"Perform a dimensional audit of the DuMux implementation and recompute the Table 1 β values including the omitted µ. Concretely: inspect the source code of the slip coupling condition (Eq. 3.18) to determine whether it evaluates (1/β)*µ*shear_rate (β in m^-1) or (1/β)*shear_rate (β in Pa·s/m). Then re-derive Eq. (4.4) with the correct factor: if the intended coefficient is the stress coefficient in Eq. (3.16), β should be µ⟨∂v_x/∂y + ∂v_y/∂x⟩/⟨v_x⟩. Recompute β for w_throat = 100 µm; if the value changes from 84,550 m^-1 to 84.55 Pa·s/m (or vice versa), rerun the coupled micromodel with the dimensionally consistent coefficient. If the reported error metrics change by more than a few percent, the central claim depends on the unspecified convention; if they are unchanged, the inconsistency is a reporting error and the verdict is unaffected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.3 introduces the slip closure in Eq. (3.16) as a balance of tangential stress: [(-µ(∇v+∇v^T)·n)·t]_FF = β_throat ([v]_FF - [v]_PNM)·t. The left side has units Pa (N/m^2). In Eq. (3.18), v_slip = (1/β_throat)[(-µ(∇v+∇v^T)·n)·t]_FF + [v]_PNM·t, so β_throat must have units Pa·s/m for v_slip to be a velocity. But the calibration formula in Eq. (4.4), β ≈ ⟨∂v_x/∂y + ∂v_y/∂x⟩/⟨v_x⟩, yields units 1/m, and all tabulated β values (e.g., 33,000 m^-1) are 1/m. Thus Eqs. (3.16)/(3.18) and Eq. (4.4) cannot both be correct: either the left side of Eq. (3.16) should be the shear rate (drop µ), or β_throat as used in the code is µ times the fitted quantity and has units Pa·s/m, in which case the reported values are off by the viscosity factor. The paper does not state which convention the DuMux implementation uses, so the central quantitative claim — relative velocity error 9.49e-3, pressure error 3.72e-4, '<10% deviation' — rests on an unspecified normalization. This is more fundamental than the transferability question: until the coefficient is unambiguously defined, the fitted β cannot be checked, reproduced, or even assigned a physical meaning.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a two-stage model reduction for steady creeping flow in a micromodel consisting of a free-flow channel over a regular porous medium. The authors first build a fully resolved 3D OpenFOAM Stokes model, validate it against high-resolution micro-PIV data, and then reduce it to a quasi-3D 2D model with a wall-friction closure. In the second stage, the porous part is replaced by a pore-network model, coupled monolithically to the quasi-3D free-flow model. The new element is a slip boundary condition at each pore-throat intersection with the free-flow interface, parameterized by a per-throat slip coefficient β_throat that is calibrated from isolated quasi-3D reference solutions. The reported results include a global relative velocity error of 9.49e-3 and a pressure error of 3.72e-4 for the coupled model against the quasi-3D solution, integral throat fluxes within about 6% of the 3D simulation, and a runtime reduction from roughly five hours on 30 cores to under five minutes on a single core.","tokens_in":20344,"tokens_out":7056,"duration_ms":71129,"significance":"If the quantitative claims can be made fully reproducible, this is a useful contribution: the per-throat slip closure removes the artificial no-slip constraint of the earlier model, improves the velocity error by a factor of about 9 in the single-throat test, and retains pore-scale detail at a small fraction of the cost of the resolved simulation. The fact that the 3D reference model is compared with micro-PIV experiments gives the reduction chain a credible anchor. The paper is also reasonably candid about limitations, including the degraded behavior for inclined throats and the large local deviations near the edges of throat openings in Fig. 19. However, the central error claims currently rest on a slip coefficient whose definition is dimensionally inconsistent between the model equations and the calibration formula, and the coupled model is validated against the same quasi-3D model class used for calibration. These issues must be resolved before the paper can be considered ready.","major_comments":[{"comment":"The slip coefficient is dimensionally inconsistent. In Eq. (3.16), the left-hand side [(-µ(∇v+∇v^T)·n)·t]_FF is a shear stress with units Pa, so β_throat in Eq. (3.18) must have units Pa·s/m for v_slip to be a velocity. However, Eq. (4.4) defines β_throat as a shear rate divided by a velocity, which has units 1/m, and Table 1 lists β in 1/m. Unless the DuMux implementation multiplies the fitted quantity by the viscosity (which is nowhere stated), the tabulated values cannot be substituted into Eq. (3.18). The central quantitative claims — relative velocity error 9.49e-3, pressure error 3.72e-4, and the abstract's '<10%' statement — all depend on this parameter. The authors must state the convention used in the implementation and make the reported numbers consistent with that convention before the results can be reproduced or checked.","section":"Sec. 3.3, Eqs. (3.16)/(3.18); Sec. 4.3, Eq. (4.4); Table 1"},{"comment":"The abstract's claim that the coupled model deviates by less than 10% from the other two model concepts is not supported as stated. The global relative L2 errors reported in Sec. 4.4 are 9.49e-3 for velocity and 3.72e-4 for pressure, but Fig. 19 shows local vertical-velocity discrepancies of up to 80% at the edges of the throat openings. The paper should explicitly state that the 10% figure refers to global norms, quantify the local deviations in the abstract, and discuss their practical consequence for quantities such as interfacial mass exchange, as is partially done around Fig. 20.","section":"Abstract; Sec. 4.4, Fig. 19"},{"comment":"The validation loop is partially circular. The throat conductances kij are obtained by numerical upscaling on a quasi-3D pillar structure, and β_throat is fitted from isolated quasi-3D reference solutions; the full coupled model is then compared against the quasi-3D solution of the same micromodel. The reported velocity and pressure errors are therefore consistency checks of the reduction chain, not independent tests against the 3D reference or the experimental data. Direct comparison with the 3D OpenFOAM fields is provided only for integrated throat fluxes in Fig. 20. The paper should either add a direct comparison of the coupled model's velocity/pressure fields with the 3D solution or explicitly restrict the claims to quasi-3D consistency.","section":"Secs. 4.3 and 4.4"},{"comment":"The transferability of a single β_throat value to the full 81-throat network is fragile and is not quantified. The calibration set gives β_throat varying by roughly a factor of 9 over throat widths from 50 to 400 µm, and the inclined-throat test changes β from 84550 to 28001 1/m for the same 100 µm width. The full coupled model uses a single value β = 33000 1/m from the power-law fit. The text states that the re-evaluated β gives 'very similar results,' but no results are shown to support this. A sensitivity study of the coupled model's global errors with respect to β_throat is needed to justify the single-parameter closure and to determine how much of the 80% local deviation in Fig. 19 is controlled by this parameter.","section":"Sec. 4.3, Figs. 12 and 14"}],"minor_comments":[{"comment":"The term 'βthoat' is a typo and should read 'βthroat'.","section":"Sec. 4, first paragraph"},{"comment":"The text near Eq. (3.9) contains 'Qis' and an incomplete definition; it should read 'Qij' with a complete sentence defining the discrete flow rate.","section":"Sec. 3.3, Eq. (3.9)"},{"comment":"The terminology switches between 'conductance' and 'resistance' for kij; the authors should use one term consistently and state the units in the text.","section":"Sec. 4.3"},{"comment":"The power-law fits should state the units of the prefactor and of wthroat, and should be reconciled with the stated unit of β; currently the plotted values inherit the dimensional ambiguity of Eq. (4.4).","section":"Fig. 12"},{"comment":"The CPU-time comparison is appropriately caveated in Sec. 5, but the abstract should mention that the 3D solver is transient and parallelized while the reduced models are stationary and serial, so the speedup is not a direct wall-clock comparison.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"The dimensional inconsistency of β_throat is likely fixable by clarifying the implementation convention, but it is load-bearing because every quantitative accuracy claim depends on it. The more substantive issue is the calibration-to-same-model validation loop; the authors should be pushed to add at least one direct comparison of the coupled model with the 3D OpenFOAM fields or with the micro-PIV data along the interface. The core idea is sound and within the scope of JCP, so I would be willing to review a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the local, per-throat slip condition at pore-throat/free-flow intersections. It fixes the no-slip artifact of the authors' earlier model for orthogonal throats. The implementation is clean, the slip term clearly improves the velocity field (a factor of 9 reduction in the isolated-throat velocity error), and the runtime speedup is credible.\n\nThe soft spots, in order of severity. First, there is a dimensional problem in the definition of beta_throat. Eq (3.16)/(3.18) require beta to have units Pa s/m, but Eq (4.4) and all tabulated values give 1/m. Either the stress balance should use the shear rate without the viscosity, or the fitted beta is intended to be multiplied by mu. The paper does not say which convention the DuMux implementation uses, so the central quantitative claims are not reproducible until this is resolved. I cannot assign a physical meaning to the fitted value as written. This is the most fundamental issue.\n\nSecond, the abstract's 'less than 10%' is too broad. Fig 19 shows local vertical-velocity discrepancies up to 80% at the edges of throat openings. The global errors are small and the mass exchange is captured, but the claim needs to be qualified to integrated or global quantities.\n\nThird, the calibration/validation loop is partly circular. beta_throat and the conductances kij are fitted to quasi-3D reference solutions, and then the coupled model is validated against the same quasi-3D model class. The comparison with 3D OpenFOAM throat fluxes in Fig 20 helps, but the coupled velocity field is never checked against the full 3D simulation or the experiment. The transferability of beta from an isolated orthogonal throat to the 81-throat network is fragile—the inclined-throat test shows beta changes by a factor of 3 for the same width—though in this geometry all interface throats are orthogonal, so it is an acceptable assumption for the main test case.\n\nNone of this kills the paper. The central idea is sound, the numerics are careful, and the authors are honest about limitations. There is real value for people building hybrid-dimensional models of microfluidic and porous-media flows. It deserves peer review, but I would send it back with a request to fix the units, qualify the accuracy claims, and ideally release code and data.","headline":"A useful local slip closure for pore-network/free-flow coupling, but a unit inconsistency in the fitted slip coefficient and an overbroad accuracy claim should be fixed before I'd trust the numbers.","tokens_in":20837,"tokens_out":3798,"would_cite":false,"duration_ms":35845,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two-stage model reduction reproduces a micromodel's pore-scale flow within 6 percent while cutting runtime from hours to minutes.","keywords":["free flow over porous media","pore-network model","hybrid-dimensional model","model reduction","slip velocity","quasi-3D Stokes flow","micro-PIV","micromodel"],"falsifier":"Take the same micromodel with the pore openings tilted by 45 degrees, run the coupled model with the coefficient fitted on perpendicular openings, and compare the velocity immediately above each opening with a fully resolved 3D simulation; if the integrated throat fluxes leave the 6 percent agreement band, the single-coefficient transferability fails.","tokens_in":19718,"feed_emoji":"🌊","tokens_out":15427,"duration_ms":135013,"temperature":0.7,"pith_summary":"Simulating free flow over a porous medium with fully resolved three-dimensional flow equations costs too much for all but small domains, so the authors ask how far the geometry can be reduced while keeping acceptable accuracy. They first collapse the 3D problem to a quasi-3D plane using a wall-friction drag term, then replace the porous pillar array by an equivalent pore network, leaving only the free-flow channel resolved. The new step is a per-throat slip condition at openings between channel and network, with a slip coefficient $\\beta_{\\text{throat}}$ fitted numerically on isolated throat geometries and extrapolated by a power-law fit. Against the fully resolved reference, the coupled model reproduces integrated throat fluxes within about 6%, keeps velocity and pressure relative errors near $10^{-2}$ and $10^{-4}$, and cuts runtime from roughly five hours on 30 cores to under five minutes on one core. The authors' claim is that this simple, local modification makes hybrid-dimensional models accurate enough for larger coupled free-flow/porous problems.","feed_headline":"Two-step model matches 3D pore-flow fluxes within 6%","feed_subtitle":"A fitted slip velocity at each pore throat reproduces full 3D flow in minutes instead of hours.","key_machinery":"The carrying mechanism is a local slip coupling condition at each pore throat that opens onto the free-flow channel. It prescribes the tangential free-flow velocity at the opening as $$v_{\\text{slip},i} = \\frac{1}{\\beta_{\\text{throat}}}\\big[(-\\mu(\\nabla\\mathbf{v}+\\nabla\\mathbf{v}^T)\\mathbf{n})\\cdot\\mathbf{t}_i\\big]^{\\text{FF}} + [\\mathbf{v}]^{\\text{PNM}}\\cdot\\mathbf{t}_i,$$ with $\\beta_{\\text{throat}}$ determined numerically from an isolated orthogonal-throat problem; this replaces the earlier condition, which yielded no slip whenever the throat was perpendicular to the interface. The rest of the reduction chain consists of the quasi-3D Stokes model with drag term $-(8\\mu/h^2)\\mathbf{v}_{2D}$, pore-network throat conductances $k_{ij}$ obtained by numerical upscaling on a reduced equivalent pillar structure, and a fully monolithic (single-system) coupling that enforces mass and pressure continuity at the interface. The empirical power-law fit $\\beta_{\\text{throat}}\\approx 4.531\\,w_{\\text{throat}}^{-1.067}$ (units of $\\mathrm{m}^{-1}$, $R^2=0.999$) extends the calibrated coefficient to other throat widths without re-solving the calibration problem; without wall friction the fit is $\\approx 7.680\\,w_{\\text{throat}}^{-1.019}$.","core_discovery":"The central claim is that a hybrid-dimensional model can stand in for a fully resolved three-dimensional simulation of a micromodel with a free-flow channel over a porous pillar array, as long as the pore-network coupling includes slip at the throat mouths. The earlier coupling forced the tangential component of the pore-throat velocity onto the free-flow boundary, which degenerates to no slip for throats meeting the interface at right angles; this paper replaces that with continuity of tangential stress, balancing the free-flow shear against a per-throat slip coefficient $\\beta_{\\text{throat}}$ fitted on a simplified single-throat geometry. For the full micromodel, the coupled model's integrated fluxes through all 81 interface throats fall within about 6% of the 3D reference, its velocity and pressure fields deviate from the quasi-3D reference by relative errors of $9.5\\times10^{-3}$ and $3.7\\times10^{-4}$, and the runtime is about 5 minutes on one core versus about 5 hours on 30 cores. Including the slip term improves the isolated orthogonal-throat velocity error by a factor of 9 over the no-slip coupling, with no added runtime; the authors also show the improvement persists, but is weaker, for an inclined throat with inflow.","pith_inferences":["Editorial inference: the single-coefficient calibration is the fragile link; the paper's own inclined-throat test shows $\\beta_{\\text{throat}}$ dropping from $84550$ to $28001$ at the same $100\\,\\mu\\mathrm{m}$ width, so applying the orthogonal no-inflow coefficient to a network with inclined or cross-flowing throats is a direct test of where the 6% flux agreement would break down.","Editorial inference: because the slip coefficient is extracted from the same quasi-3D model class used as the validation reference, the comparison does not independently confirm the slip law against true 3D flow; high-resolution velocity measurement directly at a throat mouth would close that gap.","Editorial inference: local errors concentrate at throat edges and cancel in integrated fluxes, so interfacial quantities that respond to local shear—surface deposition, dissolution, or wall reactions—should be treated with caution even where total mass exchange is accurate.","Editorial inference: generalizing $\\beta_{\\text{throat}}$ to depend on throat orientation and inflow strength, or replacing the scalar fit with a local analytical slip expression, is the natural next step toward realistic geometries where interface throats are rarely perfectly perpendicular."],"forward_implications":["The slip condition is local and adds no runtime, so the same coupling can be extended to other free-flow/pore-network geometries of this type; in the isolated orthogonal-throat test it reduces the velocity error by a factor of 9 compared with the no-slip coupling.","Integrated interface quantities are more reliable than local ones: the coupled model over-predicts some near-throat vertical velocity peaks by up to 80%, yet the total flux through each interface throat stays within about 6% of the 3D reference.","The power-law fit $\\beta_{\\text{throat}}\\approx 4.531\\,w_{\\text{throat}}^{-1.067}$ lets the calibrated slip coefficient be reused for other throat widths without re-solving the calibration problem, at least for perpendicular, no-inflow throat configurations.","The reduction shortens a full-array micromodel simulation from about five hours on 30 cores to under five minutes on one core, and the coupled model is 2.3 times faster than the quasi-3D model alone, making parameter sweeps and design optimization practical."],"supporting_citations":[{"why":"Supplies the micromodel geometry and the micro-PIV velocity data against which the fully resolved 3D reference solution is checked.","marker":"[1]"},{"why":"Provides the original monolithic coupling of free flow and pore-network model and the numerical upscaling of throat conductances $k_{ij}$ that the new slip coupling extends.","marker":"[2]"},{"why":"Introduces the wall-friction drag term that converts the 3D Stokes problem into the quasi-3D model used in the first reduction.","marker":"[23]"},{"why":"Precedent for depth-averaged single-phase flow simulation in pseudo-two-dimensional micromodels, supporting the quasi-3D reduction.","marker":"[24]"},{"why":"Extends depth-averaged flow modeling to micromodels with variable aperture, another precedent for the same reduction.","marker":"[25]"},{"why":"Lists analytical conductance expressions for non-circular capillaries; its lack of a suitable expression for this geometry motivates the numerical upscaling of $k_{ij}$.","marker":"[34]"},{"why":"Derives discrete pore-throat flux from one-dimensional Stokes flow, justifying the pressure-continuity coupling condition at the interface.","marker":"[35]"},{"why":"Supplies the classical permeable-wall slip-condition analogy on which the new per-throat slip parametrization is based.","marker":"[36]"},{"why":"Precedent for numerically determining an effective slip parameter from local flow data, though for non-creeping flow.","marker":"[37]"},{"why":"Supplies the analytical channel-flow velocity profile used to verify the 3D mesh resolution and the integrated flux computation.","marker":"[41]"}],"fun_headline_variants":["Slip at pore throats cuts coupled free-flow model error by 9x","Hybrid model with slip matches 3D flow within minutes","Pore-throat slip corrects hybrid model to 6% of full 3D","Quasi-3D with slip runs 60x faster, matches 3D fluxes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result depends on a single fitted number—how strongly the fluid slips at a pore opening—being good enough for every opening in the whole porous structure, even though it was calibrated on one simple opening with no sideways current.","fun_headline_variants_meta":{"raw":{"variants":["Slip at pore throats cuts coupled free-flow model error by 9x","Hybrid model with slip matches 3D flow within minutes","Pore-throat slip corrects hybrid model to 6% of full 3D","Quasi-3D with slip runs 60x faster, matches 3D fluxes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000337,"raw_usage":{"total_tokens":1910,"prompt_tokens":1038,"completion_tokens":872,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":785}},"tokens_in":654,"tokens_out":872,"duration_ms":7693,"temperature":1.0,"reasoning_tokens":785,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:04:07.644603+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same micromodel with the pore openings tilted by 45 degrees, run the coupled model with the coefficient fitted on perpendicular openings, and compare the velocity immediately above each opening with a fully resolved 3D simulation; if the integrated throat fluxes leave the 6 percent agreement band, the single-coefficient transferability fails.","supporting_citations":[{"cited_title":"Terzis, I","cited_arxiv_id":null,"evidence_quote":"Supplies the micromodel geometry and the micro-PIV velocity data against which the fully resolved 3D reference solution is checked."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the wall-friction drag term that converts the 3D Stokes problem into the quasi-3D model used in the first reduction."},{"cited_title":"Venturoli, E","cited_arxiv_id":null,"evidence_quote":"Precedent for depth-averaged single-phase flow simulation in pseudo-two-dimensional micromodels, supporting the quasi-3D reduction."},{"cited_title":"Laleian, A","cited_arxiv_id":null,"evidence_quote":"Extends depth-averaged flow modeling to micromodels with variable aperture, another precedent for the same reduction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Lists analytical conductance expressions for non-circular capillaries; its lack of a suitable expression for this geometry motivates the numerical upscaling of $k_{ij}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives discrete pore-throat flux from one-dimensional Stokes flow, justifying the pressure-continuity coupling condition at the interface."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the classical permeable-wall slip-condition analogy on which the new per-throat slip parametrization is based."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Precedent for numerically determining an effective slip parameter from local flow data, though for non-creeping flow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the analytical channel-flow velocity profile used to verify the 3D mesh resolution and the integrated flux computation."}],"review_version":1}