{"id":"0929ac84-fd11-4465-86f2-0ec3aad17fff","arxiv_id":"2509.17899","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"For inhomogeneous Dirichlet Stokes data, quasi-optimal and pressure-robust error estimates for the Scott-Vogelius element and convergence of the iterated penalty method are established, contingent on a compatibility condition and an unpublished Fortin-operator construction.","lead":"This paper proves error bounds for finite element solutions of the Stokes flow equations when the velocity is prescribed on the boundary, and for the Scott-Vogelius element the velocity error bound does not depend on pressure. It also analyzes an iterative penalty solver, showing convergence and that after enough iterations the method becomes pressure-robust.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 4.2's pressure estimate reverses inequality (42): the step bounding ||div e_{h,i}|| by β||∇e_{h,i}|| is invalid, so the stated IPM pressure rate does not follow.","rationale":"The paper's central theoretical contribution is the pressure-robust quasi-optimal estimate for Scott–Vogelius with inhomogeneous Dirichlet data (Theorem 3.4) and the asymptotic pressure robustness of the IPM (Theorem 4.5). The Fortin operator assumption is the foundational external dependency, and the paper is transparent that Proposition 3.3 is deferred to the unpublished companion [ET]. I did not find an internal inconsistency in the use of that assumption, so my main concern is the concrete direction error in Lemma 4.2 Step 3. Inequality (42) gives a lower bound on ||div v_h||, not an upper bound, so replacing ρ||div e_i|| by ρβ||∇e_i|| is invalid. This does not destroy the velocity contraction result or the velocity asymptotic pressure-robustness, which depend only on Step 1 of Lemma 4.2. It does, however, invalidate the stated pressure estimate (45) and the corresponding pressure estimate in Theorem 4.5 as written. The numerical experiments are consistent with the main qualitative claims and Theorem 4.6 is independently argued, but the manuscript should be revised to correct Lemma 4.2 and to either prove or cite Proposition 3.3. This supports the reader's CONDITIONAL verdict without moving it further.","tokens_in":32439,"tokens_out":17526,"duration_ms":152849,"concrete_test":"Re-derive Lemma 4.2 Step 3 with the invalid step removed: bound ρ||div e_{h,i}|| using (44) instead of (42), then compare the resulting coefficient in the θ^{i-1} bound to (ν+ρβ)/β. For ν=1, ρ=100, β=1, the claim gives 101 θ^{i-1}||∇e_1||, while the repaired derivation gives approximately 506 θ^{i-1}||∇e_1||; this algebraic mismatch settles that (45) is false as written.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Lemma 4.2, Step 3 (Eqs. (50)–(52)) bounds ρ||div e_{h,i}||_{L2} by ρβ||∇e_{h,i}||_{L2}. The only cited inequality, (42), says β||∇v_h|| ≤ ||div v_h|| for v_h ∈ V_{h,div}^\\perp — the reverse direction. Hence (45) is not proven. A valid repair combines (49) and (44): ||π|| ≤ (ν/β)||∇e_i|| + ρ||div e_i|| ≤ [ν/β θ + (1/2)√(νρ)] θ^{i-2}||∇e_1||, which is not equal to (ν+ρβ)/β θ^{i-1}||∇e_1||. So Theorem 4.5's pressure estimate is unsupported as stated. The velocity contraction (43) and the velocity part of Theorem 4.5 are unaffected. Separately, Theorem 3.4 and the IPM theorem ultimately rely on Proposition 3.3 from the unpublished companion [ET]; that dependency remains unresolved in this manuscript.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies mixed finite element methods for the Stokes equations with inhomogeneous Dirichlet boundary conditions, with emphasis on the Scott–Vogelius element. The main theoretical tools are Fortin-type operators that preserve the divergence, discrete traces, and zero mean normal traces. Under these assumptions the paper proves quasi-optimal a priori error estimates (Theorem 2.13) and, for Scott–Vogelius, a pressure-robust velocity estimate (Theorem 3.4). It then analyzes the iterated penalty method (IPM), proving convergence, a velocity contraction bound, an asymptotic pressure-robustness estimate, and monotone decay of the divergence norm (Lemma 4.2, Theorems 4.5 and 4.6). Numerical experiments in 2D illustrate the role of the compatibility condition and of mesh modifications removing (nearly) singular vertices.","tokens_in":32809,"tokens_out":7738,"duration_ms":69549,"significance":"If the results are valid, the paper is a useful contribution: it extends pressure-robust, exactly divergence-free discretizations to inhomogeneous Dirichlet data, clarifies the required compatibility condition, and provides explicit convergence rates for the IPM in this setting. The numerical experiments are well targeted and illustrate the main theoretical points. However, the central Fortin-operator existence result is deferred to an unpublished companion paper, and one key pressure estimate in the IPM analysis uses an inequality in the wrong direction. These issues must be resolved before the main claims can be accepted.","major_comments":[{"comment":"The step 'applying (42)' is invalid. Inequality (42) states β||∇v_h|| ≤ ||div v_h|| for v_h ∈ V_{h,div}^⊥, i.e., it bounds the gradient from above by the divergence. It cannot be used to conclude ρ||div e_{h,i}|| ≤ ρβ||∇e_{h,i}||, which is the reverse direction. Consequently the displayed pressure estimate (45) does not follow. A repair using the already-proved divergence bound (44) would give a pressure estimate of the form [ (ν/β) θ^{i-1} + (1/2)√(νρ) θ^{i-2} ] ||∇e_{h,1}||, not the stated (ν+ρβ)/β θ^{i-1}||∇e_{h,1}||. Since Theorem 4.5's pressure inequality relies directly on (45), that part of the theorem is unsupported as stated. The velocity contraction (43) and divergence bound (44) are unaffected.","section":"Lemma 4.2, Step 3 (Eq. (52))"},{"comment":"The existence of a Fortin operator preserving divergence, discrete traces, and zero mean normal traces with uniform H^1-stability is the load-bearing assumption for Theorem 3.4 and, through the first-iterate estimate, for Theorem 4.5. Proposition 3.3 states this existence for the Scott–Vogelius element but defers the proof entirely to [ET], an unpublished companion paper by two of the same authors. The manuscript therefore is not self-contained and the central claim cannot be verified from the submitted material. The authors should either provide a full construction and proof (e.g., in an appendix) or replace the reference with a published or otherwise publicly available source.","section":"Proposition 3.3 and Assumption 2.10"}],"minor_comments":[{"comment":"Typo: 'error stimates' should be 'error estimates'.","section":"Abstract"},{"comment":"The column alignment is confusing: each row appears to contain five numeric entries for four error columns plus N, and the values cited in the text (e.g., ||div u_h||=6.74e-12 for corrected data) do not clearly match the printed headers. Please relabel the tables or insert explicit column separators.","section":"Tables 1 and 2"},{"comment":"The numerical implementation uses the interpolation operator Ĩ_h defined on C(Ω)^d, whereas Proposition 3.3 and Assumption 2.14 concern an operator I_h on H^1(Ω)^d. The text notes this difference, but it would help to state explicitly that the manufactured boundary data are smooth enough that the C^0-based construction is covered by the same arguments.","section":"Section 3.2.2"},{"comment":"The wording 'preserves zero mean normal traces' is slightly stronger than the stated property Π_h(H^1_∼(Ω)) ⊂ X_{h,∼}. Consider rephrasing to 'maps zero-mean-normal-trace functions into the corresponding discrete subspace' to avoid possible confusion.","section":"Assumption 2.10 (iiib)"},{"comment":"The statement that tr(g_h)=tr(Π_h g) and tr(g_h)=tr(I_h g) 'imply g_h ∈ X_{h,∼}' is true only because the exact data g is assumed to lie in H^1_∼(Ω) and the operators map H^1_∼(Ω) into X_{h,∼}. This could be made explicit in the remark.","section":"Remark 2.16(a)"}],"recommendation":"major_revision","confidential_remarks":"The reverse-inequality error in Lemma 4.2 is concrete and directly affects a stated theorem; it is fixable but requires a revised constant and rate. The dependence on the unpublished companion [ET] for the central Fortin operator is a serious self-containment issue; I recommend asking the authors to include the full proof or to cite a available source. If the companion remains unavailable, the main theorems should be recast as conditional on an explicitly stated assumption rather than as unconditional results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know before you read this one. The paper does real work: Theorem 2.13 gives a clean quasi-optimal estimate for general mixed methods with inhomogeneous Dirichlet data, improving on the constrained infimum in [BGHRR24], and Theorem 3.4 gives a pressure-robust velocity estimate for Scott-Vogelius. The IPM velocity contraction is simple and correct, and the monotone divergence decay is a nice bonus. The numerics are consistent with the theory and the compatibility issue is handled carefully.\n\nThe soft spots matter. First, Lemma 4.2, Step 3, uses inequality (42) in the wrong direction. That inequality bounds ||grad v|| from above by (1/beta)||div v||; the text uses it to bound ||div e|| by beta||grad e||. That direction does not follow, so (45) is not proven and the pressure part of Theorem 4.5 is unsupported as stated. The velocity contraction (43) and divergence estimate (44) survive; this is a local error in the pressure-rate bookkeeping, not a collapse of the whole IPM argument, but it needs a repair before the pressure claims can be taken seriously.\n\nSecond, the central trace-preserving Fortin operator, Assumption 2.10 / Proposition 3.3, is deferred to [ET], an in-preparation companion by two of the authors. Theorems 3.4 and the IPM asymptotic pressure-robustness claim depend on that operator existing as stated. The paper is transparent about the dependency, so this is a burden, not circularity, but it means the headline theorem is conditional on unpublished work.\n\nWhat is genuinely good: the compatibility condition for boundary data is treated rigorously, the paper clearly separates what is proved here from what is cited, and the velocity-side analysis in Section 4.1 is sound. The stress-test note holds up; the reader's conditional verdict is fair.\n\nWho is this for: people working on divergence-free mixed methods and Uzawa-type solvers for Stokes and Navier-Stokes will want to see it, especially for the clean Theorem 2.13 and the velocity contraction. But the paper should not be accepted as-is. The pressure estimate needs a correct proof or a clearly weaker statement, and the Fortin operator claim needs to be either proved here or the dependency flagged in a way that lets readers evaluate the consequences.\n\nSend it to peer review, but ask reviewers to check Lemma 4.2 and the [ET] dependency carefully. It deserves serious referee time; it is not ready as-is.","headline":"Useful, honest conditional contribution: quasi-optimal inhomogeneous-Dirichlet estimates for Scott-Vogelius and a clean velocity contraction for IPM, but the pressure iterate bound in Lemma 4.2 reverses (42) and the central Fortin operator lives in an unpublished companion.","tokens_in":33247,"tokens_out":2446,"would_cite":true,"duration_ms":205465,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N12","65N15","65N30","65F10","76D07","76M10"],"pacs":[],"model":"deepseek-v4-flash","headline":"For Stokes flow with non-homogeneous Dirichlet data, the Scott–Vogelius element keeps its pressure-robust velocity error, provided the boundary data satisfies a zero-mean-normal-trace compatibility condition.","keywords":["Stokes problem","inhomogeneous Dirichlet boundary conditions","Scott–Vogelius element","pressure robustness","iterated penalty method","Fortin operator","divergence-free","quasi-optimality"],"falsifier":"Run the Scott–Vogelius k=4 method on a mesh with no singular vertices, impose compatible boundary data g_h ∈ X_{h,∼}, and scale the pressure loading by a factor Ra; if the H1 velocity error grows with Ra, or the IPM iterates fail to reach the predicted contraction floor, Theorems 3.4 and 4.5 would be refuted.","tokens_in":32375,"feed_emoji":"🌊","tokens_out":5312,"duration_ms":41587,"temperature":0.7,"pith_summary":"The paper extends quasi-optimal error analysis of mixed finite element methods for the Stokes problem to the case of inhomogeneous Dirichlet boundary conditions, where the boundary data must satisfy a compatibility condition. For the Scott–Vogelius element, which enforces the divergence constraint exactly, it proves that the velocity error is bounded only by the best approximation error of the velocity—no pressure term appears—so the method remains pressure-robust even with non-zero boundary data. It also analyses the iterated penalty method (an Uzawa-type iteration) and shows it converges with an explicit contraction factor and becomes asymptotically pressure-robust as the iteration count grows. A single-face correction of the boundary interpolation enforces the needed zero-mean normal trace condition, and numerical experiments confirm the theory. The key tool is a modified Fortin operator that preserves divergence, discrete traces, and zero-mean normal traces simultaneously.","feed_headline":"Scott-Vogelius stays pressure-robust under non-homogeneous data","feed_subtitle":"Trace-preserving Fortin operator and zero-mean-normal-trace boundary data remove pressure from the velocity error.","key_machinery":"The load-bearing object is a trace-preserving Fortin operator Π_h : H^1(Ω)^d → X_h that (i) preserves the divergence against the discrete pressure space, (ii) is H^1-stable, and (iii) preserves discrete traces and maps the zero-mean-normal-trace space H^1_∼(Ω) into X_{h,∼}. In the proof of Theorem 3.4 this operator is used to repair an arbitrary v_h into a function z_h that is discretely divergence-free and matches the boundary data, so the constrained best-approximation infimum reduces to the unconstrained one. For Scott–Vogelius the analysis additionally assumes div V_h = div X_{h,∼} (Assumption 3.1), which holds in 2D when the mesh has no boundary singular vertices. The iterated penalty m","core_discovery":"On its own terms, the central claim is Theorem 3.4: for the Scott–Vogelius discretization of the Stokes problem with compatible boundary data, the H1-seminorm velocity error is at most 2(1+c_F) times the best approximation error over the velocity space, uniform in mesh size, with no pressure contribution. This is the first such pressure-robust quasi-optimality result for inhomogeneous Dirichlet conditions for an exactly divergence-free element. The companion result, Theorem 4.5, shows that the iterated penalty method produces velocity iterates whose error is the sum of the same velocity best-approximation term plus an iteration-dependent term that decays geometrically with factor θ=ν/(ν+ρβ²)","pith_inferences":["Because the linear Stokes analysis can be reduced to homogeneous boundary data via a divergence-free extension, the real payoff of this framework is for nonlinear equations such as Navier–Stokes, where such a reduction fails; the trace-preserving Fortin operator is the natural tool to carry the inhomogeneous setting there.","The single-face correction of the boundary interpolation suggests a minimally invasive implementation recipe: take any existing Lagrange boundary interpolation and correct one face to restore the zero-mean normal trace, so existing codes could adopt the compatible treatment without a full rewrite.","The contraction factor's dependence on 1/β means meshes with nearly singular vertices can make the IPM stagnate in practice; the numerical experiments show this plateau, implying mesh quality and iteration count interact more strongly than the asymptotic theory alone suggests.","A testable prediction: on a sequence of meshes where a boundary vertex angle tends to zero, the number of IPM iterations to reach a fixed divergence tolerance should grow roughly like 1/β², providing a sharp experimental check of the contraction bound."],"forward_implications":["With compatible boundary data, inhomogeneous Dirichlet conditions can be imposed on Scott–Vogelius discretizations without sacrificing the pressure-robust velocity error that makes exactly divergence-free elements valuable.","The iterated penalty method converges for inhomogeneous boundary data with a contraction factor θ = ν/(ν+ρβ²), and the divergence norm of the velocity iterates decreases monotonically, giving a stopping criterion that actually terminates.","The compatibility condition is not a technicality: without enforcing zero mean normal trace, the discrete velocity is not exactly divergence-free and the error estimates lose their pressure-robust form.","Other exactly divergence-free elements, such as the Falk–Neilan and Guzmán–Neilan elements, gain the same pressure-robust inhomogeneous-boundary estimates whenever a Fortin operator with the three trace/divergence properties can be constructed.","Local mesh modification that removes (nearly) singular boundary vertices improves the inf-sup constant, which directly accelerates the IPM contraction, as confirmed by the numerical experiments."],"fun_headline_variants":["Scott-Vogelius pressure-robust for inhomogeneous Dirichlet Stokes","Trace-preserving Fortin yields pressure-robust Stokes estimates","Iterated penalty method converges pressure-robustly","Compatible boundary data key to pressure-robust Stokes","First pressure-robust quasi-optimality for Scott-Vogelius"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The results rest on the existence of a uniformly H1-stable Fortin operator that preserves the divergence, the discrete traces, and the zero-mean normal traces; for Scott–Vogelius this existence is deferred to an unpublished companion paper, and in 2D it requires meshes with no boundary singular vertices.","fun_headline_variants_meta":{"raw":{"variants":["Scott-Vogelius pressure-robust for inhomogeneous Dirichlet Stokes","Trace-preserving Fortin yields pressure-robust Stokes estimates","Iterated penalty method converges pressure-robustly","Compatible boundary data key to pressure-robust Stokes","First pressure-robust quasi-optimality for Scott-Vogelius"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000906,"raw_usage":{"total_tokens":3689,"prompt_tokens":653,"completion_tokens":3036,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":397,"completion_tokens_details":{"reasoning_tokens":2949}},"tokens_in":397,"tokens_out":3036,"duration_ms":17401,"temperature":1.0,"reasoning_tokens":2949,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T15:49:50.778253+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the Scott–Vogelius k=4 method on a mesh with no singular vertices, impose compatible boundary data g_h ∈ X_{h,∼}, and scale the pressure loading by a factor Ra; if the H1 velocity error grows with Ra, or the IPM iterates fail to reach the predicted contraction floor, Theorems 3.4 and 4.5 would be refuted.","supporting_citations":[],"review_version":1}