{"id":"5beaea70-06bc-45a7-8cae-30f42eff6ad9","arxiv_id":"2505.22026","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A single scalar PDE is derived that encodes all two-dimensional equilibrium configurations of a third-grade Korteweg fluid under the constitutive restrictions obtained from an extended Liu procedure.","lead":"This paper derives a single nonlinear elliptic equation whose solutions give all mechanical equilibrium density profiles of a two-dimensional Korteweg fluid at constant temperature under gravity. The result reduces an overdetermined pair of balance equations to one scalar PDE, with preliminary numerical solutions for a Dirichlet boundary value problem.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central reduction from (3.2) to (3.6) is asserted rather than shown, and the only numerical validation uses inconsistent Dirichlet data; an independent symbolic check is needed.","rationale":"An independent variational reading supports the scalar reduction: equation (3.6) is the Euler-Lagrange condition mu + g*y = const for the free energy corresponding to the constitutive relations (2.5), and the stress divergence satisfies grad dot T = -rho grad(mu), so solutions of (3.6) do satisfy (3.2) under the stated assumptions. Thus the mathematics of the central claim appears sound. Nevertheless, the paper as written does not show this algebra, and the numerical example in (3.11) has boundary data that are inconsistent at the corners, so the reported figures do not validate the equation. These are presentation and verification gaps rather than demonstrated mathematical errors. The reader's identified weakest assumption (separable entropy and constant temperature) is real but explicitly scoped, so it is not the most load-bearing concern. The proposal to keep a conditional verdict reflects that the central reduction needs an explicit independent check before full acceptance, but no reason has been found to reject the claim outright.","tokens_in":8344,"tokens_out":22117,"duration_ms":228287,"concrete_test":"Use a computer algebra system to substitute the constitutive relations (2.5) into (3.2) and verify that the two momentum equations are exactly the x- and y-derivatives of L = 2*rho*s1*(rho_xx+rho_yy) + d(rho*s1)/d(rho)*(rho_x^2+rho_y^2) - d(rho*s01)/d(rho) + (g/theta0)*y - kappa, up to a nonzero factor depending only on rho. Equivalently, verify that every smooth solution of L=0 satisfies both equations in (3.2), and that every smooth solution of (3.2) satisfies L=0 for some constant kappa. If the identity fails, the single-PDE reduction is incorrect; if it holds, the central claim survives and the remaining defects are the missing derivation and the inconsistent numerical boundary conditions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that every solution of the scalar PDE (3.6) identically satisfies the overdetermined equilibrium system (3.2). The paper never carries out this reduction: between (3.5) and (3.6) it says only 'after simple algebraic manipulations' and introduces an arbitrary constant kappa. Because the whole result is this equivalence, an unverified sign or an omitted regularity condition (for example, division by rho or by s1) would silently change the equilibrium set. The only numerical support is also compromised: the Dirichlet data in (3.11) assign rho0 at all four corners from the top and bottom conditions, while the left and right boundary data reach rho1 at y=d; for rho0 different from rho1 the boundary value problem has no continuous solution, so Figures 1 and 2 cannot be solutions of the stated problem. The separability assumption (3.4) and the constant-temperature condition are explicit limitations, not hidden flaws, and the non-uniform temperature case is correctly deferred. The load-bearing concern is that the central equivalence, as written, is an assertion rather than a demonstrated result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies mechanical equilibrium of a third-grade Korteweg fluid in two dimensions at constant temperature. After reviewing constitutive relations from a prior extended-Liu-procedure paper [20], it assumes that the equilibrium entropy separates as s0=s01(ρ)+s02(ε) and asserts that the overdetermined equilibrium system (3.2) reduces to the single scalar PDE (3.6) (or its dimensionless form (3.9)). This claim is the paper's central result. The paper then analyzes linear special cases, proposes a Dirichlet boundary value problem, and reports preliminary numerical solutions obtained with finite differences and a Matlab solver.","tokens_in":8579,"tokens_out":16347,"duration_ms":131304,"significance":"If the reduction is correct, the paper provides a practical tool: equilibrium density fields of Korteweg fluids satisfying Serrin's geometric condition can be computed by solving one elliptic PDE instead of an overdetermined system. Strengths include the explicit use of constitutive inputs from [20], a consistent nondimensionalization, and correct linear-case reductions. The authors are also appropriately cautious, labeling the numerics as preliminary and deferring three-dimensional and non-uniform-temperature problems. The main weaknesses are the missing derivation of the central reduction and an inconsistent numerical boundary condition.","major_comments":[{"comment":"The central reduction from the overdetermined equilibrium system (3.2) to the single scalar equation (3.6) is stated without derivation; the text between (3.5) and (3.6) says only 'after simple algebraic manipulations' and introduces an arbitrary constant κ. Since the paper's main claim is precisely that every solution of (3.6) identically satisfies (3.2), the authors should provide the complete derivation, including the role of the Serrin condition (3.3), the handling of the arbitrary function of y that arises upon integrating the x-component of (3.2), and any regularity or sign assumptions (ρ>0, s1≠0) needed for divisions. Without this, the central equivalence is an assertion rather than a demonstrated result.","section":"Section 3, Eq. (3.6)"},{"comment":"The Dirichlet data in (3.11) are discontinuous at the upper corners: at (0,d) and (1,d), the top boundary condition u(x,d)=ρ0−x^2(1−x)^2 gives u=ρ0, while the side conditions u(0,y)=u(1,y)=(ρ1−ρ0)y/d+ρ0 give u=ρ1 at y=d. Since ρ0≠ρ1, no continuous solution to the stated boundary value problem exists, so the numerical solutions in Figures 1 and 2 cannot be solutions of the stated problem. The boundary data should be made compatible (or the corner singularities should be discussed and treated) and the computations repeated.","section":"Section 3.1, Eq. (3.11) and Figures 1-2"}],"minor_comments":[{"comment":"The boundary data are written in terms of u(x,0), u(x,d), u(0,y), u(1,y), but the unknown in (3.11) is ρ; the notation should be made consistent.","section":"Section 3.1, before Eq. (3.11)"},{"comment":"The equation ρxx+ρyy+2αρ=0 is a Helmholtz equation, not a Poisson equation, since it contains the zero-order term 2αρ.","section":"Section 3.1, text after Eq. (3.10)"},{"comment":"The text says first and second derivatives are approximated by second-order and fourth-order finite differences, respectively, but the formulas displayed in (3.12) are both standard second-order central differences.","section":"Section 3.1, text before Eq. (3.12)"}],"recommendation":"major_revision","confidential_remarks":"The central reduction appears plausible and the authors are honest about the preliminary nature of the numerics. The main obstacle is the missing derivation of Eq. (3.6) and the inconsistent boundary data; both are fixable within the scope of the paper. The relation to the prior work [20] is that of a new equilibrium reduction, not a restatement, so the novelty is adequate for a journal paper if the technical gaps are closed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper takes a real step—rewriting the overdetermined 2D equilibrium system for third-grade Korteweg fluids as a single scalar elliptic PDE—and it does so modestly, with explicit limitations. The new thing is equation (3.6) and its dimensionless form (3.9); I do not see it in the cited literature. The linear reductions (m=1,n=-1 and m=-1,n=-1) check out, and the numerical solutions, while labeled preliminary, are at least plausible. The constitutive input from [20] is self-cited but used as input, not fitted; that is not circular.\n\nThe main soft spot is exactly what the stress-test says. The step from (3.2) to (3.6) is 'simple algebraic manipulations,' but the whole claim is that every solution of (3.6) satisfies the overdetermined system identically. That equivalence needs to be shown, or at least verified symbolically. There may be a division by ρ or s1 hiding in there, or a sign in the α3 term, and either would change the equilibrium set. This is load-bearing and currently unverified.\n\nSecond soft spot: the Dirichlet data in (3.11) are inconsistent when ρ0≠ρ1. The top and bottom data give ρ0 at the corners; the left and right data give ρ0 at y=0 but ρ1 at y=d, so at (0,d) and (1,d) the boundary conditions conflict. Figures 1 and 2 therefore cannot be solutions of the stated boundary value problem as written. That makes the numerics illustrative, not validating. The paper says the results are preliminary, so this is a flaw in presentation rather than a deception, but it should be fixed (e.g., re-state the boundary data with compatible corners) before anyone relies on the figures.\n\nMinor: no code, no convergence study, no residual check; the tolerance 1e-4 and grid step 0.02 are mentioned but no error analysis. For a short paper explicitly labeled preliminary, I would call this minor. The non-uniform temperature case is correctly deferred, and the separability assumption (3.4) is stated plainly.\n\nBottom line: the paper is honest, clear, and the central reduction is worth taking seriously. It deserves a referee. Whether it is fully correct depends on a symbolic check of the (3.2)→(3.6) equivalence, which a competent referee could do in an afternoon. For a reader in continuum thermodynamics or phase-boundary theory, the paper is a useful prompt; for a general math-phys reader, it is a short, self-contained look at one overdetermined system. I would not cite it yet, but I would bring it to a reading group.","headline":"Honest, compact reduction of 2D Korteweg equilibrium to one scalar PDE; the result is probably right but the key algebra is asserted and the numerics contain a boundary mismatch.","tokens_in":9119,"tokens_out":2541,"would_cite":false,"duration_ms":22306,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A10","76M20"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a two-dimensional Korteweg fluid at constant temperature, the paper shows that all mechanical equilibrium configurations are governed by a single nonlinear elliptic equation that identically satisfies the overdetermined equilibrium…","keywords":["Korteweg fluids","equilibrium configurations","overdetermined systems","nonlinear elliptic equations","third-grade fluids","capillary stress","phase boundaries","entropy principle"],"falsifier":"Numerically solve the dimensionless equilibrium equation (3.9) for a fixed parameter set on a grid, substitute the result into both components of the original equilibrium system (3.2), and evaluate the residuals; the claim that (3.6) identically satisfies the system stands or falls on whether those residuals vanish at the level of the numerical tolerance.","tokens_in":8134,"feed_emoji":"🌊","tokens_out":14032,"duration_ms":121445,"temperature":0.7,"pith_summary":"The paper investigates how a Korteweg fluid—a material whose stress depends on density and its gradients up to second order, used for capillary and phase-interface phenomena—can rest in mechanical equilibrium in a two-dimensional vertical plane under gravity. At constant temperature, the equilibrium conditions form an overdetermined system of two partial differential equations for the single unknown density field. The central claim is that, when the equilibrium entropy separates as $s_0(\\rho,\\varepsilon)=s_{01}(\\rho)+s_{02}(\\varepsilon)$, this overdetermined system is identically satisfied by any solution of the single scalar elliptic equation (3.6). The authors then specialize the free functions to power laws, rewrite the equation in dimensionless form (3.9), and show that particular parameter choices linearize it to a Poisson or Laplace equation. They present preliminary numerical solutions of a Dirichlet boundary-value problem, establishing a practical route to computing equilibrium density profiles and non-trivial phase-boundary shapes.","feed_headline":"One PDE captures 2D equilibria of Korteweg fluids","feed_subtitle":"A single scalar equation replaces an overdetermined system, making density fields and phase boundaries computable.","key_machinery":"The load-bearing object is the compatibility condition (3.3), $\\alpha_3^2-\\alpha_1\\,\\partial\\alpha_3/\\partial\\rho+2\\alpha_2\\alpha_3=0$, combined with the separated equilibrium entropy (3.4), $s_0(\\rho,\\varepsilon)=s_{01}(\\rho)+s_{02}(\\varepsilon)$. When both hold, the two equilibrium equations in (3.2) are no longer independent: they become derivatives of a single scalar expression, so the system reduces to equation (3.6). That single nonlinear elliptic equation is the mechanism carrying the argument; its dimensionless version (3.9) is the object actually solved, and the power-law specialization (3.7) turns it into a boundary-value problem amenable to finite-difference computation.","core_discovery":"The central discovery is a reduction: for a third-grade Korteweg fluid at constant temperature $\\theta_0$ in two dimensions, the overdetermined equilibrium system (3.2) collapses to the single scalar equation $2\\rho s_1(\\rho)(\\rho_{xx}+\\rho_{yy}) + \\frac{d(\\rho s_1)}{d\\rho}(\\rho_x^2+\\rho_y^2) - \\frac{d(\\rho s_{01})}{d\\rho} + \\frac{g}{\\theta_0}y - \\kappa = 0$, where $s_1(\\rho)\\le 0$ is the coefficient of $|\\nabla\\rho|^2$ in the entropy and $s_{01}(\\rho)$ is the density-dependent part of the equilibrium entropy. Any sufficiently regular solution of this equation automatically satisfies both components of (3.2). The reduction relies on the compatibility condition (3.3), which the thermodynamically derived constitutive relations satisfy precisely under the separated-entropy hypothesis (3.4). With the power-law choices $s_{01}=\\kappa_1\\rho^m$ and $s_1=-\\kappa_2\\rho^n$, the dimensionless form (3.9) is obtained; for $m=1,n=-1$ it becomes a Poisson equation, and for $m=-1,n=-1$ a Laplace equation, both admitting separable analytical solutions.","pith_inferences":["The reduction is demonstrated only for power-law choices of $s_{01}$ and $s_1$; if the compatibility condition is the only essential ingredient, the same collapse to a single equation should occur for other functional forms, which could be tested by solving (3.6) with non-power-law data and checking the residual of (3.2).","The separability assumption likely marks the boundary of the single-equation description: for non-separable equilibrium entropy, the overdetermined system should generically admit only the simple phase-boundary geometries, so a physical test would be to look for non-spherical equilibrium interfaces in a Korteweg fluid with a known non-separable entropy.","The constant-temperature hypothesis fixes $\\theta_0$; allowing $\\theta=\\theta(x,y)$ gives two unknowns and two equations, so the authors' stated future problem of non-isothermal stationary solutions is a different regime in which the single-equation shortcut may not hold."],"forward_implications":["Equilibrium density fields in two dimensions can be obtained by solving one scalar elliptic equation with boundary data instead of an overdetermined pair of PDEs.","For the parameter choices $m=1$, $n=-1$ and $m=-1$, $n=-1$, the reduced equation is linear and becomes respectively a Poisson equation and a Laplace equation, so analytical separable solutions exist.","Because the compatibility condition (3.3) is met, the classical restriction of equilibrium phase boundaries to spherical, cylindrical, or planar shapes is lifted, making more general interface geometries admissible.","The preliminary finite-difference solutions of the Dirichlet problem show that the reduced equation is numerically tractable in both linear and nonlinear regimes.","The single-equation formulation provides a starting point for three-dimensional equilibria and for comparisons with laboratory experiments, both listed by the authors as future work."],"supporting_citations":[{"why":"Supplies the thermodynamically derived constitutive relations (2.4)-(2.5) and the entropy form $s=s_0+s_1|\\nabla\\rho|^2$ that the equilibrium reduction uses.","marker":"[20]"},{"why":"Establishes the theorem on allowed phase-boundary geometries that motivates the compatibility condition (3.3).","marker":"[22]"},{"why":"Provides the overdetermined-system result combined with [22] to identify condition (3.3).","marker":"[23]"}],"fun_headline_variants":["A single PDE replaces overdetermined system for Korteweg fluids","2D Korteweg equilibria reduced to one elliptic equation","Korteweg fluids: equilibrium solved by a single scalar equation","From overdetermined to one: Korteweg equilibrium in 2D"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation requires that the equilibrium entropy split as $s_0(\\rho,\\varepsilon)=s_{01}(\\rho)+s_{02}(\\varepsilon)$ and that the equilibrium temperature be constant; if a real fluid has a non-separable equilibrium entropy or a non-uniform temperature field, equation (3.6) does not describe its equilibrium density.","fun_headline_variants_meta":{"raw":{"variants":["A single PDE replaces overdetermined system for Korteweg fluids","2D Korteweg equilibria reduced to one elliptic equation","Korteweg fluids: equilibrium solved by a single scalar equation","From overdetermined to one: Korteweg equilibrium in 2D"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000743,"raw_usage":{"total_tokens":3286,"prompt_tokens":890,"completion_tokens":2396,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":2317}},"tokens_in":506,"tokens_out":2396,"duration_ms":18274,"temperature":1.0,"reasoning_tokens":2317,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:16:58.295087+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically solve the dimensionless equilibrium equation (3.9) for a fixed parameter set on a grid, substitute the result into both components of the original equilibrium system (3.2), and evaluate the residuals; the claim that (3.6) identically satisfies the system stands or falls on whether those residuals vanish at the level of the numerical tolerance.","supporting_citations":[{"cited_title":"Gorgone, P","cited_arxiv_id":null,"evidence_quote":"Supplies the thermodynamically derived constitutive relations (2.4)-(2.5) and the entropy form $s=s_0+s_1|\\nabla\\rho|^2$ that the equilibrium reduction uses."},{"cited_title":"Serrin, The form of interfacial surfaces in Korteweg’s theory of pha se equilibria, Quart","cited_arxiv_id":null,"evidence_quote":"Establishes the theorem on allowed phase-boundary geometries that motivates the compatibility condition (3.3)."},{"cited_title":"Pucci, An overdetermined system , Quart","cited_arxiv_id":null,"evidence_quote":"Provides the overdetermined-system result combined with [22] to identify condition (3.3)."}],"review_version":1}