Pith. sign in

REVIEW 2 major objections 3 minor 24 references

Two-dimensional equilibrium configurations in Korteweg fluids

T0 review · 2 major / 3 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read 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…

desk verdict 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. read the letter →

arxiv 2505.22026 v1 pith:NALLJ4NZ submitted 2025-05-28 math-ph math.MP

classification math-phmath.MP MSC 76A1076M20
keywords Kortewegfluidsequilibriumconfigurationsoverdeterminedsystemsnonlinearellipticequationsthird-gradecapillarystressphaseboundariesentropyprinciple
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

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.

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 (2)
  1. [Section 3, Eq. (3.6)] 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.
  2. [Section 3.1, Eq. (3.11) and Figures 1-2] 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.
minor comments (3)
  1. [Section 3.1, before Eq. (3.11)] 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.
  2. [Section 3.1, text after Eq. (3.10)] The equation ρxx+ρyy+2αρ=0 is a Helmholtz equation, not a Poisson equation, since it contains the zero-order term 2αρ.
  3. [Section 3.1, text before Eq. (3.12)] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the scalar equilibrium equation is a genuine reduction, not a disguised restatement of the paper's inputs.

full rationale

The central claim—that every solution of the scalar PDE (3.6) identically satisfies the overdetermined equilibrium system (3.2)—is a mathematical reduction, not a definition or a fit. The paper combines the equilibrium condition (3.1)-(3.2) with the constitutive relations (2.5) inherited from [20] and the separability condition (3.4), and it introduces the constant κ as an integration constant. No equation is defined in terms of the target result, and no parameter is fitted to data and then renamed as a prediction: the forms s01=κ1ρ^m and s1=-κ2ρ^n in (3.7) are free choices, and the numerical solutions are explicitly labeled preliminary. The self-citation to [20] is essential input, but [20] derives the constitutive functions from an extended Liu procedure that does not include the equilibrium PDE (3.6); under the review rules this is independent support, not circularity. The Serrin/Pucci theorems are external, not self-citations. The paper's acknowledged limitations—constant temperature θ0, separable s0, and deferral of non-uniform temperature to future work—are scope restrictions, not circular reasoning. Two non-circular weaknesses are flagged: the phrase 'after simple algebraic manipulations' omits the details of the reduction to (3.6), and the Dirichlet data in (3.11) are inconsistent at the corners when ρ0≠ρ1; both are correctness/transparency concerns, not circularity.

Assumptions & free parameters 9 free parameters · 4 assumptions · 0 invented entities

The central equation (3.6) rests on the constitutive functions from [20] (an extended Liu procedure), the separability of the equilibrium entropy, and constant temperature. The power-law forms and the dimensionless parameters α, β, γ are chosen by hand for the numerical demonstrations, not fitted to data. No new physical entities are introduced.

free parameters (9)
  • m = 1 or -1 in linear cases; 1 in nonlinear examples
    Exponent in the chosen power-law form s01 = κ1 ρ^m; equation (3.9) depends on m.
  • n = -1 in linear cases; 0, 1, -2, -3 in nonlinear examples
    Exponent in the chosen power-law form s1 = -κ2 ρ^n; equation (3.9) depends on n.
  • κ1 = not numerically specified
    Positive constant in s01 = κ1 ρ^m; absorbed into dimensionless parameter α.
  • κ2 = not numerically specified
    Positive constant in s1 = -κ2 ρ^n; absorbed into dimensionless parameters α, β, γ.
  • α = 1.0 in all numerical plots
    Dimensionless parameter from κ1/(2κ2) ℓ1^2 R0^{m-n-2}; chosen by hand for numerical experiments.
  • β = -1.2 in all numerical plots
    Dimensionless parameter from -g/(2κ2θ0) ℓ1^3 R0^{-n-2}; chosen by hand.
  • γ = ±1 in numerical plots
    Dimensionless parameter from κ/(2κ2) ℓ1^2 R0^{-n-2}; chosen by hand.
  • ρ0, ρ1 = ρ0=1.4, ρ1=1.3 in numerical plots
    Boundary density constants in the Dirichlet conditions (3.11).
  • d = 3 in numerical plots
    Aspect ratio of the domain [0,1]×[0,d].
