{"id":"a4dde6b5-bad4-43de-80fc-724168fc12e7","arxiv_id":"1908.02417","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A new homogenization method approximates one-dimensional heterogeneous diffusion with an advection-diffusion equation whose effective diffusivity and velocity are fitted to match integral averages of the heterogeneous solution, reducing error in numerical tests.","lead":"This paper adds an advection term to the constant-coefficient diffusion equation commonly used to approximate transport in materials with spatially varying properties. The authors show that the extra term improves agreement with the detailed model in nine test cases, but the improvement comes from fitting the two effective coefficients to the detailed solution.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No out-of-sample test: effective coefficients are fitted to each problem via (12), (17), so the observed improvement may come from calibration rather than from the advection term itself.","rationale":"The reader's weakest assumption concerns the ad hoc nature and sufficiency of the two moment constraints. My stress-test identifies a sharper version of that concern: the constraints are not merely ad hoc, they are fitted to the exact problem used for evaluation. Because w(x) depends on the initial condition through r=f-s in Eq. (20), the effective coefficients are solution-specific, making the comparison in Table II an in-sample comparison. This does not disprove the paper's claim, and the physical motivation from Eq. (4) plus the reproducible code are real supporting evidence. However, it means the reported numbers cannot yet establish that advection, rather than calibration, improves homogenized models. A concrete fair-comparison test, using a diffusion-only model whose single coefficient is also calibrated by the same moment constraint, would settle the main ambiguity. The reader's CONDITIONAL verdict remains appropriate; no change is needed, but the stated condition should include an out-of-sample or fair-baseline check.","tokens_in":8274,"tokens_out":7430,"duration_ms":86869,"concrete_test":"For each test case A-I, replace the harmonic-mean baseline with a diffusion-only model whose single coefficient ~D_eff is calibrated using the same moment-matching philosophy, e.g. by solving the temporal constraint ∫W dx = ∫w dx in (17) with v_eff=0. Recompute the mean absolute errors (33)-(34). If the advection model's advantage largely disappears, the improvement in Table II is attributable to the extra calibrated parameter rather than to the advection term. If the advantage persists, the central claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that adding an advection term to the homogenized equation (9)-(11) improves accuracy over the standard diffusion-only model (29)-(32). The evidence in Table II is computed in-sample: for every test case, D_eff and v_eff are obtained by solving the two moment constraints (12) and (17) against the same heterogeneous problem later used to compute the errors, whereas the baseline D_eff is the harmonic average (32), which is not fitted to the solution. This confounds two effects: the presence of an advection term and the addition of a second, solution-calibrated parameter. The problem is compounded by the fact that the calibration data themselves depend on the whole initial/boundary-value problem: w(x) solves (20)-(22) with r(x)=f(x)-s(x), so the effective coefficients depend on the initial condition f(x) as well as on g0 and gL, although the conclusions state only boundary-condition dependence. No out-of-sample test, transferability check, or error bound is provided. The claim may well be true, but the reported experiments do not distinguish a genuinely better homogenized model from a two-parameter fit to each benchmark solution.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an alternative homogenized model for one-dimensional diffusion with a spatially varying diffusivity. Instead of the standard diffusion-only homogenized equation with harmonic-average diffusivity, the authors assume a constant-coefficient advection-diffusion equation, Eqs. (9)-(11), with effective diffusivity D_eff and effective velocity v_eff. The two coefficients are determined by two integral constraints: equality of the spatial averages of the steady-state solutions (12), and equality of the spatial averages of the time-integrated deviations from steady state (17). These constraints require solving two boundary value problems over the heterogeneous medium, Eqs. (13)-(14) and (20)-(22), and two over the homogenized equation, Eqs. (15)-(16) and (23)-(25), yielding the nonlinear system (26)-(28). Nine test cases are presented in Table I, and Table II reports smaller mean absolute errors for the new advection-diffusion homogenized model than for the standard diffusion-only model in all cases. The paper concludes that including an advection term improves homogenized models of continuum diffusion in one dimension.","tokens_in":8508,"tokens_out":4893,"duration_ms":53521,"significance":"If the central claim holds, the method provides a simple, practical way to capture directional, advective-like effects of spatial heterogeneity with constant coefficients, extending the author's earlier work on stochastic diffusion. The paper is clearly written and the numerical code is publicly available, which is a strength. The reported improvement is consistent across all nine test cases, including cases with time-dependent boundary conditions and non-uniform initial data. However, the evidence as presented is not fully convincing because the effective coefficients are calibrated to each benchmark problem before the errors are computed, and no out-of-sample test or error bound is provided. The significance is therefore conditional on the ability to demonstrate that the improvement is not merely an artifact of fitting two parameters to each solution.","major_comments":[{"comment":"The numerical comparison is in-sample: for every test case, D_eff and v_eff are obtained by solving the nonlinear system (26) with constraints (12) and (17) against the same heterogeneous problem that is later used to compute the errors in Eqs. (33)-(34), whereas the baseline diffusion-only model uses the uncalibrated harmonic average (32). The reported improvement therefore confounds two effects: the presence of the advection term and the addition of a second, solution-calibrated parameter. To support the central claim, the authors should include an out-of-sample test (for example, calibrate the coefficients on one set of initial/boundary conditions or on one sub-interval and then evaluate the error on a different set or sub-interval) and should also compare against a diffusion-only model whose D_eff is calibrated by a one-constraint analogue of the same procedure, so that the comparison is between models with the same number of fitted parameters.","section":"Section IV, Eqs. (26)-(28), Table II"},{"comment":"The effective coefficients depend on the initial condition f(x), not only on the diffusivity and boundary conditions as stated in the abstract and in Section V. This is evident because w(x) solves Eq. (20) with source r(x)=f(x)-s(x), and W(x) solves Eq. (23) with R(x)=f(x)-S(x); hence D_eff and v_eff obtained from (26) depend on f through these equations. Homogenized coefficients are normally intrinsic properties of the medium, and a dependence on the initial data is a qualitatively different claim. The authors should either clarify that the coefficients are intended to be problem-specific calibration parameters and discuss the implications, or modify the method so that the coefficients are determined without the initial condition.","section":"Section III, Eqs. (20)-(22) and (23)-(25)"},{"comment":"The nonlinear system (26) is central to the method, but the paper gives no existence or uniqueness results for its solution, no discussion of how the initial guess for fsolve is chosen, and no sensitivity analysis. The two constraints (12) and (17) are introduced as a modeling choice without a derivation from homogenization theory and without error bounds. In addition, no grid-convergence study is reported for N_x, N_t, or the fsolve tolerances, so the errors in Table II cannot be verified as converged. Please add the numerical details and a convergence study, and comment on the well-posedness of the system (26).","section":"Section III, Eq. (26) and Section IV"}],"minor_comments":[{"comment":"The operator defined in the text appears as Lϕ := D_eff ∂²ϕ/∂x − v_eff ∂ϕ/∂x; the second derivative should be ∂²ϕ/∂x².","section":"Section III, paragraph before Eq. (23)"},{"comment":"The word 'homgenization' in the conclusions should be 'homogenization'.","section":"Section V"},{"comment":"The phrase 'To obtain the the effective coefficients' contains a duplicated 'the'.","section":"Section IV, first paragraph"},{"comment":"In the formula for D(x_i), the expression exp(20(x−0.5)) uses an unbound variable x; it should be exp(20(x_i−0.5)) or the formula should otherwise be clarified.","section":"Table I, Case F"},{"comment":"The initial guess supplied to fsolve for the pair (D_eff, v_eff) is not reported; this information is needed for reproducibility and to assess the possibility of multiple solutions of (26).","section":"Section IV, numerical details"}],"recommendation":"major_revision","confidential_remarks":"The central idea is interesting and likely within the scope of physics.comp-ph, and the code availability is commendable. The main obstacle is the in-sample nature of the numerical evidence: the effective coefficients are fitted to each problem before the errors are computed. I would ask the authors to add an out-of-sample transfer test and a parameter-count-controlled comparison, and to address the initial-condition dependence of the coefficients. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read Carr's paper. Bottom line: it is a clean, honest numerical method with a genuine idea, but the headline claim is not yet proven because the effective coefficients are calibrated to the same problems used for evaluation. I would send it to peer review, with the expectation that the revision adds out-of-sample tests and a fair two-parameter baseline.\n\nWhat is new is the specific way of choosing D_eff and v_eff: enforce equality of the spatial average of the steady state and of the time-integrated deviation w(x). That constraint pair is not in the earlier stochastic paper [16]. The derivations are clear, the boundary value problems for w(x) and W(x) are neat, and the code is on GitHub, so the numbers are reproducible. The nine test cases cover smooth and rough diffusivities, periodic and random, time-dependent boundary conditions, and non-uniform initial data. In cases with large D'(x), the advection-diffusion model clearly beats the harmonic-average diffusion model, and v_eff has the physically sensible sign.\n\nThe soft spots are real but not fatal. The main one is in-sample evaluation. For each test case, D_eff and v_eff are computed by solving constraints that involve the same boundary and initial conditions that are later used to compute the error. The harmonic-average Deff is not fitted at all. So Table II conflates two things: the benefit of the advection term and the benefit of having a second, problem-calibrated parameter. The paper does not show that a diffusion-only model with a fitted Deff performs worse. The effective coefficients also depend on the initial condition, through r(x) in (20), which is not highlighted. There are no uniqueness proofs for the nonlinear system (26) and no numerical convergence study. These are addressable, and the method itself is not broken.\n\nWho is this for? Applied mathematicians and engineers doing one-dimensional homogenization or model reduction for diffusion in heterogeneous media, especially with strong spatial gradients in D. It is a useful contribution but not a landmark. I would accept it for peer review and ask for out-of-sample tests, a like-for-like baseline, and a short convergence study.","headline":"A clean, reproducible method that fits an advection-diffusion homogenized model to two moment constraints, but the reported gains are partly calibration and deserve out-of-sample tests.","tokens_in":9024,"tokens_out":3587,"would_cite":false,"duration_ms":37405,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B27","35K57","65M08"],"pacs":[],"model":"deepseek-v4-flash","headline":"An effective advection term in the homogenized equation makes constant-coefficient approximations of heterogeneous-medium diffusion more accurate than the standard harmonic-average model, across nine test cases.","keywords":["homogenization","effective diffusivity","effective velocity","advection-diffusion equation","moment matching","heterogeneous diffusion","boundary value problems"],"falsifier":"Construct a one-dimensional problem with the paper's Dirichlet data but an initial condition concentrated on a small interval away from both boundaries, say a sharply peaked Gaussian centred at x = L/3, and a diffusivity whose derivative changes sign; compute the effective coefficients from constraints (12) and (17), then compare the full time-dependent homogenized solution with the heterogeneous solution. If the homogenized solution is farther from the heterogeneous solution than the harmonic-average diffusion-only model in any region or time window, or if large local errors coexist with matched spatial averages, the blanket improvement claim fails for that case.","tokens_in":8034,"feed_emoji":"","tokens_out":11736,"duration_ms":120782,"temperature":0.7,"pith_summary":"The paper proposes an alternative homogenized replacement for one-dimensional diffusion with spatially varying diffusivity: instead of a pure diffusion equation with a single effective diffusivity, use an advection-diffusion equation with an effective diffusivity and an effective velocity. The two effective coefficients are fixed by matching the spatial average of the steady-state solution and the spatial average of the time-integrated deviation from steady state between the heterogeneous and homogenized models. The paper reports that across nine test cases this advection-diffusion model has a smaller mean absolute error than the standard harmonic-average homogenized model in every case. If the claim holds, the method offers a low-cost way to model transport in layered or heterogeneous media where the standard effective-diffusivity description is visibly wrong.","feed_headline":"Advection improves homogenized diffusion in heterogeneous media","feed_subtitle":"Matching two averages yields an effective diffusivity and velocity, beating the standard model in all nine tests.","key_machinery":"The engine is a pair of matching constraints applied to two boundary value problems. The first is the steady-state profile $s(x)$ of the heterogeneous model; the second is $w(x)=\\int_0^\\infty [s(x)-u(x,t)]\\,dt$, the signed, time-accumulated gap between the transient solution and steady state. The paper derives boundary value problems for $w(x)$ and its homogenized analogue $W(x)$, so both quantities are obtained without solving the time-dependent problem. Equating the spatial averages of $s$ and $S$ and of $w$ and $W$ gives two nonlinear equations in the unknown constants $D_{\\rm eff}$ and $v_{\\rm eff}$; solving that small system supplies the homogenized advection-diffusion model. In the classical rapidly oscillating limit the extra equation returns a near-zero velocity, so the model collapses toward the standard harmonic-average description.","core_discovery":"Starting from $\\partial u/\\partial t = D(x)\\partial^2 u/\\partial x^2 + D'(x)\\partial u/\\partial x$, the paper's central claim is that a constant-coefficient homogenized equation should keep the $\\partial u/\\partial x$ term: $\\partial U/\\partial t = D_{\\rm eff}\\partial^2 U/\\partial x^2 - v_{\\rm eff}\\partial U/\\partial x$. The coefficients are chosen by solving two uncoupled boundary value problems over the heterogeneous medium and enforcing that the spatial averages of the steady state and of the time-integrated transient deviation agree with the homogenized model's counterparts. In all nine test cases, covering smooth, oscillatory, random, piecewise-constant, and piecewise-linear diffusivities plus time-dependent boundary data and non-uniform initial data, the advection-diffusion model has a smaller mean absolute error than the standard diffusion-only model built on the harmonic average of $D(x)$. The effective coefficients also inherit a dependence on the boundary conditions, not just on the diffusivity.","pith_inferences":["Because the two constraints match only global spatial averages, a sharply localised initial condition or a localised sink could preserve the matched averages while producing large early-time local errors; this is a testable stress case the paper does not run.","Since $w(x)$ and $W(x)$ are zeroth temporal moments, imposing constraints on first or second temporal moments is a natural route to fixing more coefficients in two or three dimensions; the paper flags this direction in its conclusion.","The boundary-condition dependence of $v_{\\rm eff}$ suggests the effective velocity is better read as a representation of boundary-driven asymmetry in the homogenized equation than as a material property of the medium; that interpretative framing is not the paper's own."],"forward_implications":["For media with a strong spatial trend in $D(x)$, the advection-diffusion homogenized model reduces mean absolute error by a factor of several relative to the harmonic-average diffusion-only model, as in the paper's cases B, D, E, and F.","Effective coefficients are boundary-condition dependent, so a tabulated or measured effective diffusivity should not be reused for a different boundary setting without recomputing the coefficients.","The coefficient calculation requires only two steady-state boundary value problems plus a small nonlinear solve, so it stays cheaper than resolving the full heterogeneous transient.","The model accommodates time-dependent boundary data and non-uniform initial data while remaining a constant-coefficient equation, which makes it practical for heat and mass transfer applications.","In higher dimensions, additional constraints of the same type would be needed for the extra diffusivity and velocity components, the direction the paper itself points to."],"supporting_citations":[{"why":"supplies the stochastic-diffusion homogenization approach that motivates writing the homogenized equation as advection-diffusion.","marker":"[16]"},{"why":"defines the harmonic-mean effective diffusivity used as the standard diffusion-only baseline in the comparisons.","marker":"[12]"},{"why":"provides the technique that converts the time-integrated deviation w(x) into a boundary value problem without solving the transient model.","marker":"[17]"},{"why":"contributes the response-time and advection-diffusion context behind the temporal-deviation constraint.","marker":"[19]"},{"why":"supplies the asymptotic-expansion and volume-averaging background that justifies the classical harmonic-average model the new model is compared against.","marker":"[7]"}],"fun_headline_variants":["Advection term sharpens homogenized diffusion models","Adding drift improves effective diffusion equations","Advection-diffusion wins in nine diffusion tests","Two-average fit beats harmonic-mean diffusion model","Boundary-aware advection boosts diffusion homogenization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method assumes that matching two numbers, the spatial average of the steady-state profile and the spatial average of the time-integrated gap between transient and steady state, fixes the two effective coefficients and makes the homogenized solution accurate everywhere and at every time, an assumption the paper adopts rather than derives from homogenization theory.","fun_headline_variants_meta":{"raw":{"variants":["Advection term sharpens homogenized diffusion models","Adding drift improves effective diffusion equations","Advection-diffusion wins in nine diffusion tests","Two-average fit beats harmonic-mean diffusion model","Boundary-aware advection boosts diffusion homogenization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000117,"raw_usage":{"total_tokens":1045,"prompt_tokens":880,"completion_tokens":165,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":96}},"tokens_in":496,"tokens_out":165,"duration_ms":2832,"temperature":1.0,"reasoning_tokens":96,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:44:27.868871+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Construct a one-dimensional problem with the paper's Dirichlet data but an initial condition concentrated on a small interval away from both boundaries, say a sharply peaked Gaussian centred at x = L/3, and a diffusivity whose derivative changes sign; compute the effective coefficients from constraints (12) and (17), then compare the full time-dependent homogenized solution with the heterogeneous solution. If the homogenized solution is farther from the heterogeneous solution than the harmonic-average diffusion-only model in any region or time window, or if large local errors coexist with matched spatial averages, the blanket improvement claim fails for that case.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the stochastic-diffusion homogenization approach that motivates writing the homogenized equation as advection-diffusion."},{"cited_title":"Hornung, Homogenization and Porous Media, Springer-Verlag, New York, 1997","cited_arxiv_id":null,"evidence_quote":"defines the harmonic-mean effective diffusivity used as the standard diffusion-only baseline in the comparisons."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the technique that converts the time-integrated deviation w(x) into a boundary value problem without solving the transient model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"contributes the response-time and advection-diffusion context behind the temporal-deviation constraint."},{"cited_title":"Davit, C","cited_arxiv_id":null,"evidence_quote":"supplies the asymptotic-expansion and volume-averaging background that justifies the classical harmonic-average model the new model is compared against."}],"review_version":1}