{"id":"4b6dfe4a-e40b-4e2d-8c7c-6ba16a8dbc43","arxiv_id":"1908.07454","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A residual-based a posteriori error estimator is proposed and analyzed for a stabilized fully-mixed finite element discretization of the stationary Stokes-Darcy problem.","lead":"This paper derives computable error estimators for a stabilized finite element method that solves the coupled Stokes-Darcy equations, which model surface water flowing into porous media. The estimators are proven to be both reliable and efficient, meaning they can guide adaptive mesh refinement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Reliability proof of Theorem 3.4 omits the nonconforming-interface step: as written it cannot establish (40).","rationale":"The reader's weakest-assumption analysis identifies exactly the point I consider load-bearing: the proof of Theorem 3.4 does not use the stability estimate Theorem 3.3, and no other mechanism is shown for controlling the h^{-1} interface term in ||·||_h. Since the discrete solution U_h is not in the conforming space H (interface normal continuity is only penalized, not enforced), the usual continuous inf-sup argument cannot be applied directly to U-U_h; the residual with conforming test functions is blind to the interface jump. The missing split through H_h∩H is the standard repair, and Theorem 3.3 is the only stated tool for it. The paper also gives no numerical experiments, but that is secondary; the proof gap is the primary issue. I therefore agree with the CONDITIONAL verdict: the estimator is plausible and the gap is probably repairable, but the central reliability theorem is not substantiated as printed.","tokens_in":16911,"tokens_out":17278,"duration_ms":194442,"concrete_test":"Write out a complete proof of Theorem 3.4 using the tools the paper itself lists: choose W_h∈H_h∩H with ||U_h-W_h||_h^2 ≲ J_Γ(U_h,U_h) (Theorem 3.3), bound ||W_h-U_h||_h by the J_Γ term contained in Θ, apply the continuous inf-sup to U-W_h, and estimate the residual (25)-(27) with a Clément interpolation V_h that also accounts for the interface edge terms. Then check term by term that every contribution to ||U-W_h||_h, in particular h^{-1/2}||(u_{fh}-u_{ph})·n_f||_Γ and the pressure/divergence terms, is dominated by Θ+ζ with constants independent of h and δ. If any term cannot be bounded, or if the argument requires Theorem 3.3 to be proved in this paper rather than imported, the reliability claim is not established as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the reliability bound (40), but the proof of Theorem 3.4 is a one-sentence sketch and never confronts the nonconformity of U_h. The norm ||U-U_h||_h contains the interface term h^{-1}||(u_{fh}-u_{ph})·n_f||^2_Γ; since the exact solution satisfies (u_f-u_p)·n_f=0, this is essentially the stabilization term J_Γ(U_h,U_h)/δ. The only tool supplied for controlling such a nonconforming contribution is Theorem 3.3 (Eq. (24)), which bounds the H_h-distance from U_h to the conforming subspace H_h∩H by J_Γ(U_h,U_h). However, the proof of Theorem 3.4 never invokes Theorem 3.3. Instead it appeals to 'the inf-sup condition of L_h' for ||U-U_h||. That inf-sup is a conforming-space statement; for V∈H the interface term c_Γ(v_f-v_p,φ) vanishes, so L_h(U-U_h,V) with conforming test functions cannot see the jump (u_{fh}-u_{ph})·n_f appearing in the norm. A valid proof would have to split U-U_h=(U-W_h)+(W_h-U_h) with W_h∈H_h∩H, bound W_h-U_h using Theorem 3.3, and bound U-W_h by the continuous inf-sup plus the residual estimates (25)-(27). That split is absent, and the Clément step is only asserted. Consequently the derivation of (40) cannot be verified from the text.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a residual-based a posteriori error analysis for a stabilized fully-mixed finite element discretization of the stationary Stokes-Darcy problem. The scheme uses MINI elements in the fluid region, P1 and BDM1 elements in the porous region, and an interface stabilization term J_Gamma. The main results are Theorem 3.4, a global reliability bound ||U - U_h||_h <= C_rel (Theta + zeta), and Theorem 3.5, a local efficiency bound for the indicator Theta_K. The proofs rely on a Helmholtz decomposition, Clément interpolation, bubble-function localization, and three results imported from [44]. No numerical experiments are reported.","tokens_in":17301,"tokens_out":11495,"duration_ms":115387,"significance":"If the stated bounds were proven, the paper would fill a genuine gap: no a posteriori error analysis appears to exist for the stabilized fully-mixed discretization of [58]. The estimator is residual-based and does include the stabilization term, which is appropriate for the nonconforming nature of the scheme. The authors also explicitly identify the auxiliary results from [44] that the analysis needs. However, the manuscript as written does not provide a verifiable proof of the two central estimates: the reliability proof skips the nonconforming-interface step, and the efficiency proof does not localize the stabilization term. The absence of numerical verification further limits confidence. The topic is relevant and the contribution is potentially useful, but the paper requires substantial revision.","major_comments":[{"comment":"The proof of the reliability estimate (40) is not verifiable as written. The proof states that 'the inf-sup condition of L_h' gives ||U - U_h|| <= C sup_{V in H} |L_h(U - U_h, V)| / ||V|| and then says that (41)-(43), Cauchy-Schwarz, and Clément interpolation yield (40). This skips the central difficulty: the target norm ||.||_h contains the interface term h^{-1} ||(u_{f,h} - u_{p,h}) . n_f||_Gamma^2, while the test functions V in H have a continuous normal interface trace, so the inf-sup pairing with conforming V cannot see the nonconforming jump of U_h. The only tool supplied for this jump is Theorem 3.3, Eq. (24), which bounds inf_{W_h in H_h intersect H} ||U_h - W_h||_h^2 by J_Gamma(U_h, U_h), but Theorem 3.4 never invokes Theorem 3.3. A complete proof would need to split U - U_h = (U - W_h) + (W_h - U_h) with W_h in H_h intersect H, estimate W_h - U_h via (24), and U - W_h via the continuous inf-sup and the residual bounds. Because this split is absent, the derivation of (40) is unsupported.","section":"Section 3.3.1, Theorem 3.4"},{"comment":"The residual equation used for reliability is incomplete. After (25) the error is written as L_h(U - U_h, V) = L_h(U - U_h, V - V_h), and it is then expanded into the element residuals (26)-(27) with no interface term involving the discrete velocity jump (u_{f,h} - u_{p,h}) . n_f. However, the stabilized method (18) contains J_Gamma, and the Galerkin orthogonality for L_h = L + J_Gamma produces either L(U - U_h, V_h) = J_Gamma(U_h, V_h) (if L_h denotes L) or an explicit J_Gamma(U - U_h, V - V_h) term (if L_h denotes L + J_Gamma). In either reading, a term of the form delta h^{-1} <(u_{f,h} - u_{p,h}) . n_f, ((v_f - v_{f,h}) - (v_p - v_{p,h})) . n_f> must appear in the residual decomposition. Such a term is absent from (26)-(27) and from (42)-(43), so the stabilization terms in (30) are not connected to the error equation used for reliability. This is a second, independent gap in the proof of (40).","section":"Section 3.2.1, Eqs. (25)-(27)"},{"comment":"The efficiency proof does not establish the bound for the stabilization term. In (53) the authors assert h_E^{1/2} ||[rho g phi_h n_p]||_E + delta h_E^{1/2} h^{-1} ||[(u_{f,h} - u_{p,h}) . n_f]||_E is bounded by local error plus oscillation, but the proof preceding (53) only constructs a bubble test function for the piezometric-head jump [rho g phi_h n_p]_E. No test function is chosen to isolate the velocity jump ((u_{f,h} - u_{p,h}) . n_f) on Gamma, and the statement that 'by regularity Theorem 3.2 the jump of u is zero through all the edges of Omega' does not localize this quantity; the exact interface condition (3) gives (u_f - u_p) . n_f = 0, not an estimate for U_h. Since Theta_{K,p} in (30) contains the stabilization term, Theorem 3.5 cannot be considered proven for that term.","section":"Section 3.3.2, Eq. (53)"},{"comment":"The function space X_f is defined as {v_f in [L^2(Omega_f)]^d : v_f = 0 on Gamma_f}, but the weak form (8) contains a_f(u_f, v_f) = 2 nu (D(u_f), D(v_f)) and the norms used later include |v_f|_{1,Omega_f} = ||nabla v_f||. These expressions are not well-defined for general L^2 functions, and the trace and Korn inequalities quoted in Section 2.2 require H^1 regularity. The space should be a subspace of H^1(Omega_f)^d; as stated, the variational formulation (14) is not well-posed. If this is a typographical omission, it must be corrected and used consistently throughout the paper.","section":"Section 2.2, definition of X_f"}],"minor_comments":[{"comment":"The form L_h is used in (25) and in the proof of Theorem 3.4 but is never defined; please define it explicitly, for example L_h(.,.) = L(.,.) + J_Gamma(.,.), and use it consistently in the Galerkin orthogonality discussion.","section":"Sections 2.3 and 3.3.1"},{"comment":"The Galerkin orthogonality relation is asserted after an unexplained invocation of Theorem 3.2; the intermediate steps needed to subtract (14) and (18) should be written out.","section":"Remark 2.1"},{"comment":"Equation (45) claims ||nabla u_{f,h}||_K = ||nabla(u_f - u_{f,h})||_K, which is false unless u_f is piecewise affine; the intended statement is the divergence identity ||nabla . u_{f,h}||_K = ||nabla . (u_{f,h} - u_f)||_K, using nabla . u_f = 0.","section":"Section 3.3.2, Eq. (45)"},{"comment":"Several discrete-data symbols are not defined before use: f_{p,h} in (43), J_K(v_f, v_f) in the efficiency norm, and the norm ||.||_{h,w} in Theorem 3.5. Please define these quantities explicitly.","section":"Section 3.2.2 and 3.3.2"},{"comment":"There are no numerical experiments; given the proof gaps, a numerical confirmation of the estimator's effectivity would substantially strengthen the paper. The text also contains many typographical and typesetting errors that should be corrected in a revision.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper relies heavily on Theorems 3.1-3.3 from [44], a paper co-authored by the first author. This is not by itself a problem, but since Theorem 3.3 is a nontrivial stability estimate and is essential for the reliability result, the editor may wish to verify that the statement in [44] indeed covers the norm and the finite element spaces used here. The main technical gaps in the present manuscript are fixable in principle, but they require substantial rewriting of the proofs of Theorems 3.4 and 3.5."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper fills a real gap—no a posteriori estimate existed for the stabilized fully-mixed scheme of Yu et al. [58]—and the residual estimator has sensible components, including an interface stabilization term. The efficiency analysis mostly follows the standard bubble-function route. The weak spot is the reliability proof (Theorem 3.4): it is a one-sentence sketch that invokes the inf-sup condition on L_h and the Clément operator, but never confronts the fact that U_h is not in H. The target norm on the left includes the h^{-1} interface term; for test functions in the conforming subspace the coupling term cannot see the jump. The only tool supplied for controlling that is Theorem 3.3 of [44], a stability estimate from the same group, and the proof never uses it. That is a genuine gap, not a stylistic quibble. A correct proof needs the split U-U_h=(U-W_h)+(W_h-U_h) with W_h in H_h∩H, and that is absent. The efficiency proof also has an abrupt step for the stabilization term: the norm defined in Section 3.3.2 does not include the interface normal jump, so estimate (53) does not follow from the previous line. If the authors can fill these holes, the estimator would be a useful tool for adaptive simulations, but as it stands the main theorem is unverified. No numerical experiments either, which would have helped. The reliance on [44] is heavy but not disqualifying; the issue is that the imported theorem is not actually invoked where it needs to be. This is a paper for specialists in a posteriori error estimation for coupled flow; with a rigorous revision it could be a legitimate contribution. Send it to review, but the referee should ask for a complete proof of reliability and at least a 2D numerical experiment.","headline":"A real gap-filling estimator for a Stokes-Darcy scheme, but the reliability proof omits the nonconformity step and needs a serious rewrite before it can be cited.","tokens_in":17770,"tokens_out":5998,"would_cite":false,"duration_ms":58201,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74S05","74S10","74S15","74S20","74S25","74S30"],"pacs":[],"model":"deepseek-v4-flash","headline":"Residual estimator bounds Stokes-Darcy finite-element error from both sides.","keywords":["Stokes-Darcy problem","a posteriori error estimation","stabilized fully-mixed finite element method","MINI element","BDM element","Beavers-Joseph-Saffman condition","reliability and efficiency","isotropic meshes"],"falsifier":"Compute both sides of the imported stability estimate (24) on a family of uniformly refined meshes for a problem with a known exact solution. If the ratio of the infimum over conforming discrete functions to the interface stabilization term grows without bound as $h$ tends to zero, then the assumed estimate fails and the reliability proof of Theorem 3.4 lacks its essential ingredient.","tokens_in":100,"feed_emoji":"🌊","tokens_out":8420,"duration_ms":198542,"temperature":0.7,"pith_summary":"This paper develops a posteriori error estimates for a stabilized fully mixed finite-element discretization of the stationary Stokes-Darcy flow problem. The goal is a computable quantity, built from element residuals and interface stabilization terms, that measures the discretization error without knowledge of the exact solution. The paper proves an upper bound: the error in the mesh-dependent norm $\\|U - U_h\\|_h$ is controlled by the estimator $\\Theta$ plus data oscillation $\\zeta$, and a local lower bound showing each estimator term is controlled by nearby error plus oscillation. Two-sided control of this kind is what an adaptive mesh-refinement algorithm needs to decide where to refine and when to stop.","feed_headline":"Residual estimator bounds Stokes-Darcy finite-element error","feed_subtitle":"New upper and lower bounds on the discretization error let adaptive algorithms refine only where the residual is large.","key_machinery":"The argument turns on the mesh-dependent norm $\\|V\\|_h = (\\|v_f\\|_1^2 + \\|q\\|^2 + \\|v_p\\|_{\\mathrm{div}}^2 + \\|\\psi\\|^2 + h^{-1}\\|(v_f - v_p)\\cdot n_f\\|_\\Gamma^2)^{1/2}$, which weights the interface mismatch by $h^{-1}$. The estimator $\\Theta_K$ assembles element residuals, face jumps, and the stabilization term $J_\\Gamma$, and the proof uses three imported ingredients: a Helmholtz-type decomposition, a regularity result for the exact solution, and the stability estimate $\\inf_{W_h \\in H_h \\cap H}\\|U_h - W_h\\|_h^2 \\lesssim J_\\Gamma(U_h,U_h)$. Bubble functions and quasi-interpolation then convert the residual equation into the upper and lower bounds.","core_discovery":"The central claim is that for the exact solution $U$ of the weak formulation and the discrete solution $U_h$ of the stabilized fully-mixed scheme, the norm $\\|U - U_h\\|_h$ is bounded above by the residual estimator $\\Theta$ plus data oscillation $\\zeta$, while each local indicator $\\Theta_K$ is bounded below by local error plus local oscillation. This is stated as Theorem 3.4 (reliability) and Theorem 3.5 (efficiency). The estimator collects residual terms in the fluid and porous regions, interface terms from the Beavers-Joseph-Saffman condition and normal-flux balance, and the stabilization term measuring the jump of normal velocity across the interface. The proof combines a Helmholtz decomposition, a regularity result, a stability estimate for the stabilization term, and standard bubble-function and interpolation arguments.","pith_inferences":["If the imported stability estimate is the only obstruction, a natural next step is to prove or disprove it directly for the MINI–BDM1/P1 pairing; that would decide whether the reliability constant is truly independent of mesh size.","The same residual structure could be adapted to anisotropic meshes or to unsteady Stokes-Darcy problems, where the $h^{-1}$ interface term would need a time-dependent counterpart.","The paper does not run adaptive experiments; a concrete test of the theory would be to compare effectivity indices on a sequence of adaptively refined meshes against the predicted constants."],"forward_implications":["The estimator $\\Theta + \\zeta$ is computable from the discrete solution and the data, so it can drive adaptive mesh refinement without solving the continuous problem.","Local indicators $\\Theta_K$ identify which elements contribute most to the error, concentrating refinement near the fluid-porous interface and regions with large residuals.","The efficiency bound shows that stopping when the estimator is small does not miss large local errors, up to data oscillation.","Because the analysis covers dimensions two and three on isotropic meshes, the same estimator applies to practical two- and three-dimensional Stokes-Darcy simulations."],"supporting_citations":[{"why":"Defines the stabilized fully-mixed MINI–BDM1/P1 discretization whose error this paper estimates.","marker":"[58]"},{"why":"Supplies the Helmholtz decomposition, the regularity result, and the stability estimate (24) on which the reliability proof depends.","marker":"[44]"},{"why":"Provides the bubble-function and extension arguments used to prove local efficiency.","marker":"[54]"},{"why":"Provides the quasi-interpolation operator with local approximation properties used for reliability.","marker":"[17]"},{"why":"Supplies the inverse-inequality localization technique for lower error bounds in mixed methods.","marker":"[12]"},{"why":"Gives an earlier residual-based estimator for a fully-mixed Stokes-Darcy formulation that the present estimator extends to the stabilized scheme.","marker":"[28]"}],"fun_headline_variants":["Residual estimator for Stokes-Darcy: reliable and efficient bounds","A posteriori error bounds for stabilized Stokes-Darcy scheme","Adaptive Stokes-Darcy discretization proven with residual estimator","Stokes-Darcy error estimator enables adaptive refinement"],"cache_read_input_tokens":19840,"weakest_assumption_plain":"The reliability theorem depends on an estimate imported from an earlier paper: the discrete solution can be mimicked, in the norm used for the error, by a discrete function with exact interface continuity, with error no larger than the stabilization term. If that estimate is false, the stated upper bound is unsupported.","fun_headline_variants_meta":{"raw":{"variants":["Residual estimator for Stokes-Darcy: reliable and efficient bounds","A posteriori error bounds for stabilized Stokes-Darcy scheme","Adaptive Stokes-Darcy discretization proven with residual estimator","Stokes-Darcy error estimator enables adaptive refinement"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000425,"raw_usage":{"total_tokens":2116,"prompt_tokens":819,"completion_tokens":1297,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":435,"completion_tokens_details":{"reasoning_tokens":1227}},"tokens_in":435,"tokens_out":1297,"duration_ms":9772,"temperature":1.0,"reasoning_tokens":1227,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:40:22.484486+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute both sides of the imported stability estimate (24) on a family of uniformly refined meshes for a problem with a known exact solution. If the ratio of the infimum over conforming discrete functions to the interface stabilization term grows without bound as $h$ tends to zero, then the assumed estimate fails and the reliability proof of Theorem 3.4 lacks its essential ingredient.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the stabilized fully-mixed MINI–BDM1/P1 discretization whose error this paper estimates."},{"cited_title":"Nicaise, B","cited_arxiv_id":null,"evidence_quote":"Supplies the Helmholtz decomposition, the regularity result, and the stability estimate (24) on which the reliability proof depends."},{"cited_title":"Verf¨ urth","cited_arxiv_id":null,"evidence_quote":"Provides the bubble-function and extension arguments used to prove local efficiency."},{"cited_title":"Cl´ ement","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-interpolation operator with local approximation properties used for reliability."},{"cited_title":"Carstensen","cited_arxiv_id":null,"evidence_quote":"Supplies the inverse-inequality localization technique for lower error bounds in mixed methods."},{"cited_title":"Gatica, R","cited_arxiv_id":null,"evidence_quote":"Gives an earlier residual-based estimator for a fully-mixed Stokes-Darcy formulation that the present estimator extends to the stabilized scheme."}],"review_version":1}