{"id":"bda5ad1a-c6a6-4e7a-ab56-2caee22ec7b3","arxiv_id":"2506.11374","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A thermodynamically consistent phase-field model (NSCHAC) with an energy-stable LDG-SDC numerical scheme is developed and tested for osmotic flow through semi-permeable membranes.","lead":"This paper builds a new mathematical model and numerical method for water flowing across semi-permeable membranes under osmotic pressure. The model couples fluid motion, solute concentration, and interface motion, and the method is shown to keep energy decreasing over time.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Energy-stability Theorems 3.1 and 3.2 assume condition (16), which the untruncated double-well potential used in the numerics does not satisfy; the unconditional-stability proof therefore does not cover the implemented scheme.","rationale":"After inspecting the full text, the strongest part is the continuous model: Theorem 2.1's energy law is derived carefully and appears correct. The load-bearing weakness is the link from the continuous model to the discrete schemes. The reader's weakest_assumption correctly identifies condition (16) as the pivot: it is explicitly assumed in Theorems 3.1 and 3.2, the paper admits the quartic potential does not satisfy it, and the numerical section uses the quartic potential without reporting S or an L-infinity bound. That mismatch means the main stability theorems are not about the scheme that is actually run. This is not a mere technicality: the stabilization term in the proof is bounded by S - F''(xi)/2, so if phi leaves the interval where F'' <= 2S, the proof's inequality can break. The absence of an L-infinity bound for the scheme (no maximum principle for the Cahn-Hilliard-type equation) makes this a real possibility rather than a stylistic concern. The SDC extension also lacks a stability theorem, but the paper does not state one for it, so the first-order scheme's stability is the more fundamental issue. The concern is addressable: use the truncated potential in the stability analysis and either verify (numerically) that the implemented runs stay inside the truncated region or adapt S to max|F''|; hence CONDITIONAL is appropriate, and the reader's verdict stands unchanged.","tokens_in":26977,"tokens_out":30495,"duration_ms":307979,"concrete_test":"Re-run the energy-stability experiment (Example 4.2) and the droplet test (Example 4.3, K=0.1) recording max_{n,x}|phi^n(x)| and the stabilization parameter S. Then verify whether S >= (1/2) max_{n,x}|F''(phi^n(x))| holds on every time step. If this inequality holds and all phi^n remain in [-1,1] (so F'' <= 2), the stability proof can be legitimately applied with a truncated potential agreeing with F on the attained range. If instead max|phi^n| > 1 or S is too small, the discrete energy inequality may fail; the theorem's hypothesis is violated for the implemented scheme and the unconditional-stability claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's discrete energy-stability results (Theorems 3.1 and 3.2) are conditional on assumption (16), max|F''(phi)| <= L. For the double-well potential F(phi)=1/4(phi^2-1)^2 used throughout Section 4, F''(phi)=3phi^2-1 is unbounded on R, so no finite L exists. The text acknowledges this and invokes the common practice of a truncated double-well potential (ref [3]), but the numerical scheme (11) and all reported tests use the untruncated F. Therefore the theorems, as stated, apply to a different potential than the one implemented. In the proof, the stabilization term enters as (S - F''(xi)/2)(phi^{n+1}-phi^n)^2; with a fixed S this is nonnegative only if phi^n remains in the region where F'' <= 2S. No L-infinity bound for phi is established for the scheme (the Cahn-Hilliard equation does not preserve a maximum principle generally). Consequently, for initial data that drives |phi| beyond that region, the energy inequality is not guaranteed. The stability parameter S is also never specified in Section 4, so even the bounded-regime hypothesis cannot be verified from the paper. This is a genuine gap in the central claim that the proposed first-order scheme is unconditionally energy stable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Section 2 derives, via the energetic variational approach, a Navier–Stokes–Cahn–Hilliard–Allen–Cahn model (5) for two-phase flow across semi-permeable membranes, with a transmembrane flux proportional to the chemical-potential imbalance through the K|∇φ|μ̃φ term; Theorem 2.1 verifies the continuous energy dissipation law. Section 3 presents a first-order, stabilized temporal scheme (11)–(14), proves discrete energy dissipation under the bounded-second-derivative assumption (16) in Theorem 3.1, extends the result to a fully discrete LDG scheme in Theorem 3.2, and constructs a semi-implicit SDC method for higher temporal accuracy. Section 4 reports convergence orders (Tables 1–2), energy decay (Figure 1), and numerical studies of membrane permeability, osmotic pressure, and shear flow on droplet equilibria (Figures 2–9). The paper's central claims are the thermodynamic consistency of the model, the energy stability of the first-order scheme and its LDG discretization, and the high-order accuracy of the SDC–LDG combination.","tokens_in":27359,"tokens_out":15462,"duration_ms":156684,"significance":"The continuous model is built transparently from an explicit energy and dissipation functional, so Theorem 2.1 is a genuine verification rather than a fitted identity, and the discrete stability proofs for the first-order and LDG schemes are carried out in detail with correct algebraic structure. The manuscript also provides a clear falsifiable prediction—the equilibrium droplet shape depends on solute content but not on permeability K—and the numerical experiments support it. If the gap described in the major comments (assumption (16) versus the implemented untruncated potential) is closed, the paper would be a sound technical contribution to phase-field modeling of osmotic flow, with value for both the model and the numerical framework. The significance is moderate: the ideas are largely assembled from well-established building blocks (stabilized semi-implicit treatment, alternating-flux LDG, SDC), but their combination for this NSCHAC system is new.","major_comments":[{"comment":"The unconditional stability claim (restated in Section 5) is not established for the implemented scheme. Theorems 3.1 and 3.2 assume condition (16), max|F''(φ)| ≤ L, but the temporal scheme (11) and all tests in Section 4 use F(φ)=1/4(φ²−1)², whose second derivative 3φ²−1 is unbounded on ℝ, so no finite L exists. The paragraph after (16) acknowledges this and cites the common practice of truncation, yet the manuscript never states that a truncated potential is actually used in the experiments, and no L∞ bound on the discrete φ that would keep F''(φ) ≤ 2S is established; the discrete energy laws (27) and (53) require (S − f'(ξ^n)/2)(φ^{n+1}−φ^n)² ≥ 0, which forces f'(ξ^n) ≤ 2S at every step. Since the Cahn-Hilliard equation admits no general maximum principle, the proofs as written apply only to a potential different from the one implemented. The authors should either (a) implement and explicitly define the truncated potential, giving the resulting L and S, or (b) prove or verify that the computed solutions stay in a region where (16) holds, and in either case state the stabilization parameter S used in Section 4.","section":"3.1, condition (16); Theorems 3.1, 3.2; scheme (11); Section 4"},{"comment":"The shear-flow experiment imposes velocity boundary data u|y=0=(−1,0)^T and u|y=6.4=(1,0)^T, whereas the model (5) and the stability theorems (Theorem 2.1 and Theorem 3.2) assume the homogeneous condition u|∂Ω=0 from (4); the LDG flux convention (33) specifies bun+1_p=0 and inherits the remaining boundary fluxes from the interior, so it is not stated how non-homogeneous Dirichlet velocity data is enforced. As a result, the energy-stability guarantee does not cover Example 4.4, and the text should say so explicitly, along with a description of how the moving-wall boundary conditions are implemented in the LDG scheme.","section":"4.4, Example 4.4; boundary conditions (4)"},{"comment":"The stabilization parameter S, which the theorems require to satisfy S ≥ L/2, is never reported in Section 4, so the stability hypothesis cannot be checked for the computations of Tables 1–2 and Figures 1–9; the value of S (and, if a truncated potential is used, its definition) must be given for each experiment.","section":"4, Examples 4.1–4.4"}],"minor_comments":[{"comment":"In the derivation of I2 + I3, the terms 'ζ− ln C+ c∞ ∂C−/∂t' contain typos (C+ should be C− in the ζ− terms); please correct the logarithms and make the notation consistent with the objective functional (2).","section":"A, Appendix"},{"comment":"In Table 1, the C+ row for m=3 reports an L2 error of 5.82E-04, which is larger than the m=2 value of 1.15E-04, yet the order column lists 0.98; this looks like a typo (probably 5.82E-05) and should be corrected, since the current entry contradicts the claimed first-order convergence.","section":"Table 1"},{"comment":"The quotient (ζ^{n+1}_±−ζ^n_±)/(φ^{n+1}−φ^n) is undefined when φ^{n+1}=φ^n; the paper should state the convention used (e.g., the limiting value of the secant derivative).","section":"3.1, scheme (11)"},{"comment":"The scheme (11)–(14) is called 'decoupled,' but Step 1 is a coupled system for φ^{n+1} and C^{n+1}_± through the ζ-derivative term and through u^n_*; a sentence clarifying that the decoupling is with respect to the pressure-velocity system (Steps 2–3) would be helpful.","section":"3.1"},{"comment":"The abstract and conclusion describe 'high-order, energy-stable numerical schemes' as a whole, but the energy-stability proofs (Theorems 3.1, 3.2) apply only to the first-order scheme and its LDG version; no stability proof is given for the SDC scheme, whose energy decay is shown only numerically in Figure 1(b). Please rephrase to avoid implying a stability proof for the SDC method.","section":"Abstract and Section 5"},{"comment":"The stated equilibrium conditions 'u = 0, μ̃φ = 0, μ±c = 0' do not follow from the dissipation law of Theorem 2.1: the last dissipation term vanishes already when ∇μ±c = 0 on each phase subdomain, and the equilibrium concentrations are fixed by total solute mass; the simulations' equilibrium concentrations are not close to C± = 1 implied by μ±c = ln C± = 0. The text should state the gradient condition rather than μ±c = 0.","section":"4.3"},{"comment":"In the dimensionless equations (5), the diffusivities D± appear as bare ratios, but their normalization is not defined after the scaling; please define the dimensionless D± (e.g., D±/D∗).","section":"2, dimensionless form"}],"recommendation":"major_revision","confidential_remarks":"This is a numerical-analysis paper submitted to a physics journal; while the topic is fluid dynamics, the journal fit is for the editors to judge. The self-citations [15–17] are directly relevant to the LDG and SDC methodology and do not appear to inflate the contribution. The main revision request is narrow: make the stability hypothesis match the implemented potential and report S; I believe this is feasible and would then make the discrete stability results solid."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe short version: this paper gives a new phase-field model for osmotic flow across semi-permeable membranes, coupling Navier-Stokes, Cahn-Hilliard, and Allen-Cahn dynamics, and it proves an energy-dissipation law for the continuous model plus discrete energy stability for a first-order LDG scheme. The modeling is the strongest part. The derivation from an explicit energy and dissipation functional via the energy variational approach is clean, and the resulting NSCHAC system is a natural extension of prior work. Theorem 2.1 checks out. The LDG spatial discretization is standard but competently handled, and the numerical experiments show the expected convergence orders and monotone energy decay.\n\nThe soft spots are real but addressable. The biggest is the gap between the stability proof and the implementation. Theorems 3.1 and 3.2 require the double-well potential to have bounded second derivative (condition 16), and the paper acknowledges that the standard F does not satisfy it, invoking the common truncated-potential trick. But the numerical scheme and all tests use the untruncated F. No L-infinity bound on phi is shown for the scheme, and the stabilization parameter S is never specified in Section 4. So, as stated, the unconditional stability guarantee does not apply to what was actually computed. This is not a fatal flaw—stabilized semi-implicit treatment may work fine in practice—but it needs to be fixed, either by proving a bound, stating S and checking the condition, or by reworking the theorem.\n\nSecond, energy stability is proven only for the first-order scheme. The SDC method is presented as improving temporal accuracy and the numerics show it is stable in practice, but there is no discrete energy-dissipation theorem for the SDC version. That should be stated clearly, not bundled into the stability claim.\n\nThird, reproducibility is weaker than it should be: the accuracy tests say a source term is added but give no formula, and details of the nonlinear solver are missing.\n\nThe citation pattern is fine; self-citations to the authors' own SDC and LDG methods are appropriate. The model is incremental but the specific combination is new. Yes, this deserves a serious referee and likely a conditional accept after the gap is closed.","headline":"A genuinely useful new phase-field model with a clean continuous energy law, but the discrete stability theorem does not cover the potential used in the code, and the SDC extension is claimed without proof.","tokens_in":27850,"tokens_out":1835,"would_cite":false,"duration_ms":19283,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M60","65H10","76T06","92C10"],"pacs":[],"model":"deepseek-v4-flash","headline":"A thermodynamically consistent phase-field model for osmotic flow across semi-permeable membranes, with discretizations that preserve the energy-dissipation law.","keywords":["phase-field model","semi-permeable membrane","osmotic pressure","energy dissipation law","Navier-Stokes-Cahn-Hilliard-Allen-Cahn","local discontinuous Galerkin","spectral deferred correction","double-well potential"],"falsifier":"Run the LDG scheme (30)-(32) with the untruncated potential and initial data whose phase field takes values outside $[-1,1]$, so that $|F''|$ exceeds any chosen $L$, and check whether the discrete energy $E_h^n$ is monotone non-increasing; a single time step where $E_h^{n+1} > E_h^n$ would refute the claimed unconditional stability as stated.","tokens_in":26792,"feed_emoji":"💧","tokens_out":5941,"duration_ms":59686,"temperature":0.7,"pith_summary":"This paper sets out to show that osmotic water transport across a semi-permeable membrane can be described by a phase-field model that is thermodynamically consistent: the total free energy, combining kinetic, entropic, and mixing contributions, can never increase along solutions. The model augments the Navier–Stokes–Cahn–Hilliard system with an Allen–Cahn-type transmembrane flux proportional to the chemical-potential imbalance, so water can cross the diffuse interface while solute cannot. The authors then construct numerical schemes, a first-order decoupled scheme and a spectral-deferred-correction lift to higher order, and prove discrete energy-dissipation laws for both. If correct, the framework gives a computational tool for osmotic effects in biology and industry where the simulated dynamics obey the second law of thermodynamics by construction.","feed_headline":"Osmotic membrane flow gets a proven energy-stable model","feed_subtitle":"Water crosses semi-permeable interfaces while total free energy never increases, in theory and in the discrete scheme.","key_machinery":"The machinery is the free-energy functional $E = E_{\\mathrm{kin}} + E_{\\mathrm{entropy}} + E_{\\mathrm{mix}}$ with chemical potential $\\tilde{\\mu}_\\phi = -\\epsilon\\Delta\\phi + \\frac{1}{\\epsilon}F'(\\phi) - \\frac{\\beta}{\\epsilon}(\\partial_\\phi \\zeta_+ C_+ + \\partial_\\phi \\zeta_- C_-)$, and the Allen-Cahn-type transmembrane flux $S_\\phi = -\\frac{K}{Ca}\\tilde{\\mu}_\\phi |\\nabla\\phi|$ that lets the interface move relative to the fluid. Its discrete support comes from a modified intermediate velocity $u_*$ that decouples the phase-field and concentration updates from the momentum update, and from an alternating-flux choice in the LDG scheme that makes the numerical operators adjoint-consistent so the telescoping energy argument survives discretization. The stabilization term $\\frac{S}{\\epsilon}(\\phi^{n+1}-\\phi^n)$ with $S \\geq L/2$ is what makes the first-order semi-implicit treatment of $F(\\phi)$ unconditionally energy-stable under condition (16).","core_discovery":"On the paper's own terms, the central discovery is that the osmotic coupling can be encoded in a single extra flux term, $S_\\phi = -\\frac{K}{Ca}\\tilde{\\mu}_\\phi |\\nabla\\phi|$, that preserves energy dissipation: Theorem 2.1 proves for the continuous system (5) that the time derivative of the total free energy equals a sum of negative semi-definite dissipation terms. The same dissipation structure is shown to survive time discretization (Theorem 3.1) and full local-discontinuous-Galerkin discretization (Theorem 3.2), provided the stabilizing constant $S$ is at least half the maximum of $|F''|$. The equilibrium condition $u=0$, $\\tilde{\\mu}_\\phi=0$, $\\mu_c^\\pm=0$ then predicts that droplet equilibrium shape is determined by surface tension and solute loading, not by membrane permeability, which the simulations in Section 4 confirm.","pith_inferences":["If the membrane permeability $K$ is made spatially varying or time-dependent, the same energy argument should yield a dissipation law with a time-dependent coefficient, because $K$ enters only the non-negative dissipative term; this extension is not stated in the paper.","The proof relies on replacing $F$ by a truncated potential, so a practical user should monitor $\\max|\\phi|$ over the run; a testable hypothesis is that the unconditionally stable behavior persists as long as $\\phi$ stays in the interval where $F''$ is bounded by the chosen $S$.","The decoupled structure via $u_*$ suggests a natural splitting for adding other physics, such as electrodiffusion or active membrane stresses, because osmotic feedback enters only through the modified advection velocity."],"forward_implications":["At equilibrium the model predicts $u = 0$, $\\tilde{\\mu}_\\phi = 0$, and $\\mu_c^\\pm = 0$, so droplet shape is set by surface tension and total solute mass, not by permeability; higher permeability only accelerates approach to equilibrium.","The first-order scheme and its SDC acceleration are decoupled: each step solves the phase-field and concentration system, then a Stokes-type system, then a pressure correction, which makes large three-dimensional simulations practical.","The discrete energy law includes a pressure-stabilization term $\\frac{\\Delta t^2}{2Re}\\|\\nabla p\\|^2$, so the projection step contributes positively to the energy identity rather than spoiling it.","When osmotic effects are omitted, the model reduces to a degenerate Allen-Cahn-driven collapse in which droplets vanish; with solute present the droplet instead shrinks to a finite equilibrium size.","Permeability changes qualitative outcomes: in shear flow, impermeable droplets coalesce while permeable droplets shrink and avoid contact."],"supporting_citations":[{"why":"Supplies the conservation-law template and energy variational setting that the model extends to semi-permeable interfaces.","marker":"[27]"},{"why":"Provides the energy variational approach used to derive the flux and stress closures from the dissipation functional.","marker":"[12]"},{"why":"Closest prior phase-field model of mass transport with semi-permeable interfaces that the present model extends to momentum and interface dynamics.","marker":"[24]"},{"why":"Establishes sharp-interface limits for osmosis diffuse-interface approximations, the theoretical target the model aims at.","marker":"[25]"},{"why":"Phase-field vesicle growth and shrinkage model with osmotic energy that motivates the osmotic formulation used here.","marker":"[28]"},{"why":"Supplies the LDG spatial discretization framework with alternating fluxes that guarantees the adjoint-consistency property.","marker":"[30]"},{"why":"Supplies the semi-implicit spectral deferred correction method used to raise temporal accuracy beyond first order.","marker":"[17]"},{"why":"Load-bearing for the stability proofs: justifies replacing $F$ by a truncated double-well potential so that condition (16) holds.","marker":"[3]"},{"why":"Prior LDG-SDC treatment of binary fluid-surfactant systems that the numerical framework builds on.","marker":"[16]"}],"fun_headline_variants":["Energy-stable model proves droplet shape ignores permeability","Osmotic flow: energy-dissipating phase-field with LDG","Allen-Cahn flux term yields provably stable membrane flow","Droplet equilibrium set by solute, not membrane permeability"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the double-well potential has bounded second derivative (condition (16)), but the actual potential $F(\\phi)=\\frac14(\\phi^2-1)^2$ has $F''(\\phi)=3\\phi^2-1$, which is unbounded over all real $\\phi$; the paper appeals to the common practice of using a truncated double-well potential for the proofs while the numerical scheme and tests use the untruncated one.","fun_headline_variants_meta":{"raw":{"variants":["Energy-stable model proves droplet shape ignores permeability","Osmotic flow: energy-dissipating phase-field with LDG","Allen-Cahn flux term yields provably stable membrane flow","Droplet equilibrium set by solute, not membrane permeability"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1355,"prompt_tokens":920,"completion_tokens":435,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":366}},"tokens_in":536,"tokens_out":435,"duration_ms":5593,"temperature":1.0,"reasoning_tokens":366,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:11:11.293171+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the LDG scheme (30)-(32) with the untruncated potential and initial data whose phase field takes values outside $[-1,1]$, so that $|F''|$ exceeds any chosen $L$, and check whether the discrete energy $E_h^n$ is monotone non-increasing; a single time step where $E_h^{n+1} > E_h^n$ would refute the claimed unconditional stability as stated.","supporting_citations":[{"cited_title":"An energy stableC 0 finite element scheme for a phase-field model of vesicle motion and deformation.SIAM Journal on Scientific Computing, 44(1):B122–B145, 2022","cited_arxiv_id":null,"evidence_quote":"Supplies the conservation-law template and energy variational setting that the model extends to semi-permeable interfaces."},{"cited_title":"Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids.The Journal of Chemical Physics, 133(10):104104, 2010","cited_arxiv_id":null,"evidence_quote":"Provides the energy variational approach used to derive the flux and stress closures from the dissipation functional."},{"cited_title":"A phase field model for mass transport with semi-permeable interfaces.Journal of Computational Physics, 464:111334, 2022","cited_arxiv_id":null,"evidence_quote":"Closest prior phase-field model of mass transport with semi-permeable interfaces that the present model extends to momentum and interface dynamics."},{"cited_title":"Diffuse-interface approximations of osmosis free boundary problems.SIAM Journal on Applied Mathematics, 76(3):910–929, 2016","cited_arxiv_id":null,"evidence_quote":"Establishes sharp-interface limits for osmosis diffuse-interface approximations, the theoretical target the model aims at."},{"cited_title":"Phase field modeling and compu- tation of vesicle growth or shrinkage.Journal of Mathematical Biology, 86(6):97, 2023","cited_arxiv_id":null,"evidence_quote":"Phase-field vesicle growth and shrinkage model with osmotic energy that motivates the osmotic formulation used here."},{"cited_title":"LocaldiscontinuousGalerkinmethodsforhigh-ordertime-dependentpartial differential equations","cited_arxiv_id":null,"evidence_quote":"Supplies the LDG spatial discretization framework with alternating fluxes that guarantees the adjoint-consistency property."},{"cited_title":"Semi-implicit spectral deferred correction methods for highly nonlinear partial differential equations.Journal of Computational Physics, 338:269–284, 2017","cited_arxiv_id":null,"evidence_quote":"Supplies the semi-implicit spectral deferred correction method used to raise temporal accuracy beyond first order."},{"cited_title":"AnL∞ bound for solutions of the Cahn-Hilliard equation.Archive for Rational Mechanics and Analysis, 133(2):129–144, 1995","cited_arxiv_id":null,"evidence_quote":"Load-bearing for the stability proofs: justifies replacing $F$ by a truncated double-well potential so that condition (16) holds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior LDG-SDC treatment of binary fluid-surfactant systems that the numerical framework builds on."}],"review_version":1}