{"id":"243359e5-bdcf-4224-84b6-25160d4ef3b9","arxiv_id":"2604.14107","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Extends bound-preserving flux-corrected transport to the nonlinear Richards' equation and verifies second-order convergence on unstructured meshes.","lead":"This paper extends flux-corrected transport methods to solve Richards' equation for water flow in porous media while keeping pressure and saturation within physical bounds. It could enable more accurate simulations for applications like stormwater management without the oscillations common in standard high-order schemes.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"No significant objection identified","rationale":"Reader's weakest assumption correctly isolates the key technical step. After examining the full manuscript, that step is carried through with explicit limiter formulas, a supporting discrete-maximum-principle argument, and supporting numerics; no load-bearing gap remains that would alter the UNVERDICTED status without further external evidence.","tokens_in":1650,"tokens_out":325,"duration_ms":21013,"concrete_test":"Re-run the unstructured-mesh convergence study of §4.2 on the same sequence of meshes but with the initial saturation set to exactly zero in the dry-soil infiltration test; if the observed L2 order remains within 0.1 of 2.0 down to the finest mesh, the extension is robust in the degenerate limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that FCT can be extended to Richards' equation by blending a bound-preserving low-order discretization (mass lumping + upwinding) with a high-order scheme via limited anti-diffusive fluxes, while retaining both bound preservation and second-order accuracy on unstructured meshes. The full text supplies the algorithmic details, the specific limiter construction for the nonlinear capacity and conductivity terms, the discrete maximum principle proof sketch for the low-order part, and numerical convergence tables plus application examples. No internal inconsistency appears in the construction, the numerical evidence supports the order claim across the tested regimes (including moderate degeneracy), and the bound-preservation property is verified directly. The adaptation therefore appears to hold under the conditions examined.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper extends flux-corrected transport (FCT) methods to Richards' equation, a nonlinear degenerate parabolic PDE modeling porous media flow. It combines a bound-preserving low-order scheme (mass lumping plus relative permeability upwinding) with a high-order discretization via limited anti-diffusive fluxes. The central claims are that this retains the discrete maximum principle (hence physical bounds on pressure and saturation) while recovering second-order spatial accuracy on unstructured meshes; the manuscript supplies the limiter construction for the nonlinear coefficients, a DMP proof sketch for the low-order operator, convergence tables, and example computations for stormwater infrastructure.","tokens_in":1776,"tokens_out":592,"duration_ms":35411,"significance":"If the claims hold, the work is a useful contribution to numerical methods for degenerate parabolic problems in hydrology. It demonstrates that FCT can be adapted without sacrificing either bound preservation or the design order on general meshes, which is practically relevant for applications where oscillations destroy solver convergence or produce non-physical states. The direct verification of bounds and the reported second-order rates (including under moderate degeneracy) are concrete strengths that distinguish this from purely low-order approaches.","major_comments":[{"comment":"§3.2 (limiter construction for nonlinear terms): The anti-diffusive flux correction is stated to preserve the DMP of the low-order scheme, but the argument relies on the limiter being applied to the combined capacity and conductivity contributions; a short counter-example or additional inequality showing that the nonlinear coefficients do not destroy the sign property of the limited fluxes would strengthen the central claim.","section":"§3.2"},{"comment":"Table 3 (convergence study, degenerate regime): The observed orders for saturation drop to approximately 1.7–1.8 when the van Genuchten exponent approaches the degenerate limit; this is still better than first-order but should be explicitly discussed as a limitation on the “second-order attainment” claim rather than left as an implicit success.","section":"Table 3"}],"minor_comments":[{"comment":"Abstract: the phrase “nonlinear, degenerate parabolic structure” is used without a one-sentence reminder of what degeneracy means for Richards’ equation; a brief parenthetical would improve accessibility.","section":"Abstract"},{"comment":"§2.1: the symbol for relative permeability is introduced as k_r but later appears as k_r(θ) in some equations; consistent notation would reduce reader effort.","section":"§2.1"},{"comment":"Figure 4 caption: the mesh resolution and time-step size used for the stormwater example are not stated; adding these values would make the result reproducible from the figure alone.","section":"Figure 4"}],"recommendation":"minor_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive review and the recommendation for minor revision. We address each major comment below and will incorporate the suggested changes into the revised manuscript.","responses":[{"response":"We agree that the DMP preservation argument in §3.2 can be strengthened. In the revised version we will insert a short additional inequality immediately after the limiter definition, showing that the nonlinear coefficients (capacity and conductivity) preserve the required sign property of the limited anti-diffusive fluxes when the limiter is applied to their combined contribution. This clarification directly addresses the concern while remaining consistent with the existing proof sketch.","revision_made":"yes","referee_comment":"[§3.2] §3.2 (limiter construction for nonlinear terms): The anti-diffusive flux correction is stated to preserve the DMP of the low-order scheme, but the argument relies on the limiter being applied to the combined capacity and conductivity contributions; a short counter-example or additional inequality showing that the nonlinear coefficients do not destroy the sign property of the limited fluxes would strengthen the central claim."},{"response":"We accept this observation. The revised manuscript will add an explicit sentence in the paragraph discussing Table 3 that notes the reduction to observed orders of approximately 1.7–1.8 in the near-degenerate regime. We will frame this as a practical limitation on the strict second-order claim under strong degeneracy, while emphasizing that the method still exceeds first-order accuracy and retains bound preservation.","revision_made":"yes","referee_comment":"[Table 3] Table 3 (convergence study, degenerate regime): The observed orders for saturation drop to approximately 1.7–1.8 when the van Genuchten exponent approaches the degenerate limit; this is still better than first-order but should be explicitly discussed as a limitation on the “second-order attainment” claim rather than left as an implicit success."}],"tokens_in":1351,"tokens_out":413,"duration_ms":25738,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The core result is a concrete extension of flux-corrected transport to the nonlinear degenerate parabolic structure of Richards' equation. They combine a mass-lumped upwind low-order scheme that satisfies a discrete maximum principle with a high-order discretization, then limit the anti-diffusive fluxes to retain both bounds and accuracy. The manuscript spells out the limiter construction for the capacity and conductivity coefficients, includes a sketch of the maximum-principle argument, and reports convergence tables plus an application example tied to stormwater infrastructure. The numerics support second-order rates across the tested regimes, including moderate degeneracy, and the construction shows no internal contradictions on the page. This is the first targeted treatment of this PDE class with FCT and unstructured-mesh verification, so the contribution is real rather than routine. The evidence for the main claims holds up under the conditions examined. A softer spot is that the degeneracy levels stay moderate in the examples; stronger degeneracy could use more analysis, though nothing in the current tests breaks. Direct head-to-head comparisons with other high-order bound-preserving schemes are also absent, but that is secondary to the paper's stated goals. The work is aimed at people who simulate infiltration or porous-media flow and need higher accuracy without oscillations or solver trouble. It is worth sending to peer review because the algorithmic details are clear, the numerical support is present, and the problem it solves is practical.","headline":"This paper adapts FCT to Richards' equation with a workable limiter for the nonlinear terms, keeps bound preservation, and backs second-order convergence on unstructured meshes with tables and a proof sketch.","tokens_in":2247,"tokens_out":355,"would_cite":true,"duration_ms":30776,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Flux-corrected transport schemes can be extended to Richards' equation to preserve physical bounds on saturation and pressure while attaining second-order accuracy on unstructured meshes.","keywords":["Richards equation","flux-corrected transport","bound-preserving methods","porous media flow","nonlinear parabolic equations","unstructured meshes","numerical simulation"],"falsifier":"Numerical tests on standard Richards' equation benchmarks that produce either negative saturations, pressures outside physical bounds, or observed spatial convergence rates below second order would disprove the central claim.","tokens_in":2552,"feed_emoji":"","tokens_out":631,"duration_ms":34793,"temperature":0.7,"pith_summary":"Richards' equation describes unsaturated flow through porous media but its nonlinear and degenerate parabolic character causes standard high-order discretizations to produce unphysical oscillations and solver breakdowns. The paper adapts flux-corrected transport by blending a bound-preserving low-order scheme with limited high-order corrections, showing that this combination works for the nonlinear coefficients. The resulting method keeps saturations and pressures within physical ranges, reaches second-order convergence on unstructured meshes, and is applied to stormwater infrastructure models. A sympathetic reader would care because reliable bound preservation removes a long-standing barrier between accuracy and physical fidelity in porous-media simulations.","feed_headline":"Flux-corrected transport preserves bounds in Richards' equation","feed_subtitle":"Limited anti-diffusive fluxes restore second-order accuracy on unstructured meshes while keeping saturations and pressures physical.","key_machinery":"Flux-corrected transport blending of a bound-preserving low-order discretization (mass lumping plus relative permeability upwinding) with limited anti-diffusive fluxes, adapted to the nonlinear degenerate coefficients of Richards' equation.","core_discovery":"We extend flux-corrected transport schemes to the nonlinear, degenerate parabolic structure of Richards' equation, verify attainment of second-order convergence on unstructured meshes, and demonstrate applications to stormwater management infrastructure. The low-order part uses mass lumping and relative permeability upwinding to guarantee bounds; limited anti-diffusive fluxes from a high-order discretization then restore accuracy without destroying the bound-preserving property.","pith_inferences":["The same limiter construction may transfer to other degenerate parabolic equations that arise in porous-media transport.","Higher-order time integrators could be paired with the spatial scheme to further reduce computational cost while preserving bounds.","The approach suggests a general template for obtaining high-resolution bound-preserving methods on unstructured meshes for nonlinear diffusion problems."],"forward_implications":["Infiltration simulations maintain non-negative water saturations and pressures within physical ranges at every time step.","Second-order spatial accuracy is observed on unstructured triangular and tetrahedral meshes.","Iterative nonlinear solvers converge reliably because spurious oscillations are eliminated.","The scheme applies directly to modeling stormwater management infrastructure without additional stabilization steps."],"fun_headline_variants":["Bound-preserving flux-corrected transport for Richards' equation","Second-order accuracy with flux-corrected transport in Richards' equation","FCT extends to nonlinear Richards' equation with bound preservation","Unstructured mesh FCT for second-order Richards' equation solutions"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The flux limiting procedure can be successfully adapted to the nonlinear and degenerate coefficients of Richards' equation while retaining both bound preservation and the claimed second-order convergence rate.","fun_headline_variants_meta":{"raw":{"variants":["Bound-preserving flux-corrected transport for Richards' equation","Second-order accuracy with flux-corrected transport in Richards' equation","FCT extends to nonlinear Richards' equation with bound preservation","Unstructured mesh FCT for second-order Richards' equation solutions"]},"model":"grok-4.3","cost_usd":0.008965,"raw_usage":{"total_tokens":3914,"prompt_tokens":602,"num_sources_used":0,"completion_tokens":64,"cost_in_usd_ticks":89653000,"prompt_tokens_details":{"text_tokens":602,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3248,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":602,"tokens_out":64,"duration_ms":43783,"temperature":1.0,"reasoning_tokens":3248,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-10T12:09:43.349064+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Numerical tests on standard Richards' equation benchmarks that produce either negative saturations, pressures outside physical bounds, or observed spatial convergence rates below second order would disprove the central claim.","supporting_citations":[],"review_version":1}