{"id":"d7f38a27-081b-48ba-a193-2005753d40f1","arxiv_id":"2508.08090","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Global weak solutions exist for the quasi-incompressible Navier-Stokes/Cahn-Hilliard model with fractional diffusion, and as the density mismatch alpha tends to zero the solutions converge to Model H at rate alpha on the lifespan of a strong limit solution.","lead":"This paper proves that a diffuse-interface model for two fluids with slightly different densities has global weak solutions, and that as the density difference shrinks to zero the model converges to the standard matched-density Model H. The interest for PDE analysts is that this is the first rigorous incompressible-limit statement for the quasi-incompressible model in which pressure enters the chemical potential.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 4.5's improved regularity (1.14) is not established: the step from (4.7) to (4.13) replaces Λ^{s+γ/2}(ψϕδ) by ψΛ^{s+γ/2}ϕδ without estimating the commutator, an estimate needed for the α-uniform pressure controls in the incompressible limit.","rationale":"The paper's main advertised novelty is the α-uniform incompressible limit; this rests on the improved regularity (1.14) and the pressure controls of §5.2. I focused the stress test on Lemma 4.5 because it is the only place where the proof appears to skip a nontrivial analytic step: the commutation of Λ^{s+γ/2} with a time cutoff ψ. The reader is right that Theorem 5.1 is conditional on the existence of a strong solution of Model H with regularity (5.5); this is an honest limitation and standard for relative-entropy arguments, but it is not a flaw in the derivation. The commutator gap is different: it is internal and, if it cannot be filled, the uniform bound (1.14) and the whole limit argument would be unproved. I am not claiming the step is wrong—only that it is not justified in the text. The suggested test—expanding (4.7) to (4.13) with an explicit commutator bound—settles the issue. This does not change the reader's CONDITIONAL verdict; it reinforces it. The existence proof for fixed α (Theorem 1.3) is standard (implicit time discretization, compactness), and the relative-entropy structure is coherent. No circularity or fitted parameters are present. If the authors supply the missing estimate, I would regard the central claim as supported.","tokens_in":37529,"tokens_out":14016,"duration_ms":140744,"concrete_test":"Independently rewrite the step (4.7)→(4.13) with the commutator written out. Concretely, compute (Λ^{s+γ/2}ϕδ, Λ^{s+γ/2}(ψϕδ)) = (Λ^{s+γ/2}ϕδ, ψΛ^{s+γ/2}ϕδ) + (Λ^{s+γ/2}ϕδ, [Λ^{s+γ/2},ψ]ϕδ) and prove a bound of the form |(Λ^{s+γ/2}ϕδ, [Λ^{s+γ/2},ψ]ϕδ)| ≤ C(∥ψ∥_{L∞} + α∥∂tψ∥_{L1} + (1/4)∥ψ^{1/2}Λ^{s+γ/2}ϕδ∥^2) with C independent of α (e.g. via Kenig–Ponce–Vega type commutator estimates). If the only available bound contains an explicit 1/α factor, then (4.13) is false as stated, and Lemma 4.5, together with the α-uniform pressure controls in §5.2, lack justification. Also verify the cancellation in (4.9) by tracking the 1/α factors explicitly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 5.1 relies on the improved order-parameter bound (1.14) obtained in Lemma 4.5; this bound is used in Lemmas 5.3–5.4 to control the pressure terms pα f and pα∇ϕα·g uniformly in α and in (5.20)/(5.29) to pass the αpα terms to zero. In the proof of Lemma 4.5, the displayed identity (4.7) has on the left (ζα^{-1}Λ^{s+γ/2}ϕδ, Λ^{s+γ/2}(ψϕδ))_QT. The conclusion (4.13) is (Λ^{s+γ/2}ϕδ, ψΛ^{s+γ/2}ϕδ)_QT ≤ C(α∥∂tψ∥_{L1}+∥ψ∥_{L∞}). The passage replaces Λ^{s+γ/2}(ψϕδ) by ψΛ^{s+γ/2}ϕδ, i.e. discards the commutator [Λ^{s+γ/2},ψ]ϕδ = Λ^{s+γ/2}(ψϕδ)-ψΛ^{s+γ/2}ϕδ, without any estimate. Since s+γ/2 > 3/2, this is a fractional differential operator of positive order; the bilinear form (Λ^{s+γ/2}ϕδ,[Λ^{s+γ/2},ψ]ϕδ) is not controlled by the right-side terms (4.9)–(4.12) in any obvious way. In particular, J1 is bounded in (4.9) by a term of size (1/α)∥ψ∥_{L∞}; the cancellation leading to the α-independent bound in (4.13) is not shown. If the commutator term contributes even a constant multiple of (1/α)∥Λ^{s+γ/2}ϕδ∥^2, the claimed uniform bound fails, and Lemma 5.4's use of H^{s+1/2}(ϕα) in (5.28) is unjustified. This is an internal gap in the proof of a lemma advertised as the key novelty; the manuscript should either supply the commutator estimate or state explicitly which standard fractional calculus result is being used.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the quasi-incompressible Navier–Stokes/Cahn–Hilliard system (1.1) in a three-dimensional periodic domain, with unmatched densities and a fractional Laplacian in the chemical potential. It first proves global existence of weak solutions (Theorem 1.3) by an implicit time-discretization scheme and a Leray–Schauder fixed-point argument, and claims an improved order-parameter regularity, ϕ ∈ L^2(0,T;H^{s+γ/2}), stated in (1.14) as the key novelty. It then uses the relative entropy method to prove an incompressible limit (Theorem 5.1) as the density difference α → 0: weak solutions of the quasi-incompressible system converge to a strong solution of Model H, with relative energy rate α and L^2 rate √α, under a well-prepared data condition. The proof relies on non-standard uniform-in-α pressure controls (Lemmas 5.3–5.5) that exploit the claimed improved regularity of the order parameter.","tokens_in":37953,"tokens_out":6335,"duration_ms":73900,"significance":"If the proof can be completed, the paper would make a substantial contribution: it provides the first global weak-solution existence for this fractional quasi-incompressible two-phase model, introduces a genuinely new partial-damping regularity mechanism for the order parameter, and gives the first rigorous incompressible limit for the mass-averaged velocity formulation, with explicit convergence rates. The approximation and compactness architecture is standard but carefully executed, and the relative-entropy estimates in Section 5 are detailed. The paper is also honest about the conditional nature of the convergence statement: it assumes the existence of a strong solution of Model H with the regularity (5.5). However, the central improved-regularity lemma contains an unestimated commutator passage, and since that lemma is load-bearing for the pressure controls and the final convergence theorem, the main results are not fully established as written.","major_comments":[{"comment":"The step from (4.7) to (4.13) is not justified. The left-hand side of (4.7) contains Λ^{s+γ/2}(ψϕδ), while (4.13) is written with ψΛ^{s+γ/2}ϕδ. Replacing one by the other requires estimating the commutator [Λ^{s+γ/2}, ψ]ϕδ, and the manuscript gives no estimate for this term. Since s + γ/2 > 3/2, this is a positive-order fractional differential operator; the bilinear form (Λ^{s+γ/2}ϕδ, [Λ^{s+γ/2}, ψ]ϕδ) is not controlled by the bounds (4.9)–(4.12), particularly because J1 is only bounded with a factor 1/α. Without this commutator estimate, the α-independent bound (1.14) is not established. This is load-bearing: Lemma 5.4 uses ∥ϕα∥_{L^2(0,T;H^{s+1/2})} in (5.28), and Lemma 5.3/(5.20) uses the same improved integrability to pass α pα terms to zero. The authors should either supply a valid commutator estimate (with the precise function-space assumptions) or state explicitly which standard fr","section":"§4, Lemma 4.5 (Eqs. (4.7)–(4.13))"},{"comment":"The convergence theorem is conditional on the existence of a strong solution of Model H satisfying the high regularity (5.5): u ∈ H^1(0,T';H^{s+1/2}), p ∈ L^2(0,T';H^1), μ ∈ H^1(Q_{T'}), and ϕ ∈ H^2(0,T';H^s). The paper does not verify that the known local strong solutions of Model H in 3D (cf. [1,36]) attain these regularity levels. If such solutions are known only on a short interval, or only with lower regularity, then the theorem's hypothesis may be empty outside a very restrictive class. The authors should either prove or cite a local well-posedness result that yields exactly (5.5), or explicitly restate Theorem 1.6/5.1 as a conditional statement with this regularity class as part of the hypothesis. As written, the applicability of the advertised incompressible limit is not demonstrated.","section":"§5.1, Theorem 5.1 (assumption (5.5))"}],"minor_comments":[{"comment":"There are numerous typographical errors: 'pinciple part' (§3), 'F ormal argument' (§5.1), 'quasi-compressible' (§1.3), 'imcompressible' (Remark 1.7), 'Date avability' (Compliance statement), and inconsistent punctuation in displayed equations. A careful copyedit is needed.","section":"Throughout"},{"comment":"The compactness passage in the time-discretization limit uses an Aubin–Lions argument with H^1(T^3) ↪↪ H^s(T^3) written for '0 ≤ s < 1', but the symbol s already denotes the order of the fractional Laplacian with s > 3/2. This could confuse the reader; the compact embedding should be written with a different symbol (e.g., H^1 ↪↪ H^r, r < 1).","section":"§3, proof of Lemma 3.3"},{"comment":"The bound (5.28) is a fractional Leibniz-type estimate for Λ^{s-1/2}(∇ϕα·g), but no reference or proof is given. Such product estimates are nontrivial in this Sobolev range and should be stated explicitly, together with the required regularity of g.","section":"§5.2, Lemma 5.4 (Eq. (5.28))"},{"comment":"In the H_i estimates after (5.40), many constants are aggregated into C(T',D). It would improve readability to state explicitly which terms are absorbed by the dissipation and which are handled by the Gronwall term, e.g., by numbering the final estimates (5.43)–(5.54) with a table.","section":"§5.3, proof of Theorem 5.1"}],"recommendation":"major_revision","confidential_remarks":"The central gap in Lemma 4.5 is serious but may be fixable; if the authors can supply the missing commutator estimate, the paper would be a solid contribution. The conditional strong-solution regularity in (5.5) should also be checked against the existing literature, since it affects the scope of the advertised convergence result. I recommend major revision rather than rejection because the overall architecture is coherent and the gap is localized."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick take: this paper is worth taking seriously. It does two things: global weak solutions for a fractional quasi-incompressible Navier–Stokes/Cahn–Hilliard system with unmatched densities and pressure in the chemical potential, and then an incompressible limit to Model H via relative entropy, with rate alpha. The limit result is genuinely new—Abels explicitly left this direction open because of pressure regularity, and the paper gets around it with the improved order-parameter bound (1.14) and non-standard pressure controls in Lemmas 5.3–5.4.\n\nWhat the paper does well: the proof architecture is coherent and standard for this area (time discretization, fixed point, compactness, then relative entropy). The energy estimates are laid out carefully. The citation pattern is honest—they distinguish the mixture-theory model [28], the volume-averaged AGG model [34], and the compressible cases [10,31], and they identify the specific gap in Abels' work. No fitted parameters, no circular reasoning. This is a real step forward.\n\nSoft spots, in proportion: the stress-test concern about Lemma 4.5 does not land. In that lemma, ψ is a cutoff depending only on time, not space. So Λ^{s+γ/2}(ψφδ) = ψΛ^{s+γ/2}φδ exactly; there is no commutator to estimate. The step from (4.7) to (4.13) is legitimate. The genuine caveat is the one the reader flagged: Theorem 5.1 is conditional on the existence of a strong solution of the 3D limit Model H with the high regularity (5.5), and the authors simply assume this life span rather than checking it against known local well-posedness results. That makes the theorem conditional by design, and the abstract slightly oversells by not mentioning the strong-solution lifespan. It is a stated caveat, not a hidden flaw, but it should be fixed—either verify (5.5) against existing results or downgrade the abstract's wording. Minor: the proof would be easier to trust with a few more details in the discrete-to-continuous passage, but nothing there looks broken.\n\nBottom line: the math is coherent, the new result is real, and the conditional nature of the limit theorem is not a fatal objection. Definitely deserves a serious referee; I'd send it out. I would cite it if I worked on diffuse interface models, and I'd put it on the reading group list.\n\nRecommendation: engage with the paper, but push the authors to clarify the strong-solution assumption and align the abstract with the theorem.","headline":"Solid analysis: first rigorous incompressible limit for the mass-averaged quasi-incompressible model, with the main caveat being the conditional strong-solution regularity on the limit Model H, not the commutator step the stress-test flagged.","tokens_in":38598,"tokens_out":3710,"would_cite":true,"duration_ms":38538,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q35","76T06","76T99","35D30","35B25","35Q30"],"pacs":[],"model":"deepseek-v4-flash","headline":"A density-mismatched two-phase flow model is shown to have global weak solutions and to converge to the classical incompressible Model H as the density difference vanishes.","keywords":["Navier-Stokes/Cahn-Hilliard","quasi-incompressible","two-phase flows","weak solutions","relative entropy method","incompressible limit","fractional Laplacian","unmatched densities"],"falsifier":"Take a family of well-prepared initial data for which the 3D Model H strong solution blows up at a finite time T*; on any interval [0,T'] with T' > T*, Theorem 5.1 has no content, so the claimed convergence cannot be observed there. A sharper falsifier: numerically compute the left side of (5.8) for small α in a smooth test case; if the relative energy decays slower than linearly in α (e.g., like √α), then the integral estimate is not optimal, and if it grows, the theorem's bound is violated.","tokens_in":37244,"feed_emoji":"🌊","tokens_out":9303,"duration_ms":89562,"temperature":0.7,"pith_summary":"This paper establishes two results for a quasi-incompressible Navier–Stokes/Cahn–Hilliard system that models two viscous fluids of different densities in a three-dimensional periodic box. First, it proves the existence of global weak solutions, upgrading the order parameter's regularity to $\\phi \\in L^2(0,T;H^{s+\\gamma/2})$ with $\\phi \\in (-1-\\theta,1+\\theta)$ by using a partial damping effect of the capillary force. Second, via the relative entropy method, it proves that as the density-difference parameter $\\alpha \\to 0$, these weak solutions converge to the unique strong solution of the classical incompressible Model H on the strong solution's lifespan, with the squared relative energy and the $L^2$ distance bounded by $C(T',D)(\\alpha + \\text{initial relative energy})$. The significance is that this is the first rigorous incompressible-limit passage for this quasi-incompressible two-phase model, overcoming the missing uniform pressure bounds with non-standard pressure estimates.","feed_headline":"Weak solutions exist and incompressible limit holds for two-phase model","feed_subtitle":"As the density difference between the fluids shrinks, solutions converge to the classical Model H at rate √α.","key_machinery":"The proof stands on two mechanisms. Existence is built through an implicit time discretization of a twice-regularized system (parameters $\\delta$ and $\\alpha$), solved by a fixed-point argument, followed by a compactness passage; the novel step is a regularity estimate for the order parameter obtained by testing the momentum equation against $\\nabla \\Delta^{-1}\\Lambda^\\gamma(\\psi\\phi_\\delta)$, which turns the capillary force $\\phi\\nabla\\mu$ into damping of $\\Lambda^{s+\\gamma/2}\\phi$. The incompressible limit is driven by the relative entropy functional (5.10), whose coercivity controls $\\|u-u_\\alpha\\|_{L^2}^2 + \\|\\phi-\\phi_\\alpha\\|_{H^s}^2$; the remainder terms, including those involving the","core_discovery":"The paper's central claim is that the quasi-incompressible model (1.1)—where the mass-averaged velocity is not divergence-free, $\\operatorname{div}u = \\alpha \\Delta \\mu_p$, and the pressure enters the chemical potential as $\\mu_p = \\mu + \\alpha p$—is globally well-posed in the weak sense and converges to Model H in the limit of vanishing density contrast. Theorem 1.3 asserts global weak solutions in $\\mathbb{T}^3$ for arbitrary finite time, together with the improved order-parameter regularity (1.14). Theorem 5.1 then quantifies the incompressible limit: for well-prepared initial data (5.6)–(5.7) and as long as a strong solution of Model H with regularity (5.5) exists on $[0,T']$, the relati","pith_inferences":["Beyond the paper: the same relative-entropy inequality should yield a weak-strong uniqueness statement for the quasi-incompressible model itself (a weak solution and a strong solution with the same data), since the proof of Theorem 5.1 does not use the specific form of the limit system except through its regularity.","Beyond the paper: the uniform-in-$\\alpha$ pressure controls of Lemmas 5.3–5.5 resemble effective-flux-type estimates used in compressible fluid limits; they may transfer to the low-Mach-number limit of this model, connecting with the compressible diffuse-interface results the paper cites.","Beyond the paper: the confinement $\\phi\\in(-1-\\theta,1+\\theta)$ is not strict separation; a testable extension is whether the capillary damping estimate can be pushed to force $\\phi\\in(-1,1)$ when the free energy is chosen with a logarithmic singularity, which the paper identifies as desirable.","Beyond the paper: numerical experiments on smooth test cases could probe whether the $L^2$ error decays like $\\sqrt{\\alpha}$ or linearly; a faster empirical rate would indicate the abstract estimate is not sharp."],"forward_implications":["If Theorem 1.3 is correct, the quasi-incompressible model with fractional Laplacian and unmatched densities admits global weak solutions in 3D for every finite time, with the order parameter confined to a bounded neighborhood of $[-1,1]$.","Theorem 5.1 gives a quantitative justification for replacing the quasi-incompressible model by the simpler Model H when densities are nearly matched: the error in $L^2(0,T';H^1)$ for the velocity and $L^2$ for the chemical-potential gradient is $O(\\sqrt{\\alpha})$ once the initial data are well prepared.","The improved regularity (1.14) is the mechanism that makes the pressure-independent estimates possible; without it, the paper argues, the incompressible limit is out of reach under the stated assumptions.","The relative entropy inequality (5.56) implies a weak-strong uniqueness principle: any weak solution built by Theorem 1.3 coincides with the strong solution as long as the latter exists."],"supporting_citations":[{"why":"Supplies the existence strategy for weak solutions of diffuse-interface models with general densities, which the paper generalizes to the fractional-Laplacian setting.","marker":"[2]"},{"why":"Provides strong well-posedness results for the quasi-incompressible model that the present existence proof builds on.","marker":"[3]"},{"why":"Gives the original derivation of the quasi-incompressible Cahn–Hilliard model that (1.1) is a variant of.","marker":"[41]"},{"why":"Re-derives the model (1.1) with s=1 from mixture theory, the starting point for the system studied here.","marker":"[48]"},{"why":"Introduces Model H, the incompressible matched-density limit system to which the paper proves convergence.","marker":"[37]"},{"why":"Classical reference for Model H used as the limit system in the incompressible-limit theorem.","marker":"[39]"},{"why":"Gives existence of local strong solutions of the 3D Model H, on which the regularity assumption (5.5) is based.","marker":"[36]"},{"why":"Applies the relative entropy method to a compressible diffuse-interface model, the technique the paper adapts for the incompressible limit.","marker":"[10]"},{"why":"Studies the incompressible limit of a quasi-incompressible fluid model with the same parameter α, providing the comparison framework.","marker":"[28]"},{"why":"Uses relative energy for a compressible two-phase diffuse-interface model, serving as a methodological basis for the weak-strong argument.","marker":"[31]"}],"fun_headline_variants":["Weak solutions and Model H limit for two-phase flows","Two-phase flow: existence and incompressible limit proven","Density contrast vanishes, two-phase flow converges to Model H","Global weak solutions for quasi-incompressible two-phase flow","Existence and incompressible limit for viscous two-phase model"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The incompressible-limit theorem presupposes, without proof, that the three-dimensional limit system Model H possesses a strong solution on the whole interval [0,T'] with the high regularity listed in (5.5); if such a solution exists only on a shorter interval or not at all, the convergence statement is empty outside that interval.","fun_headline_variants_meta":{"raw":{"variants":["Weak solutions and Model H limit for two-phase flows","Two-phase flow: existence and incompressible limit proven","Density contrast vanishes, two-phase flow converges to Model H","Global weak solutions for quasi-incompressible two-phase flow","Existence and incompressible limit for viscous two-phase model"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001924,"raw_usage":{"total_tokens":7393,"prompt_tokens":790,"completion_tokens":6603,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":6521}},"tokens_in":534,"tokens_out":6603,"duration_ms":51408,"temperature":1.0,"reasoning_tokens":6521,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T21:41:15.354047+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a family of well-prepared initial data for which the 3D Model H strong solution blows up at a finite time T*; on any interval [0,T'] with T' > T*, Theorem 5.1 has no content, so the claimed convergence cannot be observed there. A sharper falsifier: numerically compute the left side of (5.8) for small α in a smooth test case; if the relative energy decays slower than linearly in α (e.g., like √α), then the integral estimate is not optimal, and if it grows, the theorem's bound is violated.","supporting_citations":[{"cited_title":"Abels, Existence of Weak Solutions for a Diffuse Interface Model for Viscous, Incompressible Fluids with General Densities, Commun","cited_arxiv_id":null,"evidence_quote":"Supplies the existence strategy for weak solutions of diffuse-interface models with general densities, which the paper generalizes to the fractional-Laplacian setting."},{"cited_title":"Abels, Strong well-posedness of a diffuse interface model for a viscous, quasi-incompressible two-phase flow , SIAM J","cited_arxiv_id":null,"evidence_quote":"Provides strong well-posedness results for the quasi-incompressible model that the present existence proof builds on."},{"cited_title":"Lowengrub, L","cited_arxiv_id":null,"evidence_quote":"Gives the original derivation of the quasi-incompressible Cahn–Hilliard model that (1.1) is a variant of."},{"cited_title":"Shokrpour Roudbari, G","cited_arxiv_id":null,"evidence_quote":"Re-derives the model (1.1) with s=1 from mixture theory, the starting point for the system studied here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces Model H, the incompressible matched-density limit system to which the paper proves convergence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Classical reference for Model H used as the limit system in the incompressible-limit theorem."},{"cited_title":"Giorgini, A","cited_arxiv_id":null,"evidence_quote":"Gives existence of local strong solutions of the 3D Model H, on which the regularity assumption (5.5) is based."},{"cited_title":"Abels, Y","cited_arxiv_id":null,"evidence_quote":"Applies the relative entropy method to a compressible diffuse-interface model, the technique the paper adapts for the incompressible limit."},{"cited_title":"Feireisl, Y","cited_arxiv_id":null,"evidence_quote":"Studies the incompressible limit of a quasi-incompressible fluid model with the same parameter α, providing the comparison framework."},{"cited_title":"Feireisl, M","cited_arxiv_id":null,"evidence_quote":"Uses relative energy for a compressible two-phase diffuse-interface model, serving as a methodological basis for the weak-strong argument."}],"review_version":1}