assumptions (4)
  • domain assumption Constitutive relations (2.5) from [20] are accepted as correct.
    The paper reviews but does not re-derive the extended Liu procedure results; the equilibrium reduction is built on these constitutive expressions.
  • standard math Serrin's compatibility condition (3.3) is needed for non-simple phase boundaries.
    The theorem from [22,23] is invoked in Section 3 to justify the compatibility condition.
  • ad hoc to paper The equilibrium entropy separates as s0=s01(ρ)+s02(ε).
    Equation (3.4) is introduced in Section 3 to simplify the reduction; it restricts the class of Korteweg fluids considered.
  • domain assumption Temperature is constant at equilibrium (θ0).
    Used throughout Section 3 to evaluate material functions at constant temperature and to set ε=ε(θ).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Two-dimensional equilibrium configurations in Korteweg fluids." pith.science (2026). https://pith.science/paper/NALLJ4NZ

@misc{pith2026250522026,
  author       = {Pith},
  title        = {Pith review of: Two-dimensional equilibrium configurations in Korteweg fluids},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NALLJ4NZ}},
  note         = {Machine review of arXiv:2505.22026}
}
read the original abstract

In this paper, after reviewing the form of the constitutive equations for a third grade Korteweg fluid, recently derived by means of an extended Liu procedure, an equilibrium problem is investigated. By considering a two--dimensional setting, it is derived a single nonlinear elliptic equation such that the equilibrium conditions are identically satisfied. Such an equation is discussed both analytically and numerically. Moreover, by considering a particular boundary value problem of Dirichlet type, some preliminary numerical solutions are presented.

Figures

Figures reproduced from arXiv: 2505.22026 by the authors.

Figure 1
Figure 1. Plot of the density ρ (left) and contour plot (right). The values of the parameters are (from the top): (m = 1, γ = 1), (m = 1, γ = −1), (m = −1, γ = 1), (m = −1, γ = −1). 9 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Plot of the density ρ (left) and contour plot (right). The values of the parameters are (from the top): (m = 1, n = −2), (m = 1, n = −3), (m = 1, n = 1), (m = 1, n = 0). 10 [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages

  1. [20]

    Gorgone, P

    M. Gorgone, P. Rogolino, On the characterization of constitutive equa- tions for third-grade viscous Korteweg fluids , Phys. Fluids 33 (2021), 043107

  2. [1]

    D. J. Korteweg, Sur la forme qui prennent les ´ equations du mouvement des fluids si l’on tient compte des forces capillaires par des variations de densit´ e consid´ erables mais continues et sur la th´ eorie de la capil- larit´ e dans l’hypoth` ese d’une variation continue de la densit´ e, Archives N´ eerlandaises des sciences exactes et naturelles 6 (190...

  3. [2]

    Heida, J

    M. Heida, J. M´ alek, On Korteweg–type compressible fluid-like materials , Int. J. Eng. Sci. 48 (2010), 1313–1324

  4. [3]

    J. E. Dunn, K. R. Rajagopal, Fluids of differential type: critical review and thermodynamic analysis , Int. J. Eng. Sci. 33 (1995), 689–729

  5. [4]

    Truesdell, K

    C. Truesdell, K. R. Rajagopal, An Introduction to the Mechanics of Fluids, Birkh¨ auser, Boston–Basel–Berlin, 2000

  6. [5]

    J. E. Dunn, J. Serrin, On the thermomechanics of the interstitial work- ing, Arch. Rat. Mech. Anal. 88 (1985), 95–133

  7. [6]

    J. E. Dunn, Interstitial working and a nonclassical continuum thermo- dynamics, in: J. Serrin (ed.), New Perspectives in Thermodynamics , Springer Verlag, Berlin, 1986, 11, 187–222

  8. [7]

    Truesdell, Rational Thermodynamics

    C. Truesdell, Rational Thermodynamics. 2nd ed. , Springer Verlag, Berlin, 1984

Show all 24 references
  1. [8]

    M¨ uller,On the Entropy Inequality , Arch

    I. M¨ uller,On the Entropy Inequality , Arch. Rat. Mech. Anal. 26 (1967), 118–141

  2. [9]

    V. A. Cimmelli, A. Sellitto, V. Triani, A new thermodynamic framework for second–grade Korteweg–type viscous fluids, J. Math. Phys. 50 (2009), 053101

  3. [10]

    V. A. Cimmelli, A. Sellitto, V. Triani, A new perspective on the form of the first and second laws in rational thermodynamics: Kortew eg fluids as an example , J. Non–Equil. Thermodyn. 35 (2010), 251–265. 13

  4. [11]

    V. A. Cimmelli, F. Oliveri, A. R. Pace, On the thermodynamics of Ko- rteweg fluids with heat conduction and viscosity , J. Elast. 104 (2011), 115–131

  5. [12]

    V. A. Cimmelli, F. Oliveri, A. R. Pace, Thermodynamical setting for gradient continuum theories with vectorial internal varia bles: application to granular materials , Int. J. Non–Linear Mech. 49 (2013), 72–76

  6. [13]

    V. A. Cimmelli, F. Oliveri, A. R. Pace, A nonlocal phase–field model of Ginzburg–Landau–Korteweg fluids , Contin. Mech. Thermodyn. 27 (2015), 367–378

  7. [14]

    Oliveri, A

    F. Oliveri, A. Palumbo, P. Rogolino, On a model of mixtures with inter- nal variables: extended Liu procedure for the exploitation of the entropy principle, Atti Accad. Pelorit. Pericol. Cl. Sci. Fis. Mat. Nat. 94 (2016), A2–1–17

  8. [15]

    V. A. Cimmelli, M. Gorgone, F. Oliveri, A. R. Pace, Weakly nonlocal thermodynamics of binary mixtures of Korteweg fluids with tw o velocities and two temperatures , Eur. J. Mech. B/Fluids 83 (2020), 58–65

  9. [16]

    Gorgone, F

    M. Gorgone, F. Oliveri, P. Rogolino, Continua with non–local constitu- tive laws: Exploitation of entropy inequality , Int. J. Non–Linear Mech. 126 (2020), 103573

  10. [17]

    Gorgone, F

    M. Gorgone, F. Oliveri, P. Rogolino, Thermodynamical analysis and constitutive equations for a mixture of viscous Korteweg flu ids, Phys. Fluids 33 (2021), 093102

  11. [18]

    Gorgone, P

    M. Gorgone, P. Rogolino, A thermodynamical description of third grade fluid mixtures , J. Non–Equil. Thermodyn. 47 (2022), 133–142

  12. [19]

    I–Shih Liu, Method of Lagrange multipliers for exploitation of the en- tropy principle , Arch. Rat. Mech. Anal. 46 (1972), 131–148

  13. [21]

    V. A. Cimmelli, An extension of Liu procedure in weakly nonlocal ther- modynamics, J. Math. Phys. 48 (2007), 113510. 14

  14. [22]

    Serrin, The form of interfacial surfaces in Korteweg’s theory of pha se equilibria, Quart

    J. Serrin, The form of interfacial surfaces in Korteweg’s theory of pha se equilibria, Quart. Appl. Math. 41 (1983), 357–364

  15. [23]

    Pucci, An overdetermined system , Quart

    P. Pucci, An overdetermined system , Quart. Appl. Math. 41 (1983), 365–367

  16. [24]

    A. R. Mitchell, D. F. Griffiths, The Finite Difference Method in Partial Differential Equations , John Wiley & Sons, New York, 1980. 15

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.